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Abstract 

We review results on hadron multiplicities in high energy particle col- 
lisions. Both theory and experiment are discussed. The general pro- 
cedures used to describe particle multiplicity in Quantum Chromo- 
dynamics (QCD) are summarized. The QCD equations for the gen- 
erating functions of the multiplicity distributions are presented both 
for fixed and running coupling strengths. The mean multiplicities of 
gluon and quark jets, their ratio, higher moments, and the slopes of 
multiplicities as a function of energy scale, are among the main global 
features of multiplicity for which QCD results exist. Recent data from 
high energy e"'"e~ experiments, including results for separated quark 
and gluon jets, allow rather direct tests of these results. The theo- 
retical predictions are generally quite successful when confronted with 
data. Jet and subjet multiplicities are described. Multiplicity in lim- 
ited regions of phase space is discussed in the context of intermittency 
and fractality. The problem of singularities in the generating functions 
is formulated. Some special features of average multiplicities in heavy 
quark jets are described. 
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1 Introduction 



The number of hadrons created in high energy inelastic coUisions varies from 
one event to another. The distribution of the number of produced hadrons 
(the multiphcity distribution, for short) provides a basic means to charac- 
terize the events. The multiphcity distribution contains information about 
multiparticle correlations in an integrated form, providing a general and sen- 
sitive means to probe the dynamics of the interaction. 

Here, we would like to review the main theoretical approaches and exper- 
imental results concerning hadron multiplicities. Much precise experimental 
information on multiplicities has become available in recent years from e"'"e~, 
hadron, and nucleus collisions. The experimental progress has stimulated 
additional theoretical developments. Therefore, the time seems appropriate 
for a review. The theoretical results summarized here represent an expanded 
and updated version of an earlier work 

After introducing basic definitions and notation in Section 0, we turn 
in Section ^ to describe some phenomenological approaches to multiplicity. 
The phenomenological methods are based on simplified ideas about particle 
emission, exploiting widely used distributions from probability theory (see 
e.g. Ref. [^). Among them, the negative binomial distribution is especially 
popular because it describes experimental data reasonably well for various 
reactions over a wide energy interval. 

The discovery of asymptotic freedom, in conjunction with the parton 
model used to characterize deep inelastic and e^e~ data, led to the de- 
velopment of the theory of strong interactions. Quantum Chromodynamics 
(QCD). QCD provides a means to apply perturbative techniques to hadronic 
processes with large transferred momenta, with quarks and gluons playing 
the role of partons 0-0. The current state of affairs concerning multiplicity 
in quark and gluon jets as predicted by QCD and measured by experiment 
forms the bulk of our review, presented in Sections ^|15[ 

Lattice calculations imply that hadronization - the transition of quarks 
and gluons to hadrons at the final stage of evolution - is an inherent prop- 
erty of QCD. Nonetheless, we are still unable to treat hadronization in an 
unambiguous manner since the problem of confinement has yet to be solved 
mathematically. Simplified estimates suggest either that hadronization does 
not drastically alter the parton level results or else that its effects can be 
estimated from the energy dependence of experimental observables. Phe- 
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nomenologically, the distributions of partons and hadrons are often found to 
be remarkably similar. This implies that the study of the partonic stage of 
an event is of utmost importance since the properties of high energy multi- 
hadronic events are primarily determined at that level. 

The probabilistic equations of QCD for the parton cascade can be solved 
within the perturbative approach. The parton shower framework allows the 
use of sub-series of the perturbative expansion, with terms ordered according 
to their energy behavior. The zero order approximation already involves an 
arbitrary number of produced particles. Higher order approximations lead to 
detailed predictions for multiplicity distributions which can be tested through 
comparison with experiment. Before proceeding to details in the following 
sections, however, we wish to note two problematic areas and then briefly 
comment on the main topics we address. 

First, QCD yields results on partons, and not hadrons, as already men- 
tioned. Therefore an assumption about hadronization must be made in 
any comparison of experiment with theory. A common assumption is local 
parton-hadron duality (LPHD) p], which states that parton distributions are 
simply renormalized in the process of hadronization, without changing their 
shape. LPHD originated from the idea of soft preconfinement [Q, whereby 
partons group in colorless clusters without disturbing the initial spectra. Phe- 
nomenological models of hadronization have been incorporated into Monte 
Carlo simulations of inelastic processes and in most cases support the ap- 
proximate property of LPHD. 

Second, perturbative analysis has limitations in the context of any specific 
approximation. Formally speaking, perturbation theory should be used only 
when the coupling constant is very small compared to unity. This condition 
is fulfilled by QCD only for extremely large transferred momenta, however, 
often for energies much larger than present experimental conditions. In any 
QCD process, the energies of cascading partons degrade during their evolu- 
tion, and a proper accounting for soft partons, their recoil due to interaction, 
and energy-momentum conservation laws, should be included. All these con- 
siderations are neglected in the lowest order approximation, for which only 
processes with a rather large gradient of energies and emission angles at each 
stage of the evolution are considered (the so-called double-logarithmic ap- 
proximation, DLAQ). Account is taken of soft partons and strict transverse 

^Also referred to as leading order (LO). 
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momentum ordering in subsequent terms of the perturbative series, such as 
the modified leading- logarithmic approximation (MLLA[|). Recoil effects and 
conservation laws can be incorporated at the next-to-next-to-leading order 
(NNLO), the next-to-next-to-next-to- leading order (3NL0), and higher or- 
der stages. The conservation laws are more accurately included, the higher 
the order of the perturbative approximation. 

In most cases these corrections are under control, being about 10-20% 
percent of the main terms at present energies. In spite of their rather small 
total contribution, they are often very important and change the qualita- 
tive theoretical description in regions where the corresponding functions are 
small. For example, they are crucial for the proper description of the mul- 
tiparticle production process. This manifests itself mathematically as a new 
expansion parameter equal to the product of the QCD coupling constant (or, 
more precisely, of its square root) and the rank of the moment of the distri- 
bution. Thus it is for large ranks, i.e. high multiplicities, that the corrections 
are most important. These problems are discussed in detail in Sections 

The first results on the multiplicity distributions of partons were obtained 
using the double- logarithmic approximation (for reviews, see Refs. 0, |10|). 
They yielded asymptotic Koba-Nielsen-Olesen (KNO) scaling Accord- 
ing to the KNO hypothesis |]12|, the multiplicity distribution depends only on 
the ratio of the number of particles to the average multiplicity (see Section |^). 
KNO scaling failed to be valid in the higher order approximations, however, 
i.e. the asymptotic realm is too far from present energies. At the same time, 
the width of the asymptotic KNO function was found to be much larger 
than experimental distributions. This problem was resolved [0, Q through 
consideration of the higher order effects mentioned above. The increase of av- 
erage multiplicity with energy was found to depend on the coupling constant 
in a manner which is faster than any logarithmic function and slower than 
any power-like one (if the running coupling constant is used), in agreement 
with experimental findings. One can now state that an overall agreement 
between theory and data has been achieved, at least qualitatively. 

Moreover, some qualitative predictions of perturbative QCD (pQCD) 
were found to be unexpectedly well suited for "soft" hadronic processes as 
well. Given the large value of the expansion parameter at small scales, this 
may seem puzzling. However, higher order corrections should account for 



^Also referred to as next-to-leading-order (NLO). 
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ever softer partons in a consistent manner, implying a more general cor- 
respondence between soft and hard processes than is usually considered in 
theoretical schemes. For example, the newly discovered feature of the oscil- 
lation of cumulant moments (see Section |5l^ ), prompted by solutions of the 
QCD equations, were experimentally observed in e"'"e~, hadron-, and nucleus- 
induced reactions. These oscillations were found to be extremely sensitive to 
small details of the multiplicity distributions. 

Many discussions are devoted to the value of the ratio of average multi- 
plicities between gluon and quark jets, r. The value obtained in the double- 
logarithmic approximation is r=9/4. The simplest corrections reduce this 
value by about 10%. An even larger decrease of r is predicted by the exact 
solution of the equations for the generating functions in the case of a fixed 
coupling constant, by the analytic higher order approximations with a run- 
ning coupling constant, and by computer solutions of these equations. These 
results are discussed in Sections A proper accounting for the boundary 
between the perturbative and non-perturbative regions appears to be cru- 
cial for this quantity. The energy dependence of the average multiplicities, 
i.e. their slopes, is considered in Section |[ The ratios of the slopes and of 
the corresponding curvatures are found to be less sensitive to higher order 
corrections than the ratio of multiplicities itself. Some attempts to account 
for non-perturbative effects in evolving jets are briefly described in Section ^ 

The behavior of the moments of multiplicity is strongly influenced by 
the nature of the singularities of the generating function, which are not yet 



known. Approaches to this problem are discussed in Section 10 



Experimental measurements of multiplicity as they relate to tests of the 



QCD predictions are described in Section |Tl|. Especially noteworthy are 
recent data from e"'"e~ annihilation experiments operating on the Z° peak, 
namely the four experiments ALEPH, DELPHI, L3 and OPAL at the LEP 
storage ring at CERN and the SLD experiment at the SLC collider at SLAG. 
The simplicity of the hadronic production process in e"'"e~ events, along with 
the large data samples collected by these experiments, has resulted in an 
unprecedented level of experimental precision for multiplicity related quanti- 
ties, leading to the possibility, for example, of examining multiplicity in small 
phase space windows. Recent progress in distinguishing quark and gluon jets 
by these experiments has made it possible to study their properties sepa- 
rately. In particular, the multiplicity distributions for each set have recently 
been analyzed for the first time p 
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The evolution of jets can be studied by resolving subjets, and by de- 
termining the subjet multiplicity rate as a function of the resolution scale. 
Theoretical and experiment results on subjet multiplicities are described in 
Section |1^. This is followed in Section |T3| by a brief discussion of particle 
multiplicity in e"'"e~ three-jet events. 

It is of interest to study multiplicity not only in total phase space but 
also in small subregions. These studies are generally focused on the inter- 
mittency phenomenon and on the fractality of particle distributions within 
a selected phase space volume (for a recent review, see ||16|), related to a 
relative widening of the multiplicity distribution for smaller phase space vol- 
umes. Intermittency gives rise to an increase of multiplicity moments in a 
power-like manner as the phase space window decreases. Such tendencies 
have been experimentally observed. Quantum Chromo dynamics provides 
a qualitative description of the increase of the moments, relates the inter- 
mittency exponents (or fractal dimensions) directly to the QCD anomalous 
dimension (i.e. the coupling constant), and clearly delineates the region of 
applicability of the regularities, indicating the scales at which one should 
consider the coupling constant to be running or at which it can be treated as 
approximately fixed. Theoretical and experimental aspects of intermittency 
and fractality are described in Section 

A deeper understanding of specific features of multiplicity can be gained 
if inclusive distributions of particles and their mutual correlations are con- 
sidered. The quantum mechanical origin of the interacting partons reveals 
itself in various interference effects. They lead to the hump-backed plateau 
of rapidity distributions, to correlations of partons in energies and azimuthal 
angles, to the string (or drag) effect in three-jet events, and to interference 
phenomena in the production of heavy bosons and lepton pairs at large trans- 
verse momenta. We do not describe these results here, instead referring the 
reader to monographs 0, |^ and to a recent review |17|. Nonetheless, we 
discuss one interference effect, namely the suppression of the forward pro- 
duction of accompanying particles in processes with heavy quarks. Inclusion 
of this topic is justified here because it directly affects the relation between 
the mean multiplicities in heavy- and light-quark jets. A summary of theo- 
retical and experimental results on multiplicity in heavy quark jets is given 
in Section |15[ 

To keep this review to a reasonable length, we do not describe the in- 
teractions of hadrons, nuclei or polarized quarks except in passing. These 
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topics merit their own review. Our main purpose is to present a coherent, 
updated overview of the theoretical and experimental status of multiplicity 
in high energy hadron jets. We apologize to the authors of papers whose 
contributions have not been mentioned. These omissions are unintentional. 



2 Definitions and notation 

The multiplicity distribution is defined by the formula 



Eoo ■ 
n=0 n '^inel 



(1) 



where (Tn is the cross section of an n-particle production process (the so-called 
topological cross section), Uinei is the inelastic cross section, and the sum is 
over all possible values of n so that 

oo 

E^n = l. (2) 

n=0 

It is often more convenient to represent the multiplicity distribution by its 
moments, i.e. by another set of numbers obtained from it by a definite algo- 
rithm. All such sets can be obtained from the so-called generating function 
defined by the formula 

oo 

G{y,z) = Y.PM{l + z)\ (3) 

n=0 

which substitutes an analytic function of z in place of the set of numbers 
Pn{y) at a fixed energy y. 

In what follows, use will often be made of the (normalized) factorial 
moments Fq and cumulants Kq determined from the generating function 
G{z) by the relations 

Y.nPnn{n-l)...{n-q + l) 1 d'iG{z) 



{Y^nPnUY {nY dzi 

1 dniiG{z) 



K — 
^ dzi 



z=o, (4) 
z=o, (5) 
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where 



n 



n 



n=0 



is the average multiphcity. The expression for G{z) can be written as 



g=0 ^• 



oo 



g=l 



(6) 



(7) 



(8) 



The distribution P„ and its ordinary moments Cq can be derived from the 
generating function G{z) using the formulas 



1_ d-G{z] 

rj 



a, 



1 



2=0- 



(9) 



(10) 



{n)i dzi 

All the moments are connected by definite relations that can be derived 
from their definitions given above. For example, the factorial moments and 
cumulants are related to each other through the identities 



9-1 



P \ ^ rim TV- TP 

/ , ^ q—l-'^^q—m-'^ mi 



(11) 



m=0 



which are nothing other than the relations between the derivatives of a func- 
tion and its logarithm at the point where the function itself equals unity. 
Here 



r(g) 



1 



m!(g — m — 1)! T {m + 1)T {q — m) mB{q,m) 



;i2) 



are the binomial coefficients, and F and B denote the gamma and beta func- 
tions, respectively. There are only numerical coefficients in the recurrence 
relations ([TT|) . Therefore iterative solution yields all cumulants if the factorial 
moments are given, and vice versa. In this sense, cumulants and factorial 
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moments are equally suitable. The relations between them for the lowest 
ranks are 

Fi = Ki = l, 
F2 = K2 + I, 
F3 = KS + 3K2 + I, 

F5 = 7^5 + 57^4 + 10/^3^2 + 10/^3 + 15/^2+107^2 + 1. (13) 

The physical meaning of these moments can be seen from their definitions 
if they are presented in the form of integrals of correlation functions (for a re- 
view, see Ref. []TB[)- Let the single symbol y represent all kinematic variables 
needed to specify the position of each particle in the phase space volume Q. 
A sequence of inclusive g-particle differential cross sections d'^a / dyi . . . dyq 
defines the factorial moments in the following manner: 

If r d'^a 

^1 = / dyi... dyg- — . (14) 

(Tinei{n)'' Jn Jn dyi . . . dyg 

Therefore, the factorial moments include in integrated form all correlations 
within the system of particles under consideration. The corresponding ex- 
pressions for cumulants are more complicated as can be seen from rela- 
tions (p!T| ) and (|13]). Therefore we do not present them here but instead 
remark that analogous relations exist in quantum field theory where formu- 
las similar to eqs. (^ and (^ define the complete set of Feynman graphs 
(both connected and disconnected) and the subset of connected diagrams, 
respectively (see e.g. Ref. |^). Thus, factorial moments represent integral 
characteristics of any correlation between the particles whereas cumulants of 
rank q represent genuine g-particle correlations not reducible to the product 
of lower order correlations.^ To be more precise, all q particles are connected 
to each other in the gth cumulant and cannot be split into disconnected 
groups. One can say they form a g-particle cluster which is not divisible into 
smaller clusters, in analogy with Mayer cluster decomposition in statistical 
mechanics. 

■^This interpretation is valid, however, only for moments with a rank smaller than the 
average multiplicity at a given energy (for more details, see Ref. ]l6|). 
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It is a common feature of distributions in particle physics that their facto- 
rial moments and cumulants increase rapidly as the rank q increases. There- 
fore it is often convenient to consider their ratio introduced in Ref. IM: 



= ^, (15) 

which does not share this feature. At the same time the Hg ratios exhibit 
all the qualitative features of the cumulants. As will be discussed below, the 
ratios Hq appear in a natural manner as the solution of the QCD equations 
for the generating functions \T^\. 



From eq. (^, one can derive the relation between the factorial moments 
Fq and the ordinary moments Cq with the same or lower ranks. The two 
types of moments are found to differ by terms which depend only on the 
mean multiplicity so that, for example, 

F^--^^^-C._-(nrK (16) 

Thus ordinary and factorial moments have a different energy dependence, 
since {n)~^ is energy dependent. In general, factorial moments differ from 
ordinary moments by lower order correlation terms suppressed by the inverse 
power of the mean multiplicity to the corresponding power. Thus ordinary 
and factorial moments coincide asymptotically. 

The generating function contains the same physical information as the 
multiplicity distribution. This information is also contained in the unnor- 
malized moments and their ratio. For normalized moments, the average 
multiplicity at a given energy also needs to be specified. Note that higher 
rank moments emphasize higher multiplicity events, i.e. multiplicities n > q 
contribute to factorial moments of (integer) rank q, as seen from eq. (^). If 
the distribution is truncated at some n = nmax, all factorial moments with 
rank q > rimax are zero, whereas they are positive for smaller q. The cumu- 
lants may be either positive or negative if the distribution is truncated. 

An emphasis on low multiplicity events can be made by considering 
derivatives of the generating functions at z=-l. For example, the so-called 
combinants [1^, |T^ , defined by 



_ 1 rinGjz) 

- d dz- ^^^^ 
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measure correlations within the low multiplicity part of the distribution. 

Of potential interest would be moments with an exponentially damped 
high-multiplicity tail, defined by 



^ 1 d^G{e^-r 
{n)i dzi 



z=-l 



(18) 



(see the expression for the ordinary moments Cq given by eq. (0)). Events 
with multiplicity n =q are enhanced in these moments. This type of moment 
has not yet been studied, however. 

To this point we have assumed that the rank of a moment is a positive 



integer. However, the definitions (^, (|^) and (10) can be easily generalized 



to non-integer moments pO|, by re-expressing the factorial moments as 



which is valid for any real value of q. This result can also be obtained using 
a more general approach, by forming the derivative of any real order q of the 
generating function, i.e. by using fractional calculus. Fractional moments 



have not yet found much application in particle physics (see section p. 3. 2 
however). Therefore, we refer the reader to Refs. [|l|, [|2l|-[|23| for further 
discussion. 

The method of generating functions is a particular realization of the more 
general method of generating functionals. The latter approach considers the 
probability to find a certain n-particle distribution inside the phase space vol- 
ume, in place of P„ in eq. (|), and therefore utilizes functions zijji) instead of 
the variable z, where the z{yi) depend on the phase space variables Ui. The 
corresponding distributions and correlation functions are obtained by taking 
the variational derivatives of the generating functional over the probing func- 
tions z{yi). The moments of the multiplicity distribution are related to the 
integrals of the correlation functions, as seen from eq. (iMl) . For more details, 



see Refs. [16, 17] 



3 Phenomenology 
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3.1 KNO scaling and F scaling 

A principal phenomenological issue is the energy dependence of multiplicity 
for different colliding particles and nuclei. One of the most successful as- 
sumptions about the shape of the multiplicity distribution at high energies 
is the hypothesis that its energy dependence is determined by the average 
multiplicity in such a manner that P„ may be represented as: 

1 Tl 

Pn = T-yf^T-O- (20) 

(n) {n) 



This property is called KNO scaling after the names of its authors ||T2[, who 
proposed it on the basis of the Feynman plateau of rapidity distributions. 
Earlier, a similar relation was obtained in the framework of conformal field 



theories KNO scaling is usually considered to be an asymptotic property, 
i.e. valid in the limit {n) —>■ oo. The multiplicity distribution is restored from 
the asymptotic function f{x) through the Poisson transform 



Pin) 



' f^^^M^^-{n).^^_ (21) 



n\ 



The function f{x) is unspecified. 
The normalization conditions are 



poo roo 

/ f{x)dx = 1, / xf{x)dx = 1. (22) 
Jo Jo 



It is clear that the ordinary moments of the KNO distribution, eq. (poD, do 
not depend on energy and are just functions of their rank q: 

roo 

Cq = x'^f{x)dx = independent of E. (23) 
Jo 

In contrast, the factorial moments of the KNO distribution do depend on 
energy, as follows from eq. (|l^) for F2, tending to constant values only at 
asymptotically high energies. Constancy of the factorial moments is called 
F scaling. Therefore F scaling coincides with KNO scaling asymptotically. 
The generating function depends on the average multiplicity in both cases 
as seen from the definitions (|^), (0) and (H). 

If Feynman scaling is approximate, KNO scaling is violated, suggesting 



the replacement of eq. ( pOD by other functions. In particular, it has been 
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proposed to utilize the set of relations with arbitrary functions /; 



_ k 

which reduces to KNO result for k=0. 

Another proposal is based on the scaling properties of systems undergoing 
a second order phase transition, and yields [|^, ^ 

which reduces to KNO case for 5=1. It has been argued that solutions of this 
type appear for QCD equations when they are modified by non-perturbative 
terms (see Section ||). It has also been proposed to consider so-called log- 
KNO scaling pSf, which has the general form 

(26) 




which amounts to a scale- and power-transform of multiplicity. 

In eq. (|20|) , a continuous function / is used to approximate the discrete 
distribution P„. This procedure is justified for {n) ^ 1. The procedure to 
restore the discrete distribution P„ from the continuous function / is 

drn J m_\ 
. {my \{m) )' ^ ^ 

Note that (m) is slightly shifted compared to (n), i.e. (m) ~ (n) + 0.5. 



3.2 Entropy scaling 

A convenient means to study the energy dependence of multiplicity in an 
integrated form is provided by the so-called information entropy |20|, defined 
by 

oo 

^=-^P„lnP„. (28) 

n=0 

Information entropy possesses some peculiar properties (see Refs. 
viz.: 
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The entropy of k statistically independent sources equals the sum of 
the entropies of the individual sources: 5" = Y.m=i ^m- 

The entropy is invariant under an arbitrary distortion of the multiplic- 
ity scale, in particular if a subsample of particles is chosen such as the 
charged particles. 

There exists a relationship between the entropy S, the average mul- 
tiplicity (n), and the KNO function /(x), as long as the KNO form 
eq. (^) is valid: 



S - \n{n) = / fix) In ---dx. (29) 
Jo fix) 



fix) 

• For functions satisfying the normalization conditions p^), the right 
hand side of eq. (^) does not exceed unity. Therefore the following 
inequality is valid: 

m{n) m(n) 

From eq. ( POD we see that KNO scaling is equivalent to scaling of the ratio 
S/ ln(n) in the limit of asymptotics, and that the approach to this limit is 
very slow, of order (9(l/ln(n)). In contrast, the ratio of the entropy 5* to 
the total available rapidity range, namely S/ \r^{^/s/m■,^), has been shown to 
exhibit scaling at experimentally accessible scales. The observed constancy 



of S / \n{^/s / m,,) is called entropy scaling If correct, entropy scaling 

implies a violation of KNO scaling unless the average multiplicity increases 
as a power of energy. 

If subsamples of particles are chosen, it is necessary to employ the proper 
mean values of the subset multiplicities and the correspondingly normalized 
KNO functions in eqs. (|29|) and ( pOD (for more details, see [0). The values 
of entropy are the same for the subsamples as for the full sample, as noted 
above. 



The above treatment can be generalized using the so-called Renyi 
information entropy of gth order: 

1 oo 

^. = T^1^E(W (31) 

y n=0 
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Entropy S is obtained from eq. ( |3lD by setting q=l. Such a generalization is 
reminiscent of the moments approach to multiphcity described in Section ^ 
The relation between Ig and the KNO function is 

I, - ln(n) = ^ In r{z)dz. (32) 

Asymptotically, KNO scaling therefore corresponds to scaling of the ratio 
Ig/ ln(n), analogous to the result discussed above for the ratio S/ \ia{n). 



The entropy approach is also valid for smaller regions of phase space ||33| , |34 | 
and can be formulated in terms of the fractals discussed in Section Here 
we merely note that a multi-fractal character of the multiparticle production 
process requires violation of KNO scaling, and predicts a widening of the 



KNO function f{x) and an increase of its maximum as energy increases [33 



3.3 Conventional distributions 

The purely phenomenological study of multiplicity is based on the use of well 
known functions from probability theory. Here we consider three functions 
for which analytic expressions for the generating functions and all moments 
can be derived . These will serve as a starting point for our discussion 

of the QCD results presented in Sections ^^j. 



3.3.1 The Poisson distribution 

The presence of correlations in a process is conventionally described by the 
difference between its typical distribution and the Poisson distribution, which 
is written as 

P„ = ^exp(-(n)). (33) 
The generating function is (see eq. (|D) 

G(^)(2) =exp((n)z), (34) 
where, from eqs. (^) and (^, 

Fg = l, Kg = Hg = 6g,. (35) 

Therefore the measure of correlations can be defined as the difference be- 
tween Fq and 1, or, equivalently, between Kg (or Hg) and 0. The Poisson 
distribution exhibits exact F scaling and asymptotic KNO scaling. 
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3.3.2 The negative binomial distribution and its generalizations 

The negative binomial distribution (NBD) deserves special attention. The 



NBD has been used to describe experimental measurements of multiplic- 
ity in a wide variety of processes and over a large energy range. In particular, 
the NBD qualitatively describes the distribution of multiplicity in almost all 
inelastic, high energy processes, except for data at the highest available en- 
ergies. The NBD is given by 




r(n + i)r(fc) 



(36) 



where k is an adjustable parameter with the physical meaning of the number 
of independent sources. The Bose-Einstein distribution is a special case of 
the NBD with k=l. The Poisson distribution is obtained from eq. in 
the limit k—^oo. 

The generating function is 

G^^^^\z)= (l-^) (37) 
The integer rank moments are 

_ r(fc + 9) 
- g^. (39) 

For a fixed value of k, the rate of increase of the factorial moments with q 
is larger than exponential. As q increases, the cumulants at first decrease, 
reaching a minimum at g ~ fc, after which they increase. The cumulants are 
all positive. The product of several generating functions of negative binomial 
distributions with different parameters also yields positive cumulants since 
the unnormalized total cumulant is just the sum of unnormalized individual 
cumulants. Similarly, the ratio Hq is always positive. The Hq moments 
decreases monotonically and are proportional to q~^ at large q. 
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Fig. |1] shows the behavior of InFg, In Kg and InHg as a function of q for 
k=5 and 10. The function P„ becomes narrower as k increases. This explains 
the slower rise of Fg for k=10 compared to k=5. The dependence of Kg on k is 
seen to be more pronounced than for Hg. These properties are characteristic 



of the NBD and not of QCD (see e.g. Section [Sp . 

The negative binomial distribution with fixed parameter k exhibits F 
scaling and asymptotic KNO scaling. For large n, the KNO function with 
fixed k behaves as 

m = (41) 

where x=n/{n). This distribution is wider than the Poisson distribution for 
all k. The generating function eq. (|37D is singular for z=k/ (n) ^ if (n) — > oo 
and k is constant. Therefore, it is necessary to approach the singularity when 
calculating the derivatives of the generating function at z=0, i.e. the factorial 
moments (see eq. (^)). The singularity moves closer to z=0 as the energy 
increases. 

The generalized negative binomial distribution can be represented as the 
Poisson transform of the generalized gamma distribution ^8\, given by the 
special scaling function (see eq. (PO)) 



/(x) = JJ^X^'''x''''-^ exp(-[Ax]'^) (42) 



with three adjustable parameters. The Poisson transform of eq. (42) leads to 



a multiplicity distribution which can be expressed in terms of different Fox 



functions in various intervals of the parameter fi, as was studied in Ref. [§£ 
The notion of fractional moments is quite useful when considering the gen- 
eralized NBD because its integer order moments are formally equivalent to 
the fractional moments of the ordinary NBD (see Ref. [^). 

The so-called modified negative binomial distribution with the generating 
function ^ 

has also been used to describe data [^. The number of adjustable param- 
eters is again three. They are related to the average multiplicity through 
the relation {n)=k{x ~ ^)- The negative binomial distribution corresponds 
to A=0. Such fits of multiplicity usually yield better results than the NBD, 
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20 q 



Figure 1: Moments of the negative binomial distribution ||3^ for k=5 and 
10, calculated for integer values of rank q: (a) InFg , (b) In Kg , (c) InHg. 
The curves are drawn to guide the eye. 
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which is no surprise given the larger number of parameters. On the other 
hand, it is difficult to give a physical interpretation to these parameters. The 
cumulants derived from eq. (E3|) are 



which change sign at each subsequent integer q for the range of parameters x 
and A given in . This differs from QCD results described below. Moreover 
negative probabilities appear at large n which is forbidden. 

A three parameter fit also has been attempted in the framework of the 



Saleh-Teich distribution with the claim that it perfectly describes the 
observed shape of the distribution at high multiplicities, including the shoul- 
der (see Section |375|) . Moments of the multiplicity distribution have not been 
studied using this approach, however. 

A combination of two NBDs has also been used in some fits. This also 
enlarges the number of adjustable parameters. 

We will also mention the log-normal distribution which describes the 
Gaussian distribution of the logarithms of the variable. The discrete proba- 
bility distribution for charged particle multiplicity n is defined by 

Ar r+''" dn' , [ln(n' + c) - 

Pn{i2,(7,c) = N ^— exp(-^^ — — ^ '-), (45 

Jn n' + c Za"^ 

where /i, a and c are adjustable parameters and is a normalization factor. 
The integration parameter is 6n=2 for full phase space. In restricted rapidity 
intervals, where both even and odd multiplicities contribute, 6n=l. This 
distribution corresponds to a scale invariant stochastic branching process 
(for details see Refs. H, H). 



3.3.3 The fixed multiplicity and Gaussian distributions 

We next consider the example of fixed multiplicity. Fixed multiplicity demon- 
strates how severely experimental selection criteria can infiuence moments. 
Events with a given multiplicity have often been studied in the past, the 
so-called semi-inclusive events. In this case one deals with the distribution 

Pn = Snno {no = COUSt) . (46) 

The generating function is 

G(^)(z) = (l + 2)-o. (47) 
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Since {n)=no, one obtains 




1 < g < no, 



(48) 




Fg = 0, q> no, 
no)^-''(g-l)! = (-no)^-T(g) 
:(-l)i-%o5(g,no-g + l). 



(49) 
(50) 
(51) 



The factorial moments decrease monotonically with q until q=nQ. All fac- 
torial moments with rank larger than uq are zero. Thus the Hg ratios are 
defined for q < Uq only. The integer order cumulants alternate in sign for 
successive ranks, being positive at odd values of q and negative at even val- 
ues: therefore the cumulants "oscillate." As q increases, the amphtude of 
the oscillations decreases until q=nQ, after which the amplitude increases. 

Oscillations of the cumulants are also observed in QCD but with a dif- 
ferent periodicity (see Section p.2| ). Furthermore, in QCD, the factorial mo- 
ments increase rapidly with q. 

Fig. shows the Fg, Kg and Hg moments of the fixed multiplicity distri- 
bution for no =10. The insets show In \Kq\ and In 

For fixed multiplicity, the oscillations of the cumulant moments are en- 
tirely a consequence of the event selection procedure. This is easy to recognize 
when one selects, for example, 10-particle events from a set of events which 
are distributed according to the Poisson distribution. Then we obtain alter- 
nating sign cumulants at integer ranks instead of cumulants equal to zero. 
The amplitude of the oscillations preserves information about the original 
distribution as long as the normalization is not changed. 

A similar behavior of cumulants is observed for the Gaussian distribution 



If cr oc (n), the Gaussian exhibits KNO scaling. Its generating function is 



For 0" ^ 0, the generating function tends to the fixed multiplicity case, i.e. the 
cumulants are positive for odd ranks and negative for even ranks. For cr^ > 
(n), the cumulants are positive for even ranks and negative for odd ranks. 




(52) 




(53) 
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Figure 2: Moments of the fixed multiplicity distribution |^6|] for no=10, cal- 
culated for integer values of rank q: (a) Fg , (b) Kg ,\n\Kq\ , (c) Hg , In 
The curves are drawn to guide the eye. 
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3.4 Some models 



The distributions discussed above are well known from probability theory. 
In particle physics, they have been used in the context of cluster or clan 
models and in the multiperipheral cluster model P| with one particle ex- 
change. 

At first sight, the graph-theoretic description of multiparticle produc- 
tion seems entirely different between e^e~ and hadron initiated processes. 
In the former case, the main graphs are tree-like, with a highly virtual ini- 
tial time-like parton. In the latter case, one used to consider a sequence 
of multiperipheral-type graphs with low space-like virtualities and rather 
complicated topologies. A more general and unifying picture emerges by 
considering strings to be stretched between the color charges (cf. the Lund 
model PSl, HH], dual topological model E^, H3], and the quark-gluon string 



model [^, 0) and by introducing models for final particle clusterization 
(the multiperipheral cluster model [^, clans 0, etc.). The multiplicity 
distributions in these models are not usually described by a single analytic 
formula but are formed from a combination of distributions. For example, 
the multiperipheral model with a single ladder (chain, string) is based on the 
Poisson distribution for the particle emission centers (resonances, fireballs, 
clusters, clans, etc.). In general, the resulting distribution is obtained by con- 
volution of the Poisson distribution, which specifies the number of sources, 
with the decay distribution of the sources. If one chooses a logarithmic dis- 
tribution to describe the decay multiplicity, the distribution of final particle 
multiplicity is given by the NBD. This model does not incorporate energy- 
momentum conservation, however, which would almost certainly modify the 
simplicity of this result. 

The simultaneous creation of several ladders yields a more complicated 
distribution. Particles produced in different chains are treated as uncorre- 
lated. For each configuration, the multiplicity distributions of different lad- 
ders are convoluted. Averaging over all multiple scattering processes yields 
the final multiplicity distribution. 

The multiplicity is sometimes approximated by a sum of negative bino- 
mials with different parameters, resulting in a distribution with "shoulders" 
or "quasi-oscillations" imposed on a smooth curve. (The possible relation 
between these oscillations and those of the Hq moments discussed below in 



Sections |5.2| and |11.5| is considered in Refs. ||52|, 0.) In particular, the mul 
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tiplicity distribution within a single jet in e^e^ annihilation is sometimes 
approximated by a single NBD, while the superposition of several jets yields 
a distribution with a shoulder |CT in qualitative agreement with experiment 
(see Section [STSp. However, the solutions of QCD equations do not support 
the hypothesis of superimposed NBDs, as we discuss below. The detailed 
study of semi-phenomenological models is usually performed using Monte 
Carlo computations, involving assumptions about hadronization. It is possi- 
ble that the discrepancy between the QCD solutions and the results of the 
phenomenological approaches is related to these assumptions. 

In hadron-induced reactions, one generally describes multiparticle pro- 
duction as proceeding through a set of clusters created through the exchange 
of some mean transferred momentum between the colliding constituents. 
The general kinematic relation between the primary energy s, the cluster 
masses Sj, and the effective transferred momentum is easily obtained by 
iteration of the formula s = S1S2/PI [0. For equal cluster masses and trans- 
ferred momenta, this relation can be written as s/p1 ~ (si/pf)'^ or, in terms 
of mean multiplicities, as 

{n{y)) = ^— ^(n(l/,,)) ^ —^{n{m)), (54) 

where ?/ = In s, yi = In Sj, and Lp = \np^ (with s, Sj and p^ appropriately nor- 
malized by, say, 1 GeV^), and N= {n{y)) / {n{yi)) is the number of clusters. 
The traditional multiperipheral approach assumes constant values of and 
Lp and therefore constant {n{yi)). Thus {n{y)) oc y, i.e. the mean multi- 
plicity is predicted to increase logarithmically with energy. A faster than 
logarithmic rise is observed experimentally, however. According to eq. (|5^), 
this implies that the ratio of Si to p^ tends to unity at high energies, an effect 
which is also observed by experiment. The model of classical fields as applied 



to nucleus-nucleus collisions (see Refs. 0) should yield similar features 
if interpreted in terms of clusters. 

Another important consideration in the comparison of theory with hadron- 
hadron data is that the phase space is relatively unoccupied in typical events, 
i.e. for events with multiplicities close to the mean, because the phase space 
volume is proportional to the primary energy ^/s whereas the mean mul- 
tiplicity increases more slowly, as exp(cA/logs) according to perturbative 
QCD or as a power of logs according to some phenomenological models. At 
present, we describe these data using peripheral- like ladder type processes 
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with reggeons or reggeized gluon exchange. The phase space is elongated 
along the axis of the impinging hadrons and fixed in the transverse direction. 
For the tail of the multiplicity distribution at very high multiplicities, the 
phase space becomes saturated. Which dynamics governs this region is not 
yet clear. Surely, the longitudinal momenta will decrease compared to typ- 
ical events because of conservation laws. However it is unclear whether the 
transverse momenta are the same or larger than in than in typical events. 
This latter possibility has been suggested by some theoretical speculations. 
In the former case, one should employ non-perturbative approaches while 
in the latter case perturbative solutions can be found. Additional experi- 
mental and theoretical effort is required in this area, which could result in 
the observation of new physical effects or at least a deeper understanding of 
the dynamics. The Fermi statistical model |5^ and the Landau hydro- 
dynamical model (or their modifications), popular in the 1950s, could 



apply here. A statistical approach to the description of very high multiplicity 
processes based on the generating function technique is proposed in ||6T| . 

The treatment of nucleus-nucleus interactions is even more complicated, 
involving an extremely large number of exchanged ladders between numer- 
ous colliding nucleons, yielding a much wider multiplicity distribution than 
hadron-hadron collisions. This feature is in qualitative agreement with ex- 
periment. In central heavy ion collisions, particle production is described 
by multi-ladder exchange in the reggeon approach. This implies that the 
number of partons emitted by colliding nuclei is very large, and that their 
combined action can be treated as a source of some classical field which cre- 
ates the final-state particles . The multiplicity distribution is a Gaussian 



type 1^ and exhibits KNO scaling, i.e. aoc (n) in eq. (|52D. The cumulants 
change sign at each subsequent rank. However, such a classical field model 
is found to be applicable only to the tail of the observed multiplicity distri- 
bution. 

Despite its greater complexity, collisions involving nuclei share some fea- 
tures with the simpler hadron-hadron and e"'"e~ collisions, as we discuss in 
the following sections. 

3.5 Experimental phenomenology 



When KNO scaling was first proposed [|12| , experimental data were available 
at comparatively low energies only. The available data from hadron collisions 
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approximately satisfied the KNO scaling condition. A slight deficiency of the 
fits at low multiplicities was cured by use of the correct relation between dis- 
crete and continuous distributions . It was only when UA5 data became 
available in 1985 with energies of several hundred GeV, that violation of 
KNO scaling was clearly observed. The negative binomial distribution was 
then used to fit the UA5 data. However, a shoulder structure was present in 
the data which could not be reproduced using a single NBD. 

The situation in e^e~ annihilations followed a similar path, with sugges- 
tions at low energies that KNO scaling was valid followed by observations 
at higher energies that it was not. A shoulder structure in the multiplicity 
distribution, similar to that observed in hadron collisions, was also observed 
here. The shoulder is ascribed to the multi-source nature of the processes, 
either from additional ladders in hadron reactions or to multi-jet emission in 
the e"'"e~ events. 

A universal energy dependence of average multiplicities in reactions ini- 
tiated by different particles was proposed |6^, 3. It was argued that the 
scale for hadron-initiated reactions should be the so-called effective energy, 
obtained by subtracting the energy of leading particles from the total event 
energy. The effective energy corresponds to the notion of inelasticity coef- 
ficients in cosmic ray studies. The mean multiplicities of e"'"e~ and hadron 
reactions were found to coincide up to ISR energies (a/s~50 GeV) if the 
effective energy was used in the latter case. This correspondance has not yet 
been checked for higher energies, however, where it is known that multiplicity 
sensitive quantities such as rapidity differ between the two processes. 

With the advent of high statistics e"'"e~ experiments at the T? energy, it 
became possible to analyze multiplicity not only in full phase space but in 
limited phase space regions. As an example. Fig. ^ shows results from the 
ALEPH Collaboration at LEP for the distribution of charged particle multi- 
plicity in hadronic Z° decays |65]. The data are shown for full phase space and 
in limited phase space intervals defined by rapidity |Y| < 2.0 and |Y| < 0.5. 
The data are compared to the results of the negative binomial and log-normal 
distributions. Overall, the log-normal result is in somewhat better agreement 
with data than the NBD. The shoulder structure, especially pronounced in 
the data for multiplicities n ~ 30 in the |Y| < 2.0 rapidity bin, is not well de- 
scribed by either model, however. Concerning factorial moments, the NBD 
underestimates the experimental results whereas the log-normal distribution 
overestimates them |^6|, Further analysis of factorial moments reveals 
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disagreement between the data and the log-normal distribution both in small 
intervals and full phase space |67|. This is an example of how moment anal- 



ysis can reveal small differences between theory and data in a more clear 
manner than a direct fit of multiplicity. Neither the modified negative bino- 
mial distribution nor the so-called pure birth stochastic model ||68| improve 
the situation. 




Figure 3: The measured charged particle multiplicity distribution in e^e 



hadronic Z° decays ^5|, for small (|Y| < 0.5, left) and intermediate (|Y| < 
2.0, center) rapidity windows, and in full phase space (right), in comparison 
to various models. 



The predictions of the JETSET [3£] and HERWIG |7C] parton shower 
Monte Carlo event generators are included in Fig. ^ These two programs, 



along with the ARIADNE parton shower Monte Carlo |71[, are widely used 
to interpret e"'"e~ data (see, for example, Ref. [0). The Monte Carlo results 
are seen to describe the data quite well, presumably because they account 
for hard gluon emission relevant for the shoulder. Monte Carlo studies based 
on the quark-gluon string model and the dual topological model similarly 
provide a rather successful guide to experimental distributions from hadron- 
and nucleus-initiated reactions. 
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As a final remark we note that other phenomenological distributions have 
been used to describe experimental data on multiplicity, such as the so-called 
modified gamma distribution [ |75| . 

4 Equations of Quantum Chromodynamics 

Multiparticle production in Quantum Chromodynamics arises from the in- 
teractions of quarks and gluons. The interactions lead to the creation of 
additional quarks and gluons and eventually to the formation of hadrons. 
The most characteristic features of QCD processes are determined by the 
vector nature of the massless gluons and by the dimensionless coupling con- 
stant. Gluons carry color charge and therefore emit gluons in addition to 
quark- ant iquark pairs. The development of quark and gluon jets is described 
by evolution equations. The main parameter of the evolution is the angle 
of divergence of the jet or, more precisely, its transverse momentum. The 
essential difference between Quantum Electrodynamics and Quantum Chro- 
modynamics - due to the triple gluon coupling - can be traced to the lowest 
order graphs in which two colliding electrons (or quarks) emit a single photon 
(or gluon). In the latter case there exists an additional graph with a gluon 
emitted by the exchanged gluon. The interference of this graph with the 
others leads to transverse momentum ordering, in contrast to pure angular 
ordering as in Quantum Electrodynamics (see e.g. Ref. W^ ). 

The subsequent emission of gluons and quarks within the jet fills in the 
internal regions of cones defined by the previous emissions (transverse mo- 
mentum ordering). This remarkable property of "angular ordering" can be 
exploited to formulate a probabilistic scheme for the development of the jet 
as a whole. In this case the evolution equations are reminiscent of the well 
known classical Markovian equations for "birth-death" or "mother-daughter" 
processes. (For a detailed discussion of this approach, based on the coherence 
phenomenon, see Refs. 0, [75[|). 

The system of two equations for the generating functions Gp and Gg of 
quark and gluon jets, respectively, is (with A, B, C = F, G) 0, ^ ^ 

GA{y,z) = e-^^^y^ 

\ ry rl 

+ dxexp[-WA{y) + WAiy')] 
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X ^K^''ix)GB[y' + Hx),z]Gc[y' + Hl-x),z], (55) 
Zn 

where 

y = HpQ/Qo) = H2Q/Qo) = ln(Q/A), (56) 

with p the initial momentum, the angle of divergence of the jet (the ini- 
tial jet opening angle), assumed here to be small, Q the jet virtuality, Qq 
a constant which defines the limit of perturbative evolution, and as the 
strong interaction coupling strength or "constant." A is the so-called scale 
parameter of QCD, related to Qq by A=(5o/2.Q The scale parameter A is 
strongly related to Aj^ |j79|, Q but is not, in general, the same. We shall 



consider both running and fixed as- An effective infrared safe coupling con- 
stant is sometimes used as a substitute for the phenomenological parameter 



Qq (see HI)- Here, we use the more traditional approach, however. 

The first term on the right hand side of eq. (|55|) corresponds to the prop- 
agation of the primary parton without interaction, and is described by the 
form factor exp[—WA{y)]- The second term describes the creation of two jets 
B and C which carry proportions x and 1 — a; of the primary energy after 
their production at the vertex K^'-' , which is reached by the primary par- 
ton after it evolves to scale y' without splitting, as dictated by the factor 
exp[—WA{y) + WA{y')]- Iteration of eq. (|55D generates all tree graphs of the 
perturbation series, ordered by transverse momentum. 

Considering just terms of order {as log^ s)", where s ~ is the cms en- 
ergy squared, and then summing, leads to the so-called DGLAP evolution 
equation (see Ref. |8^) with kernels explicitly written below. The approx- 
imation in which only these terms are taken into account is the double- 
logarithmic approximation (DLA). Adding further terms with lower pow- 
ers of logs yields the modified leading-logarithmic approximation (MLLA). 
Adding yet lower powers of logs, one obtains the next-to-next-to-leading or- 
der approximation (NNLO), and so on. Even the lowest order approximation 
includes processes with an arbitrary number of partons and is described by a 
sub-series of the perturbative expansion with different powers of the coupling 
constant. Thus the classification of approximations is not based simply on 
powers of as but rather on the product of a§ with the corresponding power 
of log s. This explains why terms from different sub-series are combined when 

^ Note that the relationship between Qq and A is not unique; other relationships have 
been chosen (see e.g. Refs. U^) and they are often treated as independent variables. 
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the conservation laws (e.g. energy conservation) are included, as discussed in 
more detail below (see Section ^.2|) . 

Multiplying both sides of eq. ( |53| ) by exp[wA{y)] and then differentiating 
by y, we eliminate the form factors and obtain the final system of equa- 
tions i, a 

G'a = rdxKSixh'o[GGiy + \nx)GGiy + \na-x))-GGiy)] 
Jo 

+ rif f'dxK^{x)^'o[GF{y + lnx)GF{y + ln{l-x))-GG{y)], (57) 
Jo 

G'p = dxK^ix)^^o[GGiy + \nx)GF{y + ln(l - x)) - Gpiy)], (58) 
where G'{y) = dG/dy, Uf is the number of active quark flavors, and 



7o = (59) 



TT 



The kernels of the equations are 



K^(^x) = {l-x)[2-x{l-x)], (60) 
K^{x) = + {1 - xn (61) 



K^(x) = ^ 



1 X 

1 + - 

X 2 



(62) 



with N^=3 the number of colors and CpHN^ - l)/2Ar^=4/3 in QCD. 

The variable z has been omitted in eqs. ( p7D and (pH]). It should be 
remembered, however, that derivation of the equations for moments relies on 
the expansions (0) and (H) of the above equations. 

A typical feature of any field theory with a dimensionless coupling con- 
stant, and of QCD in particular, is the presence of singular terms at x — >^ in 
the kernels of the equations. These singularities imply an uneven sharing of 
energy between the newly created jets and play an important role in the jet 



""The system of eq. (6.29) in Ref. |g| is the same as the system of eqs. ( ply and (|58|) in 
the present work, using a sHghtly different notation; the notion of the generating function 
in this context was first introduced in Ref. [M as is also mentioned in Ref. M . 
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evolution, giving rise to a larger average multiplicity compared to the case 
of equal sharing of energy (the non-singular case). 

The system of equations ( p7D and ( |5^ ) is physically appealing but not 
absolutely exact. This last point is clear since, for example, there is no 
four-gluon interaction term in eqs. (^) and (^8]). Such a term does not con- 
tribute a singularity to the kernels, justifying its omission from the lowest 
order approximation. Despite this lack of exactness, the modified series of 
perturbation theory (i.e. with three parton vertices only), is well reproduced 
by eqs. (|57|) and (^8]) up to terms including higher order logarithmic correc- 
tions. As shown in Ref. |0, the neglected terms contribute at the level of the 
product of, at least, five generating functions. The physical interpretation of 
the corresponding graphs leads to treatment of the "color polarizability" of 
jets. 

There are complications associated with the definition of the evolution 
parameter, with preasymptotic corrections, etc. (see, e.g., Ref. |8^). For 



example, the limits of integration in eqs. (^) and (|58|) are in principle con- 
strained by a restriction on the transverse momentum: 

kt = x(l - x)pQ' > go/2. (63) 

This condition originates from the requirement that the formation time of 
a gluon (tform ~ k/k^) be less than its hadronization time (thad ~ kR'^ ~ 
k/Q"^). This leads to the requirement that the arguments of the generating 
functions in eqs. (|57|) and (|58D be positive, which in turn implies that the 



limits of integration in eqs. (^71) and (58) should be e ^ and 1 — e ^ rather 



than and 1. These limits tend to and 1 at high energies. Therefore it 
seems reasonable to use integration limits of and 1 in eqs. ( pTj) and (|58|), 
to study their analytic solutions with high accuracy, and then to account for 
the neglected terms by considering them to be preasymptotic corrections to 
these solutions. 

This issue of the limits of integration has physical significance. With lim- 
its of e~y and 1 — e~^, the partonic cascade terminates at the virtuality Qo, 
as seen from the arguments of the multiplicities in the integrals. With limits 
of and 1, the cascade extends into the non-perturbative region with virtual- 
ities smaller than Qo/2, i.e. Qi ^ xpQ/2 and Q2 ~ (1 ~x)pQ/2. This region 
contributes terms of the order of e~^, power-suppressed in energy. It is not 
known whether the equations and the LPHD hypothesis are valid below the 
cutoff Qq. 



32 



5 Gluodynamics 



It is natural to begin our study of QCD with the case of gluodynamics, in 
which there are no quarks and only the interactions of gluons are considered. 
Gluodynamics exhibits all the qualitative features of QCD evolution while 
being more transparent. In gluodynamics, the system of equations ( ^Tl) and 
( |58D reduces to the single equation 

G'{y)= dxK{x)-f^[G{y + In x)G{y + \n{l - x)) - G{y)], (64) 
Jo 

with G{y) = GG{y) and K{x) = Kq{x). Eq. (|6^) is a non-linear, integro- 
differential equation with shifted arguments in the nonlinear term which ac- 
count for energy conservation in the triple QCD vertex, when one gluon splits 
into two. As mentioned above, kt ordering influences the limits of integration, 
and instead of and 1 in eq. (0) we should insert and 1 — e~^. Use of 
these latter limits leads to preasymptotic power corrections which we neglect 
for the moment and estimate later. Even though momentum conservation is 
not, strictly speaking, included, it is taken into account in an approximate 
manner through a combined action of energy conservation and transverse 
momentum ordering. 

In the lowest order (double-logarithmic) approximation, only the most 
singular terms in the kernel K{x) and within the square brackets are retained. 
Thus K ^ 1/x and ln(l — x) —>■ (therefore one neglects energy conservation 
in the asymptotic limit), while 7q is chosen to be constant. 

5.1 Approximate solutions of equations with fixed cou- 
pling constant and the shape of the KNO function 

In Section || we describe the exact solutions of the QCD equations for the 
moments of multiplicity, assuming a fixed coupling strength. We show that 
the properties of gluon jets do not change appreciably if quarks are included. 
Therefore it is instructive to examine the more transparent approximate so- 
lutions of eq. (|6^) . 

Formally, the assumptions of the double-logarithmic approximation are 
equivalent for the three factors in the integrand of eq. (0) since non-leading 
terms are neglected. In Quantum Chromodynamics with a fixed coupling 
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strength, F scaling is favored over KNO scaling (see Section |^). The dif- 
ference between F and KNO scaling is usually neglected, however, since the 
calculations are often performed for asymptotically high energies. Preasymp- 
totic corrections for the second moment of the multiplicity distribution are 
discussed in Ref . , and for higher moments in Ref . . 

The generating function for the lowest order solution of eq. (|6^) is inde- 
pendent of both the energy and coupling constant (see Refs. 0, l^)- In this 
case eq. ( |64D reduces to the differential equation (for more details see eq. (^) 
below) 

[\nG{y)r = ll{G{y)-l). (65) 

The corresponding KNO function /(x) defined by eq. (|20|) decreases expo- 
nentially at large values of x: 



fix) ~ 2C{Cx - 1 



3Cx 



+ ...) exp(-Cx); Ca; > 1, 



where C ~ 2.552. For small values of x its behavior is 



/(x) ~ x ^ exp(— In^ x/2). 



(66) 



(67) 



The appearance of asymptotic KNO scaling and its independence of the cou- 
pling constant in the lowest order approximation are by themselves a success 
of perturbative Quantum Chromodynamics [0. The shape of the scaling 
function (^) does not fit experiment, however. Experiment favors shapes 
which are much narrower than predicted by eqs. (|66|) and (|67D. The correc- 



tions of the modified leading-logarithmic approximation yield a function that 
reducing the width of /(x) and introducing a dependence 
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is less wide 
on as- 
Detailed studies of the solutions of gluodynamics have been presented in 
many papers 0, |14| ||36|, |Q, In most cases only moments 

of low rank are considered, i.e. the average multiplicity and its dispersion. 
The approximations in these papers do not always allow higher rank moments 
to be treated, leading to unphysical results such as negative factorial moments 



of the fifth rank (see Ref. 



which are forbidden by the definition, eq. (^). 
This illustrates the importance of consistency when treating terms of the 
same order in QCD equations. The role of conservation laws in the shifted 
arguments of the generating functions is extremely important, providing large 



corrections. It was shown in Ref. |T3| that conservation laws can be taken into 
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account. However, in Ref. the running property of the couphng constant 
was disregarded, the non-singular terms in the kernel were neglected (as 
were some other terms), and the difference between the coupling constant 70 
(eq. (|59D) and the QCD anomalous dimension 7, defined as 

{n)=exp{J%{y')dy'), (68) 

was neglected as well. 

Eqs. (0) and (|58D possess exact solutions for a fixed coupling constant, as 
discussed in Section ^j. To find these solutions requires the numerical solution 
of some algebraic equations (see Section |^). Therefore we first consider the 
approximate approach of Ref. flSl because it yields an analytic expression 
for the KNO function which reveals the importance of the conservation laws, 
while differing from eq. (|66D by predicting a narrower width, in much better 
agreement with experiment. 

First, one obtains a recurrence relation for the factorial moments by sub- 
stituting expression (|^) into eq. (0) and equating the coefficients of z'^ on 
both sides: 

9-1 

(g - q-')F^ = 7 E C^Bijl, j{q - I) + 1)F,_;F,. (69) 
1=1 

This system of equations can be solved^ using the initial conditions Fq=Fi=1. 
The inset in Fig. §(a) shows the ratio of the resulting moments to their 



values at large ranks q: 

[1(1+7)]^ 2gr(g + l 



r(l + 7g) C<i 
For large values of 57 one obtains 



where 



^^^r(Ul\, (71) 



^=(1-7)-'. i) = cll^^. (72) 



6 



The exact solution of the system of equations for quark and gluon jets in case of the 



fixed couphng constant is given in Section 3.2 



35 




2 4 6 8 10 12 q 



Figure 4: (a) The ratio of factorial moments derived from the KNO curve in 
Fig. Ilto the asymptotic result of eq. ([70|); the inset shows the corresponding 
ratio between the results of eqs. (|69|) and ([70|) . (b) The ratio Hq obtained 
from the modified KNO curve in Fig. ^ (solid line) compared to its NBD 
counterpart with k=7 (dashed line). 



36 



The asymptotic results for Fq at large ranks determine the asymptotics of 
the KNO function f{x) at high multiplicities: 



f{x) ^ '^jl^^^eM-{Dxn (/i- » 1. (73) 

Thus the tail of the distribution at large multiplicities is far more suppressed 
than in the double-logarithmic approximation, eq. (|66D . Using the value 
of 7Ri0.4 at presently accessible energies for which /i=(l — 7)"^ ~ 1-6, one 
obtains a Gaussian-like suppression rather than an exponential falloff as in 
eq. (|66D. Thus, we conclude that the conservation laws drastically reduce the 
width of the multiplicity distribution. 

This is demonstrated in Fig. H, where the modified QCD distribution, 
i.e. the distribution including energy conservation, is compared with the low- 
est order QCD result. The modified distribution also accounts for some 
corrections at low multiplicities [|13|. The modified distribution (thick solid 



line) is seen to be much narrower than the DLA distribution (thin solid line). 
Also shown in Fig. ^ is the result of a fit of the modified distribution by the 
NBD, with parameters k=7 and (n)=30. 

Making use of the modified QCD curve in Fig. |^, the genuine factorial 
moments can be calculated, i.e. the moments including the corrections at low 
multiphcities. Their ratio to the asymptotic solution, eq. (|70|) , is shown in 
the main part of Fig. ^(a). Comparison of the two curves in Fig. ^(a) reveals 
the importance of the corrections at low multiplicities. Using the genuine 
factorial moments, one can compute the cumulants and the Hq ratios. The 
latter are shown in Fig. 0(b) and compared with the Hq moments obtained 
from the fitted NBD curve in Fig. ^ A dramatic difference is seen between 
the behavior of the modified QCD solution and the NBD. The Hq ratios in 
the former case are seen to oscillate as the rank increases. In contrast, the Hq 
ratios of the NBD decrease monotonically, as was already seen in Fig. 0(c). 
These oscillations of the Hq ratios are reminiscent of the oscillations observed 
for the fixed multiplicity distribution. Fig. 0(c). Oscillations with a differ- 
ent periodicity are obtained from the solutions to the full QCD generating 
equations, as discussed below in Sections |^ and ^. 
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Figure 5: The modified KNO function (tliick solid curve) for 7=0.4 is 
much narrower than the lowest order distribution (thin sohd curve). The 
negative binomial distribution with k=7 is shown for comparison (dashed 
curve). 
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5.2 Higher order approximations with a running cou- 
pUng constant 



Eq. (0) for the generating function in gluodynamics can be solved in a 
more general approximation |]14|, by accounting for all terms of the kernel 



K including the non-singular ones, by incorporating the running coupling 
constant, and by using a Taylor series expansion of the generating functions 
in the nonlinear term at large y: 

G{y + 6) ^ G{y) + G'{y)e + \G"{y)e' + ... (74) 

This approach illustrates the distinction between the various assumptions 
and their relative importance. 

Each derivative in eq. (ff^ produces a derivative of (n) which is propor- 
tional to the anomalous dimension 7 according to eq. (|68D. In turn, the 
anomalous dimension 7 can be represented by a perturbative series with 70 
as the leading term. Therefore, higher derivatives in the Taylor series give 
rise to terms of a higher order approximation in the perturbative expansion. 
The Taylor series leads to a recombination of different order terms in this 
expansion. It is important to note that energy conservation, accounted for 



by X and 1 — a; in eqs. (|57|) and (0), is an essential and integral aspect of the 
higher order perturbative corrections. At first sight, it may seem odd that 
the order of the approximation can be changed when energy conservation is 
included. However, this becomes clear if one considers that energy sharing 
among partons results in terms with lower powers of log s at the same power 
of i.e. in a shift to a higher approximation as seen from eq. ([7^) . 

The running coupling strength as should be considered in the two-loop 
approximation when working in NNLO. It is given by the formula 



27r 



o^s{y) = i^,[l-^,--^]+0{y~\ (75) 




with 

,3„ = ft = 17iV|-.v(5iV. + 3C.)^ 

where n/ is the number of active quark flavors (in gluodynamics nj=0), and 
the difference between the QCD scale A and (5o/2 is neglected in asymptotics. 
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In MLLA, one can use the one-loop expression, i.e. the first factor of the 
product in eq. (|75D only. 

Note that - in contrast to QED - the QCD series is an expansion in 
the product of the coupling strength and energy dependent (Ins) terms, 
rather than an expansion in the coupling strength alone. Therefore higher 
order corrections appearing in the QCD expansion arise not only from higher 
loop graphs but also from a proper accounting for conservation laws, as 
mentioned above. In this review, we show that the effects of the latter 
- resulting in a limitation of the available phase space - are much more 
important than the two-loop terms of the coupling strength. This explains 
e.g. why NNLO calculations of the multiplicity ratio r between gluon and 
quark jets yield rather different results depending on the accuracy with which 
energy conservation is included, as we discuss below in Sections and ^ 

Inserting expression ([Tip for the generating function into the nonlinear 
term of eq. (|6lD, dividing both sides of (0) by G{y) and differentiating with 
respect to y, we obtain 



[\nG{y)]" = ^',{y)[{l + 4hh',{y)){G{y) 

oo oo 

+ E(-ir/^nG'(")(y)+ E (-1)' 



\m+n 



hn 



1) - 2hG'{y) 

/(^(m)^(n)\ ' 



n=2 



m,n=l 



G 



(77) 



where 



hr 



h = 


11/24, 


h2 = 


(67 - 67r^)/36, 


h = 


(47r2 - 15)/24, /i4 = 13/3, 


h — 


|2-2-"-3-"-C(n)|, 


(in) = 


oo 

^ m"", n > 2, 

m=l 


1 

m\n\ 


/^da;ir(x)ln"xln"(l-x 

JO 



(78) 



(79) 



The first term on the right hand side {G{y) — 1) without its pre-factor is the 
well known expression of the double-logarithmic approximation, eq. (|65|) . 
The additional higher order contribution 4/i^7q is due to the derivative of Jq, 
and is considered in [93| but not in [Q. Similar corrections appear in the 



40 



subsequent terms of eq. ([771). In particular, the term proportional to hi cor- 
responds to the modified leading-logarithmic approximation, while the term 
proportional to /i2 together with the 4/if7Q correction to the first term, men- 
tioned above, corresponds to the next-to-next-to-leading order correction. 
Actually, the transverse momentum kt given by eq. (|63D should be used as 
the argument of as in place of y, but it can be shown that such a replacement 
leads to yet higher order corrections than those considered here. Note that 
the use of the preasymptotic limits of integration, i.e. and 1 — e~^, do 



not lead to a differential equation which is as simple as eq. (|77|) but result in 
power-like corrections. 

A straightforward solution of eq. ( [TTD appears problematic even if only 
the terms with hi and /i2 are included in addition to the double-logarithmic 



terms. The solution is very simple for the moments of the distributions [|14 
however, since G{z) and \nG{z) are the generating functions for factorial mo- 
ments and cumulants, respectively. Using formulas (|^ and (Q) and assuming 
asymptotic F scaling, one obtains the product 57 (and its derivatives) at 
each differentiation in eq. (|7^) because the average multiphcity is the only 
^/-dependent term which remains. The coefficients of z'^ on both sides should 
be equal. Hence, one obtains in NNLO: 

rr _ 7g[l - 2hiq^ + Ahhi + h2{qh' + ^70] ...^ 



The anomalous dimension 7 is defined by eq. (68). The condition Fi=Ki=l 
determines the relation between 7 and 70: 

7 ^ 70 - \hi^l + ^(4/^2 + 15/i?)7o' + (^(7o), (81) 

which shows that the increase of the average multiplicity with energy is 
slower in the modified leading-logarithmic approximation than in the double- 
logarithmic approximation, since the term with hi is negative (see eq. (|68D). 
However, the higher orders reverse this situation. The MLLA and NNLO 
terms in eq. ( |8TD practically cancel each other for present values of 70. There- 
fore the series (0) essentially yields the asymptotic result. Only at very large 
energies when 70 ~ do subsequent terms in the series become smaller than 
the lower order ones (if only the first term in the expression for 05 is consid- 
ered) . However, one should keep in mind that in NNLO the value of 70 itself 
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acquires a correction, as given by the second term in the brackets in eq. ([75|). 
This correction is negative and drastically reduces the value of the coefficient 
in front of 7q in eq. (|HT|), resulting in a much better convergence of the series. 
This is discussed below in connection with the NNLO and 3NL0 corrections 
to the expression for the energy dependence of mean multiplicities. 

The running property of 70 was taken into account in expression 
according to eq. (|75|) : 

7o^-/^i7o + (^(7o), (82) 

which leads to 

7' ^-/ii7o'(l-/^i7o) + 0(7o')- (83) 

The lesson we learn from eq. (|80D is that in all correction terms (which 
contain hi, h2, . . .), the expansion parameter 7 appears in a product with 
the rank, namely the product 57, which becomes large at high ranks, i.e. at 
high multiplicities. Therefore, for high multiplicity events one should take 
into account ever higher order terms in 7. This problem was mentioned a 
long time ago and is discussed in some detail in Ref. fO^, but has only 
recently been analyzed. 

As already mentioned, the double-logarithmic formulas are obtained from 
eq. dT?! ) by setting hi=h2=0, 7=70, and 7'=0. In this case, 

H, = q-\ (84) 

which is similar to the asymptotic form of the negative binomial distribution 
with k^2, corresponding to an extremely wide multiplicity distribution (see 
eq. (|66D). In contrast, experimental data typically yield values of k in the 
range from 3.5 to 100. We stress, however, that the NBD result for k=2 is 

which shows that the "genuine" correlations are about twice as large in the 
NBD as in QCD. Therefore the QCD predictions for a single jet do not sup- 
port the hypothesis of the negative binomial distribution even in asymptotics, 
not to mention other distinctive features of eq. (|5D| ) and its generalization. 

It is easy to recognize that eq. (Q) is just an asymptotic form (for ?/ — >■ 00) 
of eq. (BUD ^^"^ fixed rank q. The value of 7 tends to zero as y"^^"^ at 



large y. In this limit one obtains eq. (84) from eq. (BQ). The approach to 
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the asymptotic limit is very slow, however. The omitted terms are of order 
(9(ln~^/^ s). Therefore they should be taken into account at present energies. 

Note that the preasymptotic difference between assumed F scaling and 
KNO scaling does not alter our conclusions because it is of order C((n)~^), 
i.e. it decreases with energy faster than the terms and con- 

sidered. 

A more interesting feature of the solutions is the qualitative behavior of 
Hq as a function of rank at a given energy, i.e. when 70 is kept constant in 
eq. (|80|) . According to eq. (0), the Hq moments are a smoothly decreas- 
ing function of q in the double-logarithmic approximation. In the modified 
leading-logarithmic approximation in which the linear hi term is kept but 
the /i2 term and higher order terms are neglected in eqs. ( [77| ) and (^Dj), Hq 
acquires a minimum at 

^ +1 + 0(70). (86) 



/ii7o 2 

Note that the minimum position is shifted to larger rank values at higher 
energies due to the decrease of 70, and to smaller ranks for lower energies or 



for smaller phase space windows if higher order corrections are neglected 
(see Section |14D. At present energies, with a value of 70 of about 0.45-0.5, 
we obtain 

qmin ~ 5 (±1). (87) 

The uncertainty of ±1 for this result is due to the 0{jq) correction in eq. (|86|). 

In MLLA, the ratio Hq crosses the abscissa axis, is negative at the min- 
imum, and tends to zero from below as ~ —q~^ as the rank increases. If 
one includes the NNLO terms (proportional to /i2), the ratio Hq is shifted 
up by an amount /i27o independent of q. The ratio then exhibits a second 
zero, crosses the abscissa axis again, and tends asymptotically to a positive 
constant /i27o- The location of the minimum shifts to slightly larger values 
of q because of the term with h\ and the two-loop correction to as- If the 
the exact limits of integration and 1 — are used instead of and 1, 
the minimum moves back to smaller ranks, however. Thus we conclude that 
the location of the minimum as approximately given by eqs. (|86D and (|87D is 
a rather stable prediction of gluodynamics. Furthermore we see that eq. (^) 
predicts an oscillation of the Hq ratio analogous to that found using the more 



approximate methods of Section p]T 
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One could object that the expansion parameter 57 entering eq. ( pOD is 
large at present energies. This, however, is similar to the situation with the 
expansion of, for example, cosx in a Taylor series, for which the higher order 
terms approximate the function quite well even for large values of x. As 
more terms of the Taylor expansion of cos x are included, the approximation 
improves at yet larger values of x. 



A computer solution of eq. (u7i) |96] also yields an oscillating behavior for 



Hq when account is taken of the higher order terms with hs and hn, see Fig. |[ 
This is consistent with the discussion in the previous paragraph. However, 
the position of the first minimum is shifted to g ~ 4 in the approximation of 



Ref. ||96[, illustrating the sensitivity of Hq to the various assumptions. In par 



ticular, the term with h^, important for the proper limit in supersymmetric 



(SUSY) QCD as shown below, is not included in the results of Ref. . 

Although the assumptions of the calculation influence the quantitative 
details of the oscillations, the main qualitative features are stable. An implicit 
cutoff of the multiplicity distribution at some maximum value due to the 
restrictions of finite energy produces similar oscillations, although with a 
smaller amplitude |^ |9^. Our present goal is to discuss the qualitative 
properties of Hq. More precise expressions are presented in the following 
sections. We emphasize that the amplitudes and periodicity of the oscillations 
are different in the more precise treatment than in the results shown in Figs. ^ 
and |. 

In this section, we have demonstrated that conservation laws and higher 
order perturbative corrections lead, in the framework of gluodynamics, to a 
substantial reduction in the width of the multiplicity distribution and to a 
qualitative change in the behavior of the cumulant and factorial moments, 
namely to oscillations of their ratio Hq as a function of q which disappear 
only at extremely high energies. 



6 Perturbative solutions of QCD equations 
with a running coupling constant 

A discussion of Quantum Chromodynamics, in which quarks are created as 
well as gluons, leads back to the system of coupled equations (^) and (^) 
for the generating functions. The structure of these equations does not differ. 
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Figure 6: The ratio Hq as a function of q reveals oscillations in higher order 
perturbative QCD [Q . The results shown here correspond to the Z° energy. 
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however, from the equation for gluodynamics, eq. 



Therefore we shall 
W 
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not present the results for the full solutions (see e.g. Refs. ^ 
but will merely describe them. 

In analogy to gluodynamics, one obtains a system of coupled recurrent 
equations for factorial moments and cumulants when the Taylor series ex- 
pansion (ff^) and the formulas and (H) are used. This system of equations 



has been solved analytically for average multiplicities, accounting for correc- 
tions up to C(7o), i.e. the 3NL0 approximation |98[ 



Analytic calculations 
of the second and third factorial moments have been performed in the same 

Numerical solutions are available for yet higher mo- 
The properties of gluon jets do not change appreciably 



approximation 
ments 19 
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m, El- 

compared to Section^, i.e. the gluon jet cumulants and factorial moments are 
very similar to those determined in gluodynamics. In particular, the ratio Hg 
of gluon jets exhibits a minimum at the same value g ~ 5 found in Section ^ 
The factorial moments of quark jets are determined to be larger than those 
of gluon jets, i.e. the parton multiplicity distribution of quark jets is wider 
than that of gluon jets even though the average multiplicity is smaller. This 
fact was known long ago ||101|| - |[l03[| . The quark jet cumulants and Hg ratios 



are also found to oscillate. The first minimum of the quark jet oscillation is 
predicted by the analytic solutions to be located at g~7-8, but it moves to 
g:^5 if the non-asymptotic limits of integration are incorporated [|104|| . 

To apply these results to the two jet process of electron-positron annihila- 
tions, it is necessary to relate the e"'"e~ generating function to those of quarks 
and gluons. Bearing in mind the Feynman diagram for the production of two 
quark jets in e"'"e~ annihilations at a very early stage, one obtains 



(88) 



with further corrections (see, e.g., Ref. [Q). Thus the zeros of the quark 
jet and e^e~ cumulants coincide since the logarithms of their generating 
functions are proportional to each other (see eq. (|^)). 



6.1 Average multiplicities and their slopes 



We now describe the QCD results for average multiplicities ||9^, |100|| in more 
detail. The equations for the average multiplicities in jets are obtained from 
the system of equations (^7^ and ( |5^ ) by expanding the generating functions 
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in z and keeping the terms with g=0 and 1 in the series (|^ and (|^). They 
read 



(ncd/))' = y dx-il\K%[x){l(nG{y + Inx)) + {nciy + ln(l - x)) - {ndy))) 

+nfK^{x){{nF{y + Inx)) + {npiy + ln(l - x)) - (^^(y)))], (89) 

My))' = I dx-ilK'f{x){{nG{y + \Tix)) + {nF{y+Hl-x))-{nF{y))). (90) 



From these equations one can predict the energy evolution of the ratio of mul- 
tiphcities r and the QCD anomalous dimension 7 (the slope of the logarithm 
of the average multiplicity in a gluon jet), defined by 

The perturbative expansions of 7 and r are 

7 = 7o(l - ai7o - a27o - d^ll) + O(7o), (92) 

r = ro(l - ri7o - r27o - ^37o) + 0{^t). (93) 

The asymptotic value of r is rQ=Nc/CF=9/4:, as was first obtained in ||105|| . 
Using a Taylor series expansion of (n) at large y in eqs. (|9]) and (pO|), in 



conjunction with ( p2D and (0), the coefficients and can be determined 
analytically. Their numerical values are given in Table [l| for various numbers 
of active quark flavors and for SUSY QCD. In SUSY QCD where nf=Nc=CF, 
all Tj are identically zero. 

The value of the coefficient ri was first reported in Ref. ||106|] and later in 



Ref. [g^]- The coefficient r2 was first presented in Ref. ||8^. The r2 results 



in Ref. do not account for energy conservation in the equations for the 



generating functions, however, and as a consequence are much smaller than 
the r2 values given in Table The larger values of r2, i.e. including energy 
conservation, were first presented in Ref. P3. Similar results to those in 



Ref. 1 93 1 were subsequently obtained in the framework of the Lund dipole 



model ||107| with cutoff fractal triangles. Large NNLO corrections are also 



advocated in Ref. using somewhat different reasoning. The energy con- 
servation effects responsible for the large values of r2 are closely related to 
the oscillations of higher rank moments discussed below. The 3NL0 values 
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ri 


r2 


^^3 


ai 


(32 


as 


3 


0.185 


0.426 


0.189 


0.280 


- 0.379 


0.209 


4 


0.191 


0.468 


0.080 


0.297 


- 0.339 


0.162 


5 


0.198 


0.510 


-0.041 


0.314 


- 0.301 


0.112 


S 











0.188 


-0.190 


-0.130 



Table 1: Numerical values of the perturbative corrections up to order 7q 
for the multiplicity ratio r and the QCD anomalous dimension 7, based on 
integration limits of the generating functions from to 1 (see text), nj is 
the number of active quark flavors while S refers to supersymmetric (SUSY) 
QCD. 



of 03 and were first presented in [03 where one can also find analytic 



expressions for all and rj with 2 < 3. 

According to eq. (^) and the values of in Table 0, r at higher orders is 
much smaller ^3| than in the double-logarithmic approximation. On average, 
it is smaller by about 20%. At the Z° resonance, the subsequent terms in 
eq. (0) diminish the theoretical value of r compared with its asymptotic 
value by approximately 10%, 13%, and 1% for ?t,/=4, see Section p.l.2| and 

Fig. ig. 

Inserting eq. (^) into eq. (|68|), the energy behavior of the average multi- 
plicity in gluon jets at 3NL0 is determined to be ||100|| 

{nc) = Ky--^^' exp{2c^ + ^[2a2C^ + §(ln2|/ + 2)] 

Po 

+ f[«3C^-^(ln2y + l)]), (94) 

where c={4:Nj PoY/'^={2B)-^^^ and K is a normalization constant. For ^/=0 
and B=hi this reduces to the corresponding result of gluodynamics. The 
pre-exponential term and the first term in the exponent correspond to the 
MLLA expression ||108| , |109|| (for a coupled system of quark and gluon jets, 



these terms were first found in |^). The second term in the exponent, 
proportional to c/ y/y, is the NNLO correction. This term can be also written 
as a small negative correction to the pre-exponent (see Table |T|) . The third 
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term in the exponent, proportional to c^/y, is the 3NL0 result. The role of 
the 3NL0 term in the gluon jet multiplicity is not important compared to the 
lower order terms because of the smallness of a^. Note that the NNLO and 
3NL0 corrections are almost constant and somewhat compensate each other 
at currently accessible energies. As a consequence, the MLLA expression for 
gluon jets is a good approximation to the higher order result, eq. (p^. 

The solid curves in Fig. ^ show the behavior of the average gluon jet 
multiplicity as a function of energy (i.e. y) as calculated in NNLO using a 
running coupling constant. The different curves show the results for ^/=3, 
4 and 5. The parameter Qo=l-3Aj^, with A(^=175 MeV for %=5, and 
taken in the proportions 63:100:130 for n/=5:4:3, has been chosen as in ||110| 



For purposes of comparison, we also show the energy dependence of the 
mean multiplicity for a fixed coupling constant. This is indicated by the 
dashed curves. At low energies, the multiplicity for fixed coupling increases 
rather slowly. At higher energies, the rate of increase exceeds that found 
using the running coupling constant. This is reasonable since the coupling 
strength in the case of fixed coupling has been evaluated at a rather high 
scale, namely yzo ~ 6.6 corresponding to the mass of Z°, yielding a relatively 
small (fixed) coupling strength. In actuality, the coupling strength should 
increase during the course of evolution of the jet, while the number of active 
flavors should decrease. The two trends somewhat compensate each other. 

The slope of r, namely r' = dr/dy, is extremely sensitive to higher order 
perturbative corrections. The role of higher order corrections is increased 
here compared to r because each nth order term, proportional to 7o , gets an 
additional factor of n in front of it when differentiated, the main constant 



term disappears, and the large value of the ratio r2/ri becomes crucial 



(95) 



where the relation 7q ^ — _B7q(H-_Bi7q) has been used, and where B=f3o/d>Nc 
and Bi=(5i/ 4:Nc(3o. The factor in front of the brackets is very small even at 
present energies: i^rgri 0.156 and 70 ~ 0.5. Nonetheless, the numerical 
estimate of r is unreliable because of the expression inside the brackets. 
Each differentiation leads to a factor as (or 79), i.e. to terms of higher order. 
The values of ri, r2 and r-^ given in Table |l| for nf=4 yield 2r2/ri ^4.9 and 
(3r3/ri) + Bi^l.5. Thus the first correction in eq. (|95D (proportional to 
7o) is more than twice as large as unity at present energies, while the next 
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Figure 7: The energy {y=\n{Q / Qq)) dependence of the average multiphcity 
of gluon jets, for a running (sohd curves) or fixed (dashed curves) couphng 
The results are shown for different numbers of active flavors, 



constant 
nf= 3, 4 



93 
and 5. 



correction (proportional to 7q) is about 0.4. Therefore even higher order 
terms are necessary before the perturbative result for r' can be considered to 
be reliable. 

The slope r' is for a fixed coupling constant. In the case of the running 
coupling, r' evolves rapidly with y in the lowest perturbative approximation, 
ranging according to eq. (^) from 0.06 at the Z° to 0.25 at the T. 

Next we consider the ratios of slopes r^^^ and curvatures defined as 

Their ratios to r can be written as 

T T 
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r z^rr + rr — 2r . . 

p^ = ^ = ^ (^fTTp ■ ^ ^ 

Since r' oc7o, the asymptotic values of r, r^^^ and r*^^) coincide, equaling 2.25. 
Moreover, the values of r, r^^^ and r^^^ coincide in MLLA. The first preasymp- 
totic corrections are very small. They are of order O(7o) (NNLO) with a 
small factor in front and contribute about 2-4% at present energies: 

pi = 1 - Bn-f', ^ 1 - 0.077o', (99) 

P2 = l-25ri72^ 1-0.147^. (100) 

However, this favorable situation does not persist when the next terms are 
considered. For pi the numerical values of the higher order terms are 

pi = 1 - 0.077o'(l + 5.3870 + 4.217^) {ri f = 4), (101) 



and similarly for p2 (see Ref. |^8[). All the correction terms are much smaller 
than unity. However, the numerical factors inside the brackets are so large 
that at present values of 70 ~ 0.5 the subsequent terms are larger than the 
first one and the sum of the series is unknown. The series can be summed 
using the simplest Fade approximant (i.e. by assuming the steady decrease 
of the unaccounted terms in proportions determined by the ratio of two last 
calculated terms), resulting in 



pi = 1 - 0.077 



^ ^ 5.3870 



1 - 0.7870 



(102) 



In contrast to pi, the perturbative corrections to the ratio of slopes r^^^ 
are small. The lowest order O(7o) correction to r^^^ is the same as for r 
but higher order corrections are smaller because they are negative both in 
r and pi which define r^^) = r/pi. This explains why experimental values 
of r(^) are similar to values of r calculated in the MLLA approximation (see 
Section |ll]), whereas experimental values of r are about 25% lower than 
the MLLA prediction. Yet more substantial cancelations of higher order 
terms should occur for r^'^\ It would be interesting to test this prediction 
experimentally. 

An interesting feature of eq. ( p,01| ) is that pi contains terms up to C(7o) 
whereas r in the numerator of eq. (^) is known to O(7o) only. This situation 
is related to the fact that r'/7 ~ 0(7q). 
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The value of pi also determines the ratio of the slope of the logarithm of 
the average multiplicity of a quark jet (7^) to that of a gluon jet (7): 



jF = (In(n^))' = Pi7 = 7 - -• (103) 

r 

Once again we see that the logarithmic slopes of quark and gluon jets are 
equal in MLLA. They differ in higher orders in such a manner that 7f < 7 
since both r and r are positive. Our failure to obtain a precise estimate 
of pi implies that we cannot reliably evaluate either. Thus the ratio of 
logarithmic slopes pi is more sensitive to higher order perturbative terms 
than the ratio of slopes r^^\ 

A high sensitivity of multiplicities in quark jets to higher order terms 
is especially noticeable if one writes down the expression for their energy 
dependence, analogous to eq. (Pl) for gluon jets. Taking into account per- 



turbatively the ratio r{y) one obtains |10Q 



{np) = -y-'^^^' exp(2cVy + ^[n + 2a2C^ + §(ln 2y + 2)] 

+-[a,c' + ^ + r2 - ^(ln2y + 1)]). (104) 

Especially striking is the role of the last term which cannot be neglected even 
at the energy of the Z°. 

Another representation of this dependence is the direct use of formula (|9lD : 



ro(l - ri7o - r27^ - r^-^l) 



However one easily observes that formulas ( |104| ) and (|105|) differ since there 
is no term with in eq. ( |1U4| ). This is an important point which emphasizes 
that the notion of the order of approximations is different in the definitions 
of 7 and r. The MLLA term of r (proportional to 70), given by ri, appears 
in the NNLO term of 7i? (proportional to 7q) in combination with 02. This 
situation is a consequence of the fact that r ~ C(7o). Therefore the notion 
of orders is somewhat a question of convention. We conclude that it is im- 
proper to use the term with ra in eq. ( p.05| ) until the 04 contribution to 7 
is known. By extension this implies that if the MLLA formula is used to 
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describe multiplicities, as is common practice, the DLA result r = ro should 
be inserted into all formulas to be self consistent within the perturbative 
approach, leading to the same energy dependence for quark and gluon jets. 
The NNLO term with ri may be included in 'fp only if the 02 contribution 
has been incorporated in 7, and so on. 

Therefore within the present accuracy of O(7o) corrections, the pertur- 
bative QCD approach fails in the precise determination of the logarithmic 
slope of the quark jet multiplicity •jp, even at the Z°, and can be trusted 
only at much higher energies. 



6.2 Widths of the distributions 

Next we consider higher order corrections for the moments of multiplicity. 
The normalized factorial moments of any rank q can be obtained by differen- 
tiating the generating functions according to eq. or, equivalently, by using 
the series (0) and collecting the terms with equal powers of z on both sides 
of eqs. and (|58|). The QCD equations for the moments of multiplicity 



are given in Refs. B3, 85, pq, 99 



The normalized second factorial moments of gluon and quark jets, F2, are 
given by 

The F2 moments define the widths of the multiplicity distributions, being 
related to the dispersion = (n^) — (n)^ by 

= (F2 - + {n) = K2{nf + {n). (107) 

The perturbative expansions of F2 up to 7q are 

= ^(1 - /i7o - Mo - /37o), (108) 

F( = (1 + y)(l - 0170 - 027o - Mo)- (109) 



Inserting expressions (|108D and ( |109| ) into the QCD equations yields pre- 
dictions for the coefficients /, and (pi. The asymptotic (70^0) values of 
and F2 can also be determined by equating the leading terms in 70 on 
both sides of the equations. We have presented their explicit expressions in 
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/i 


/2 


/3 


01 


02 


03 


3 


0.364 


-0.0279 


0.795 


0.637 


-0.276 


2.12 


4 


0.358 


-0.0457 


0.740 


0.631 


-0.286 


2.04 


5 


0.352 


-0.0629 


0.689 


0.625 


-0.295 


1.95 


S 


0.313 


0.310 


-0.120 


0.313 


0.310 


-0.120 



Table 2: Numerical values of the perturbative corrections up to order 7q for 
the normalized second factorial moments of gluon (/j) and quark (0j) jets. 

is the number of active quark flavors while S refers to supersymmetric 
(SUSY) QCD. 



eqs. ( |108|) and (|109|) to simplify later notation. Analytic formulas for /j and 



0j are given in Ref. [Q. Their numerical values for different numbers of 
active flavors and for SUSY QCD are listed in Table |[ From these values it 
is seen that the series ( |108| ) and ( 109 ) change in sign at each term and that 



higher order terms are more important for the width of a quark jet than for 
a gluon jet. 

The asymptotic values of and F.[ are 0: 

FSas = l, C. = l + y = i- (110) 

At lower energies the widths are slightly smaller because of the increase of as, 
leading to larger values of the first corrections (with coefficients /i and (pi) 
in eqs. (|0|) and ([T09|) . 

The MLLA prediction for the energy dependence of the second factorial 
moments is shown in Fig. using nf=4 and the values of as in Ref. [Q. The 
curves demonstrate that the widths are almost Poissonian at low energies, 
i.e. they approach unity as the scale decreases. The small dots at the ends 
of the curves show the effects of a change in n/ to 3 at Q=10 GeV and 
to 5 at Q=90 GeV: it is seen that the results are almost insensitive to the 
number of effective flavors. Note that n/ is the only free parameter in these 
results. Cutoff of the integration at £:=e~^ ~ g-27r//3oas fj^om below and at 1—e 
from above as a consequence of preasymptotic effects is also not particularly 
important at present energies, as shown by the dashed curves in Fig. |^. The 
preasymptotic corrections slightly increase the widths and reduce the slope 
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Figure 8: The energy behavior of the second factorial moments of quark 
(F) and gluon (G) jets ||8^. The scale of the energy Q is GeV. The limits 
of integration are chosen to be either and 1 (solid curves) or e and 1 — e 
(dashed curves). The dots at the ends of the curves demonstrate that the 
results are insensitive to variation of the effective number of flavors (see text). 
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with respect to Q. 

Unfortunately, higher order corrections do not improve the theoretical 
estimates. The NNLO terms are positive, while the 3NL0 corrections are 
negative and so large they lead to sub-Poissonian widths (F2 < 1) for as'=0.2. 
Indeed, the 3NL0 corrections are just as large as those in MLLA. The origin 
of the large values of /s and (f)^ can be traced to rather large contributions 
of integrals of the form In" xdx oc n!, i.e. to the region of very soft gluons. 
Thus the cutoff at and 1 — e"^ should be quite important for the higher 
order terms. This is analogous to the situation with renormalons (see, e.g., 
Ref. pT| ). 



The 3NL0 corrections are overestimated in the method of the Taylor 
series expansion because of the assumption ?/^|lnx| which is invalid for 
soft gluons. For example, the kt dependence of the coupling strength is 
transformed so that 



"50C r~- 1 ^ -]■ (111) 

y + \nx[l- x) y \ V / 

The second term becomes infinitely large for x—^0. The cutoff at x=e~^ 
yields a factor of 2 only. Thus the above expansion implies a presumption 
about the behavior of the coupling strength in the non-perturbative region 
as well as its modification at the limits of the perturbative one. 

The slopes of the widths are especially sensitive to higher order terms be- 
cause each of them is multiplied by the factor n when differentiating 7q . Thus 
the 3NL0 term in the perturbative expression for the slopes of the widths 
is about 3 times larger than the MLLA term. This suggests that a precise 
quantitative estimate of these slopes is not possible. At present energies, 
F^'{MLLA)^OM and Ff (MLLA) ^0.092 for as=0.2. The asymptotic 
value of the ratio of slopes can be estimated to be 

(^2) as _ 2^ ^ Q 43 (112) 

which surely coincides with the MLLA prediction for this ratio. Thus the 
second factorial moment of quark jets approaches its asymptotic value faster 
than the corresponding quantity for gluon jets. 

We stress that all slopes and curvatures in pQCD are related to the 
running property of the QCD coupling constant since they are proportional 
to its derivatives which are zero for a fixed coupling strength. 
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h 


h2 


h 


91 


92 


93 


3 


0.986 


-0.342 


2.49 


1.61 


-1.58 


7.74 


4 


0.972 


-0.380 


2.36 


1.60 


-1.59 


7.54 


5 


0.957 


-0.417 


2.25 


1.59 


-1.60 


7.34 


S 


0.844 


0.722 


-1.09 


0.844 


0.722 


-1.09 



Table 3: Numerical values of the perturbative corrections up to order 7q 
for the normalized third factorial moments of gluon (hi) and quark (gi) jets. 

is the number of active quark flavors while S refers to supersymmetric 
(SUSY) QCD. 



One can easily check (see Ref. [^) that the relations of SUSY QCD, 
where nf=Nc=CF, are valid for the coefficients given above (e.g., F^=F2 
etc.). The asymptotic SUSY values of F2 equal 4/3 for both gluon and quark 
jets. 



Third moments of the multiplicity distributions have been calculated 1^5 
in a manner analogous to that described above for the second moments. The 
results are 

3 3 
F,^ = h,{l -Y.h^ll)■, F[ = g,{l -Y.9^%), (113) 

i=l i=l 

where /io=9/4, go=l + ro + r^/i, and the values of the coefficients hi and gi 
are listed in Table ^ In the asymptotic limit, the third moment of quark jets 
is about twice as large as the third moment of gluon jets. 

Comparing and (pi to hi and gi, one concludes that the corrections 
increase for higher moments, even in MLLA. Moreover, at present values of 
7o ~ 0.5, the corrections are rather large. Note the similarity in the structure 
of the corrections for the second and third moments. The corrections alter- 
nate in sign, and the third coefficients are larger than the first ones. This is 
an indication of a sign-alternating asymptotic series, and Borel summation 
could be effective. 

In SUSY QCD the asymptotic values of the third moments equal 9/4. The 
first correction, given by /ii(SUSY)=5fi(SUSY)=0.844, is about the same size 
as for ordinary gluon jets and is similar to the correction for second moments. 
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However, the NNLO and 3NL0 terms for moments in SUSY QCD differ 
drastically from those for ordinary jets, both in absolute value and sign, as 
seen from Tables |^ and |^. This demonstrates the sensitivity of the results to 
the value of tq, which is different in the two cases. 

A similar procedure can be used to determine yet higher rank moments. 
However, the series will probably be badly convergent. 



7 Computer solutions 

The analytic approach described in Sections ^ and ^ accounts for energy 
conservation in an approximate manner. Energy conservation can be in- 



cluded more accurately by solving eqs. (|57D and ( pSf ) numerically, i.e. by 
implementing a computer solution. The computer solution also allows the 
preasymptotic limits of integration to be directly incorporated, i.e. the limits 
e~y and 1 — rather than and 1. Thus the non-perturbative region of x 
values near and 1 can be avoided. 



A computer solution of eqs. (^Tj) and (^) is presented in Refs. ||104| , ^ 



Thresholds for new quark production are approximated using a smoothed 
function. The coupling constant is considered in the one-loop approxima- 
tion. A deficiency, shared by the analytic calculations, is the neglect of strict 
conservation of transverse momentum and a corresponding limitation on the 
available phase space. This deficiency effectively corresponds to a modifica- 
tion of eqs. (|57| ) and (|58D. 

The principal result of the computer solution is that the value of r is re- 
duced by more than 10% compared with the 3NL0 prediction. The computer 
solution for r versus energy scale Q is shown in Fig. ^. The experimental 



results shown in this figure are discussed in Section |Tl. 

A perturbative cutoff Qo=507 MeV is found from a fit to the measured 
mean particle multiplicity in quark jets [^. With increasing resolution, 
more and more jets are resolved, and predictions for the evolution of the 
subjet multiphcity (see Section |T^) can be obtained as well. The subjet mul- 
tiplicities are found to agree quite well with experiment using this same value 
of Qq ||78|. Correlation functions of high ranks up to g=18 in quark and gluon 



jets were also determined, fully accounting for energy conservation ||104|| . A 
quantitative description of the existing data on factorial moments, factorial 
cumulants, and their ratio Hq was obtained for single quark and gluon jets 
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Figure 9: Computer solution |78[ for the ratio of mean hadron multiplicities 
between gluon and quark jets (solid curve) in comparison to data. The 
experimental result at yc~3 x 10~^ is based on e"'"e~ T? events from the 
OPAL Collaboration ||112| , |113|| . The results at y, 

ITTl 



recorded by the CLEO Collaboration 



0.01 are from T decays 
The dashed curve shows the 



analogous prediction [|78[ for subjet multiphcities. 
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and for two-jet e+e events, using the same parameters determined for the 
average subjet and particle multiphcities, see Fig. |10] and Section |11.4 . 



8 Exact solutions of QCD equations with a 
fixed coupling constant 

The QCD results described above demonstrate that conservation laws and 
the non-singular terms of the kernels play a more important role up to MLLA 
than the running coupling constant, which becomes essential in NNLO only. 
Thus the running property of the coupling constant does not alter the qual- 
itative predictions of QCD for multiplicity. It is possible to obtain an exact 
analytic solution of eqs. (|5^) and (pSf), i.e. without resorting to a perturbative 



approximation, if the coupling constant is fixed I^B], |117|| . No other assump- 
tion is necessary. One can obtain the general solution for the moments of 
any rank in this approximation. The effect of the running coupling constant 
can be assessed by varying the value chosen for 70. This analytic approach 
provides further insight into the behavior of the solutions. 

In the following we first discuss the lowest rank moments obtained in 
the fixed coupling approximation, and compare them with the results of the 
previous sections. Following this we discuss the solutions for higher rank 
moments. 



8.1 First moments and the ratio of average multiplic- 
ities in gluon jets to those in quark jets 

The equations for average multiplicities, i.e. unnormalized moments of rank 
q=l, are obtained by substituting the series (|^ for the generating functions 
into eqs. (|57D and ( p8|) and equating the terms linear in z, with the conditions 
Fq=Fi=^q=^i=1. The factorial moments of quark jets are denoted by 
and their cumulants by If the coupling constant is fixed, the average 
multiplicities behave like 

(^g(2/)) oc exp(7y), oc r"^exp(7?/), (114) 
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Figure 10: Computer solution for the ratio Hg of cumulant over factorial 
moments as a function of rank q (solid curve) ||104|| in comparison to exper- 
imental results from e"'"e~ decays from the DELPHI Collaboration ||115 



as analyzed in [|116 
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where the anomalous dimension 7 and the ratio r are constants. This behav- 
ior follows directly from eqs. (^) and (|58|) as can be seen using the relations 

(^g,f(2/ + lna;))/(nG,F(z/)) = a;^, (115) 

{ncAy))' = l{nG,F). (116) 

Then eqs. (0) and (|58|) can be rewritten as a system of algebraic equations 
for the variables 7 and r: 

7 = 7o[< + n;r(Mf - Mq^)], (117) 

7 = 7o'[L2-Lo + rLi], (118) 

where 



Mf = 


dxK^lx^ + {1 - xy - 1], 

Jo 


Mf = 


f dxK^lx"/ + {1 - xy], 

JO 


< = 


['dxK^ = M[ij = 0)/2, 
Jo 


Li = 


[ dxK^x^, 
Jo 


L2 = 


[ dxK'f(l-xy, 
Jo 


Lo = 


I dxK^ = Li{j = 0). 
Jo 



The coefficients Mj and Lj can be expressed in terms of Euler beta functions 
and psi functions and depend only on 7. 

From this, one can estimate the corrections induced by the preasymptotic 
limits of integration (see e.g. Section ^). For example, the leading term in 
gives rise to a correction of order 0(— exp(— 7?/)/7), about 10% of the main 



term at the Z energy. Computer calculations support this estimate [7S]. The 



corrections decrease with energy in a power-like manner similar to higher 
twist effects. It is possible that other effects of the same order, e.g. due to 
instantons, could influence the multiplicity distributions. No definite answer 
to this question has yet been given. 
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It should be stressed that 7 does not equal 70 even in gluodynamics 
(n/=0) because Mf differs from 7"^. An approximate equality between 7 
and 7o is valid for 70 <^ 1 but the perturbative expansion for 7 in the fixed 
coupling approximation, 



7 ~ 70 - /ii7o + 2(^1 + ^2)70 + 0{^t) 



(119) 



differs from the corresponding formula (81) for the case of the running cou- 
pling constant, from which one sees that the first correction is twice as large 
in the former case as in the latter. This explains why the average multiplic- 
ity at low energies increases more slowly for fixed coupling than for running 
coupling (see Fig. |^), and why the increase in multiplicity for fixed coupling 
follows a power law at extremely high energies only. 

The ratio r appears linearly in eqs. ( p.l7| ) and ( |118| ), allowing these two 
equations to be written as 



r(7) = 6(7) 



1_ 



0(7) 



-1 



(120) 



where 



r(7) 



[121) 



a(7) 

6(7) 
c(7) 

d{l) 



^(1) - ^(7 + 1) + 5(7, 1) - 25(7 + 1,2) 
25(7 + 2,1) + 5(7 + 2, 3) + 5(7 + 3, 2) + ^ 



[5(7 + 3,1) + 5(7 + 1,3)], 



12 QNc 

^[5(7,1) -5(7 + 1,1) + ^5(7 + 2,1)], 

^--^1) - ^(7 + 1) - 5(7 + 1, 1) + ^5(7 + 1, 2) + ^ 



N, 



The beta functions are just inverse polynomials of 7 but the above notation 
is less cumbersome. 
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Thus the non-hnear, integro-differential equations ( pT] ) and ( p8D can be 
reduced to the algebraic expressions (|12CI|) and ( |121| ), which yield 7 and r as 
a function of 70 and nj. The dependence of the solution on nj is very small. 

The anomalous coupling 7 can be related to 70 by the simple fitted for- 
mula 

7 = 0.077 + 0.6270, (122) 
or by a more theoretically motivated one based on eq. ( p.l9| ), 

7 = 0.9770 - 0.487o' + 0.27o^ (123) 

fitted in the range of 70 from 0.48 to 0.6. Thus 7 changes very slowly with 70. 
The alternating signs of the correction terms correspond to the tendencies 
observed for the running coupling constant, see Table In itself, the value 
of 7 is not of much interest. Even though it is related to the energy depen- 
dence of the average multiplicity, it is known that the power law increase of 
mean multiplicity for fixed coupling is replaced asymptotically by a slower 
dependence ~ exp (ln^/2^) for running coupling (see Fig. |^). The more re- 
alistic behavior provided by the running coupling constant was discussed in 
the previous sections. 

The ratio r is of more interest because the energy dependence of the 
average multiplicities in gluon and quark jets cancels. Thus the fixed coupling 
prediction for r can be expected to be rather general. A computer solution 
of eqs. ( |120| ) and (|121|) yields the results shown in Fig. |1^. The dependence 



on nj is seen to be very mild. More important, the dependence of r on 70 is 
even weaker than that of 7, and its average effective value is 

r = 1.84 ±0.02. (124) 

This value of r is much lower than the prediction from the double-logarithmic 
approximation and corresponds well with the estimate of r from the higher 
order approximate solution of QCD equations with a running coupling con- 
stant, discussed in Section ^ This suggests that yet higher order corrections 
to eq. ( p3|) will be small. Introduction of the exact (i.e. preasymptotic) lim- 
its of integration would reduce the value of r in eq. ( |124| ) yet further, as 
mentioned in Section 0. 

In a realistic process, the virtuality in a jet degrades as the partons evolve 
toward hadrons, presumably with an associated change in the value of n/. 
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Figure 11: r versus 70 in fixed coupling QCD for n/= 3, 4 and 5 [|117| . 



In tlie framework of the calculations with a fixed coupling strength, the 
dependence on can be addressed with the aid of formulas (|5 



and ([75|). 

The dependence on n/ is found to be very small. Thus, as a jet of partons 
evolves toward lower Q, we need not be concerned with a change in the 
number of active flavors. The parton multiplicity will depend on the evolving 
virtuality through a mild variation of 7, but not enough to invalidate the 
fixed coupling approximation. Certainly, in the ratio of the multiplicities, 
such dependences cancel, yielding a stable result for r. This conclusion is 
supported by the perturbative solutions of the equations with the running 
coupling as has already been discussed (see Table |ip. 



8.2 



Widths and higher rank moments of multipHcity 
distributions in gluon and quark jets 



Relations (|115|) and ( |116| ) allow the solution for higher moments of multi- 
plicity to be obtained in the fixed coupling approximation for any rank q. A 
system of coupled recurrent equations is obtained by substituting eq. (|^) 
into eqs. ( |57D and ( pH]) and by comparing the coefficients of z"^ on both sides. 
These equations can be solved by iteration. The results will not be given 
here (see Only the final analytic expressions for the moments of rank q 

as they relate to the lower rank moments will be presented. For this purpose. 
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let us introduce 

The solution of the equations is 

/, = [a,^,(/,0) +6,T,(/,0)]A-\ (126) 

0, = [c,^,(/,0) +c?,r,(/,0)]A7\ (127) 



where 



S, = EWZfif.-i + nfN[,M,^i], (128) 
1=1 

9-1 

T, = Y.L,M.i, (129) 
z=i 

a„ = ^ + Lo,o — Lq„, (130) 
7o 

6, = UfM^, (131) 

= Lqfi, (132) 

d, = ^-Mf + n^<o, (133) 
7o 

Ag = aqrfg - 6gCq, (134) 

= V'(l) - ^(g7 + 1) + 5(97, l)-2fi(g7 + 1,2) 

- 25(g7 + 2,1) + 5(g7 + 2, 3) + 5(g7 + 3, 2) + ^, 

< = ^[5(^7 + 3,1) + 5(^7 + 1,3)], 

iVS = i?[/7,(g-/)7 + l]-2i?[/7 + l,(g-07 + 2] 
+ E[/7 + 2,(g-07 + 3], 

Ki = ^(S[/7 + 3,(g-07 + l]+i?[^7 + l,(g-07 + 3]), 
L,,, = ^[B[h + l,(q-l)^]-B[l^ + l,{q-lh + l] 



N, 



c 



+ l5[/7 + l,(g-/)7 + 2)]. 
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The above expressions may seem cumbersome but their structure is simple. 
They generahze the formulas of the preceding section to any q. The formulas 
of gluodynamics are obtained by setting n/=Ci?=0 if one retains, in and 
A^^;, the leading terms B{q'y, 1) = l/g7 and B[lj, {q — /)7 + 1], respectively. 
Using the results for 7 and r from the preceding section at a given value of 
7o and Uf, one can determine F2 and $2? then increase g by 1 to determine 
F3 and $3, etc. 

The evolution parameter y disappears from the formulas. A posteriori, 
this means that F scaling, with all dependence on y hidden in the average 
multiplicities {ncriy)), is valid for fixed coupling. This leads to a self con- 
sistent system of algebraic equations where all quantities, including Fg and 
$g, are independent of energy. It should be stressed that F scaling is only as 
accurate as the main equations of the fixed coupling approximation. In fact, 
one should speak about asymptotic F scaling since the limits of x-integration 
in eqs. ( p7D and (|58| ) are asymptotic. 

The results of the calculation, when expressed in terms of Fg and $g, are 
shown in Fig. |l^(a) for 7o=0.48 and nf=5. The solutions are seen to increase 
rapidly with q, more so for $q than for Fg. Since these are normalized factorial 
moments, this once again demonstrates that the multiplicity distribution of 
quark jets is wider than that of gluon jets. The results are insensitive to 
the value of rij. The dependence on the coupling strength is also very mild. 
Indeed, the results are insensitive to whether the coupling strength is fixed 
or running. The corresponding results for cumulants are shown in Fig. [T^(b). 

The ratios of factorial moments to cumulants for gluon and quark jets are 
defined here as: 

= ^, (135) 

Vg = (136) 

The results of the calculations for Hg and rjg in fixed coupling QCD are 
presented in Figs. |T3| and [T^, respectively. 



It is of interest to compare the QCD results in Figs. pr^ - |l^ to those of the 
phenomenological models discussed in Section ^. By comparison of Fig. |T^ 
with Figs. |l] and |^, the QCD results are seen to be more similar to the NBD 
than to fixed multiplicity. In fact, Fg in Fig. 0(a) can be approximated 
by the corresponding NBD result with ^=5. While this characterization is 
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convenient, fits of liiglier rank moments by tlie NBD are not appropriate. 
For example, tlie ratio Hq of the NBD decreases monotonically with g, as 
shown in Fig. |I|(c). In QCD at present energies, the Hq moments become 
negative and then oscillate. These oscillations become even more pronounced 
if the preasymptotic limits of integration are used. Thus the QCD and NBD 
predictions for higher moments differ drastically from each other. 

The negative binomial distribution fits the second and third moments 
well and provides a good overall description of experiment. It does not de- 
scribe features of the data such as the shoulders of the distributions or the 
oscillations of the Hq ratios, however. The fixed multiplicity distribution, 
in contrast to the NBD, does yield oscillations. However the Hq moments 
change sign at each subsequent integer value of q in this case, yielding re- 
sults for the moments and for the period of oscillation in disagreement with 
experiment (see Section [11.51). 



The sensitivity of Hq to the detailed shape of the multiplicity distribution 
is clearly demonstrated by its various qualitative forms in the DLA, MLLA, 
NNLO, and higher order solutions of eqs. ( pT]) and (0). At large ranks the 
moments behave in radically different manners in these different approxima- 
tions, as discussed in Section |5^ . Moreover, the behavior of Hq depends 
strongly on slight variations of the particular shape of the factorial moments 
at low values of q. One can demonstrate how easy it is to obtain oscillations 
of the fixed multiplicity type imposed on the double-logarithmic behavior 



from the following exercise. It is known [^g, |118|| that the large q behavior 
of the factorial moments Fq in the double-logarithmic approximation is (see 
eq. del)) 

F, - (137) 

If one adds a preasymptotic term by replacing the factor 2q by 2g + 1 in the 
numerator (which restores the condition Fq=1 but not Fi=l), one obtains 
an additional term in Hq which imposes oscillations of the fixed multiplicity 
type on the monotonous decrease of the form and the ratio Hq becomes 

- ,(2g + l) ' 

where the second term in the numerator appears because of the newly added 
preasymptotic term. 



70 




Figure 14: Ratio rjq of the quark jet distribution in fixed coupling QCD for 
70=0.48 and n/=5 (a) rjg, (b) ln|?7g|. 
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The above examples illustrate the sensitivity of Hq to the approximations 
made in solving the set of equations for the generating functions. The dif- 
ferent qualitative behavior obtained for different approximations in QCD is 
illustrated in Figs. H, 
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and |1J, while the difference between QCD and 
the phenomenological distributions is illustrated by comparing these figures 
with Figs, m and ^. The experimental cutoff of the tail of multiplicity at 
finite energies induces additional oscillations which have been shown to be 
rather small 1^^, as already mentioned in Section |5.2| . The numerical solution 
of equations which incorporate energy conservation both in the generating 
functions and in the limits of integration yields an extremely good fit of ex- 
perimental data even for high moments [|104|| (Section |11.5| and Fig. |10D, but 
also poses some questions about the potential importance of neglected power 
corrections, as mentioned above. 

Unfortunately, there is still no clear understanding of the physical origin 
of the oscillations, i.e. of their periodicity, amplitude, and dependence on 
rank q (it seems that the amplitude increases and the period decreases with 
increasing q). It has been suggested that the oscillations are related to the 
experimental fact of the shoulder structure of the multiplicity distribution, 
interpreted as originating from the superposition of jets with different flavors 
[TT9| 



and topologies 



However, the insight provided by the solution of 
fixed coupling QCD suggests otherwise. In QCD, the jet initiated by a single 
light flavored parton exhibits Hq oscillations, as long as account is taken of 
conservation laws. In the NBD framework one needs a superposition of at 
least two NBDs, interpreted as the contributions of both light- and 6-quark 
jets, to describe the oscillations ||120|| . An analysis of the light quark jet data 



in Ref. demonstrates that the Hq moments of single light quark jets do 
oscillate, however, with no admixture of 6-quarks ||121j| . This supports the 
QCD prediction but argues against the NBD conjecture. 

The fixed coupling solution demonstrates that the running property of the 
coupling constant is not important for the oscillations. More significant is the 
vector nature of massless gluons which controls the singularities of the kernels 
in the QCD equations. This has been demonstrated ||122|| by comparing 
the QCD results with those of the A0g theory (i.e. for scalar mesons in six 
dimensions). This model exhibits asymptotic freedom and contains a non- 
singular kernel in the evolution equation. There is no minimum at q=5 in 
this theory. Oscillations appear only at much larger ranks than in QCD. 
Perhaps the behavior of Hq will provide some insight into how the equations 
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for the generating functions can be generalized, including the fine effects of 
the interaction of "color monsters" |8^, p.l8|| . However, this is a difficult 
problem to solve. 

Note that the above oscillations occur for integer values of q and are not 
related to oscillations of fractional moments. The latter would impose lower 
period harmonics on the oscillations from integer ranks. 



9 Non-perturbative modifications of QCD equa- 
tions 

The non-perturbative aspects of jet evolution are less clear at the moment. 
They are usually hidden within the hypothesis of local parton-hadron duality 
or in Monte Carlo hadronization models. The influence of the QCD vacuum 
condensate on jets cannot be estimated within the perturbative approach. 
The first proposal to phenomenologically account for its effects within the 
equations describing QCD jet evolution is given in Refs. ||123| , p.24|| . 



The action of the vacuum condensate on partons in jets was assumed to 
be similar to the influence of the medium on electron-positron pairs in elec- 
tromagnetic cascades. As is well known, electromagnetic cascades evolve to 
a large number of particles, of the order of the ratio of the initial energy to 
the electron mass, if the ionization processes in the medium are neglected. 
Because of ionization, low energy electrons stop and exit the cascade. The 
ionization process is described by a non-scaling term in the evolution equa- 
tions. 

A similar term was added to the QCD equations ||123|, |124||. Later on. 



when angular ordering was properly incorporated, it was shown that this 
non-scaling term could also be included in the equations for the generating 
functionals ||125|| . This implies that some energy is spent on the formation 



of non-perturbative strings, diminishing the number of partons in the cas- 
cade. The ratio of partons with and without account of the non-scaling 



[non-perturbative) term is, in the lowest order approximation ||125|| 



where k denotes the typical energy spent for soft non-perturbative interac- 
tions, so that k/p is a small parameter. It is also possible to estimate the 
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energy distribution of partons in the jet [|125|| and the behavior of moments of 

The most interesting property of eq. 



the multiphcity distribution 126 



is its non-scahng behavior, namely it does not depend on the product pQ (i.e. 
on y), but on p and separately. Also, it depends on the product of the 
perturbative and non-perturbative scales. 
A similar approach is advocated in Refs. 



26, 27]. A non-perturbative term 



is added to the QCD equation for the generating function, and is interpreted 
as an analogue of an "inactivation process" implying hadronization. This 
leads to modifications of KNO scaling (see eq. (|25|) ). 

A somewhat different approach to the same problem is formulated in 
Refs. ||127| , |128|| where non-perturbative effects are included in the effective 



Lagrangian, and a Monte Carlo scheme is developed. The general idea is to 
construct an effective field theory that embodies both partonic and hadronic 
degrees of freedom. The effective action is composed of three parts, with 
quark and gluon fields at small distance scales, an effective low energy La- 
grangian with vacuum condensate fields representing hadronic states at large 
scales, and a term which couples these different degrees of freedom at interme- 
diate scales. The final goal is to describe the transition from the quark-gluon 
interaction stage to hadrons within a unified scheme of strong interactions. 

All such proposals aim at the description of the hadronization stage and, 
in particular, of the relation between the numbers of partons and final-state 
hadrons. Unfortunately, it is still necessary to introduce ad hoc phenomeno- 
logical assumptions to achieve this goal. 



10 Singularities of the generating functions 

To this point, we have considered 2 to be a subsidiary variable in the gen- 
erating function G{y, z), which is set to a constant value after calculation of 
the moments. It is of interest to study the behavior of G{y, z) in the com- 
plex z plane. This interest stems from the fact that the singularities of the 
generating function are located close to the point z={] where the moments 
are calculated. For example, the singularity of the generating function of 
the negative binomial distribution is located at z=k/{n), as mentioned in 



Section |3.3.2[ For k=2 this is just a pole of second order at z=2/ (n), which 
tends to z=0 at high energies. 

In QCD the singularities of the generating function are known in DLA 
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only (see [0, p. 137). In this approximation, eq. ([77| ) reduces to eq. (]65|). 
The leading singularity governing eq. (^) can be written as 

G{z,y)^{z-zo{y))-^ (140) 

since differentiating In G twice produces a second order pole on the left hand 
side which must be matched by a singularity of the same order on the right 
hand side. Since KNO scaling is valid in this approximation, the factorial 
moments do not depend on y asymptotically, and the generating function 
depends on the product of z and (n). The singularity is located at 

Zq^"^ ~ C ^ 2.552, (141) 

{n) 

i.e. close to the singularity of the negative binomial distribution with A;=2, 
as it should be. This singularity tends to z=Q at high multiplicities as for 
the NBD. 

The overall structure of the singularities is more complicated in DLA 
QCD than for the NBD, however, because there are poles of second and 
first orders, logarithmic branch points, etc. A more complete DLA result for 
G{z,y) near the singularity is 

Oi.y)^j^.^y,^n-^.Om. (142) 

In a certain sense this result is not very helpful since the DLA is a rather crude 
approximation to reality at present energies. This result shows, however, that 
z=0 is still quite far from the singularity since the pole contributions cancel 
at this point. There is no cancelation in the derivatives of G, however, which 
increase like factorials because of the leading singularity. It is not yet known 
how this structure is modified by higher orders. Attempts to go beyond DLA 
are connected with the Taylor series expansion of G discussed above. Then 
eq. ( |65D acquires additional terms on the right hand side proportional to 
derivatives of G{z, y). Each derivative exhibits a stronger singularity than G 
itself. Therefore such an approach cannot be applied near the singular points. 
This also demonstrates why the DLA has difficulties in practice. However, 
attempts to employ eq. (^) rather than eq. (p5| ) have not yet been done 
explicitly. 



75 




Figure 15: Complex zeros of the charged particle multiplicity distribution, 
obtained [ |133|| from e"'"e~ Monte Carlo events at Ec.m=l TeV, where u = z+l. 



In experiment, only the so-called truncated generating function has been 
determined 

G'^{y,z)= Pn{l + zr, (143) 

n=0 

where the sum extends up to the maximum multiplicity n^nax available at the 
energy scale y. Therefore the truncated generating function is a polynomial 
in z and has no singularity at finite \z\. It possesses rimax complex conjugated 
zeros. When rimax tends to infinity these zeros should converge to the singu- 
larity location of G{y,z). The properties of zeros (the rate of convergence, 
etc.) can be used ||129|| to obtain some knowledge about the singularities of 
the generating function. 

The problem becomes even more intriguing if an analogy with statistical 
mechanics is invoked (see, e.g. [0, |130|| ). In this case z is interpreted as 
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fugacity, and G{y^ z) and C^^ly, z) play the roles of the canonical and grand 
canonical partition functions, respectively. Therefore their properties in the 
complex z-plane can be considered to reflect the statistical properties of 
multiparticle production. In statistical mechanics the motion of zeros to the 



real z axis, "pinching it" asymptotically, implies a phase transition ||131 



132|| . Experiment indicates that such a pinch occurs as the energy increases. 



for various reactions initiated by particles and nuclei. The zeros tend to 
lie close to a circle of unit radius in the complex z plane, as is illustrated 
in Fig. |l^ using e"'"e~ annihilation Monte Carlo events generated at a cms 



energy of 1 TeV ||133||. Similar results have been obtained using simulations of 



nucleus-nucleus collisions generated according to the dual parton model ||134 



Moreover, the quantity analogous to the free energy exhibits its maximum 
in QCD |9^. This fact is closely related to the minimum of the Hq ratios at 



g ~ 5 discussed above. The interpretation in terms of a phase transition is 
the same here. 

Nonetheless, we cannot rely too much on the analogy with statistical 
mechanics to obtain physics conclusions. The properties of the zeros dis- 
cussed above have been shown ||135|| in purely mathematical terms to be a 
consequence of the falloff of the multiplicity distribution at large n and its 
flattening with increasing energy. This behavior forces all zeros to the unit 
circle in the complex 2;-plane, centered at z=— 1, with the point z=0 re- 
moved. Experiment favors these tendencies. The physics implications of this 
behavior are therefore not entirely clear. 



11 Measurements of mean multiplicity, slopes, 
and higher moments 

We now explicitly turn our attention to experimental results. In this, we 
concentrate on data which test the QCD predictions discussed in Sections ^ 
^. We emphasize data from e^e~ annihilations. The relative simplicity of 
hadronic events from e"'"e~ collisions allows a level of precision and conclu- 
siveness difficult to achieve in other types of reactions. For example, ep and 
pp collisions are characterized by initial-state strong interactions and final- 
state remnants from the colliding hadrons. These features - absent in e^e~ 
events and the calculations - make a comparison to theory uncertain. The 
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comparison between experiment and theory is yet more uncertain for nuclear 
collisions. The available measurements of average multiplicities in ep, pp and 



pp collisions are compiled in Ref. |80 



To test the predictions in a meaningful manner, the experimental defini- 
tion of jets should match the theoretical one. The theoretical definition of 
jets is based on the production of a virtual quark- ant iquark qq pair (for quark 
jets) or gluon-gluon gg pair (for gluon jets) from a color singlet point source. 
The inclusive multiplicity of these qq or gg events defines the multiplicity of 
"two-jet events." A single quark or gluon jet corresponds to a hemisphere in 
these events. Thus there is no selection of a specific event topology, i.e. jets 
are not selected using a mathematical algorithm such as the k± [[136| , |137|| or 
JADE | |138| | jet finders, in contrast to common experimental practice. Jets 
defined according to the inclusive theoretical prescription are called "unbi- 
ased." 

It has proven difficult to measure gluon jet multiplicity in an unbiased 
manner because gg production from a color singlet point source is a process 
which is practically unobserved in nature. This difficulty impeded experimen- 
tal progress in the field for many years. One channel where the experimental 
selection criteria match the theoretical definition is the decays of bb bound 
states to 7gg: events like this have been studied through the selection of 
T — >■ 7gg— 'y+hadrons events. Another possibility is rare e^e~ hadronic 
annihilation events in which the quark jets q and q from the electroweak 
decay of the intermediate or virtual 7 are approximately collinear: the 
gluon jet hemisphere against which the q and q recoil in these events cor- 
responds to a nearly unbiased gluon jet [p.39| , |14CI|| . The two photon process 
e"'"e~ — >■ e"'"e~gg^ e^e~ +hadrons has also been suggested as a source for high 
energy unbiased gluon jets ||141|| : this possibility has yet to be explored ex- 
perimentally, however. 

In contrast to gluon jets, it is easy to measure quark jet multiplicity in 
an unbiased manner, using the charged particle multiplicity in hemispheres 
of e^e~^hadrons events. e"'"e~ hadronic annihilations therefore provide a 
natural source for quark jets, studied by many experiments. 

For the results discussed here, the data have been corrected for detector 
response and initial-state photon radiation, treating all charged and neutral 
particles with lifetimes greater than 3 x 10~^° s as stable. Hence charged 
particles from the decays of Kg and weakly decaying hyperons are included 
in the definition of multiplicity. 
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Before proceeding to the main results, presented in Sections |11.1[ - |11.6| , it 
is of interest to recount the original experimental efforts to measure gluon jet 
multiphcity, and in particular the multiplicity ratio r={nG) / {np) (see also 
Ref. [|140|| ). The JADE Collaboration selected three-jet qqg events in e+e" 
annihilations ||142|| . The gluon jet was identified by assuming it had the lowest 
energy of the three jets in an event. Quark jets were the higher energy jets in 
the same events. The UA2 Collaboration studied gluon jets produced in two- 
Mi 



jet gg final states from pp collisions 



These gluon jets were compared 
Til , 



For both 



to lower energy quark jet data from an e"'"e~ experiment 
studies, qualitative indications were reported that the gluon jet multiplicity 
was slightly larger. The interpretation of these results was unclear, however, 
since the gluon and quark jets were not identified using the same criteria in 
either study. 

The HRS Collaboration chose a different strategy, selecting three-fold 
symmetric e"*"e~ qqg events in which the quark and gluon jets were pro- 
duced with about the same energies and inter-jet angles ||145|| . Thus in this 
case the gluon and quark jet identification criteria were equivalent. The prob- 
ability that the gluon jet had a larger multiplicity was tested by assuming 
a Poissonian multiplicity distribution and independent production of each 
of the three jets. A value r=1.29lQ;4g "was derived, the first quoted result 
for r. Although the largeness of the uncertainty precluded any definite in- 
terpretation, this result was considered at the time to indicate that r was 
much smaller than the naive expectation r^2 (see Section |11.6|) . The DEL- 
PHI Collaboration later applied a similar analysis technique to a much larger 
data sample and obtained r ^ 1 ||115|| , however, thus demonstrating that the 



method was not sensitive to gluon and quark jet differences. 

The OPAL Collaboration selected quark and gluon jets in one-fold sym- 
metric qqg e"'"e~ events in which the angle between the highest energy jet 
and each of the two lower energy jets was the same, namely 150° ||146|| . These 
events later became known as "Y events" since they are shaped like the let- 
ter Y. This analysis utilized b-quark tagging to identify a high purity sample 
of gluon jets which were then compared to a sample of about 50% quark and 
50% gluon jets selected using the same criteria. The result r=1.27±0.07 was 
derived from these data, the first positive observation of a larger multiplicity 
in gluon jets compared to quark jets. This result was later confirmed by 
other studies [|147| , |148|| . Although the results of this analysis represented a 



considerable success, establishing a value of r significantly larger than unity 
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Experiment 


Q (GeV) {rich), gg events 


CLEO ||114|| 
CLEO |114 
CLEO 114| 
CLEO pOl 
OPAL |113|| 


4.5 4.88 ±0.10 
5.5 5.28 ±0.10 
6.5 5.65 ±0.12 
10.3 9.339 ±0.090 ±0.045 
80.2 28.56 ±0.36 ±0.62 



Table 4: The mean charged particle muhiphcity, (nch.), of unbiased "two-jet" 
gg events. For the OPAL data and the CLEO measurement at 10.3 GeV, the 
first uncertainty is statistical and the second is systematic. For the CLEO 
measurements below 10 GeV, only the statistical uncertainty is given. 



for the first time, the results were found to depend strongly on the jet finding 
algorithm used to select the events, see Refs. ||147| , |149|| . Since the jets were 
defined using a jet finder rather than employing a hemisphere definition, they 
could not be used for a quantitative test of the analytic results presented in 
Sections ^||, unlike the unbiased gluon jet measurements discussed below. 



11.1 Mean multiplicity 



Experimental measurements of the inclusive charged particle multiplicity of 
unbiased gluon jets as a function of energy scale Q are shown in Fig. |16[ 
The three data points at scale Q ~ 5 GeV are derived from the hadronic 
component of T(1S)— > 7gg events ||1 14|| . The virtuality Q is given by the 
invariant mass of the hadronic system. Similarly, T(3S)— >■ 7X62(10.27)^ 7gg 
events provide the result at Q ~ 10 GeV ||150|| , with the scale given by the Xb2 
mass. The measurement at Q ~ 80 GeV |[L12| , |113|| is based on hadronic Z° 



decays: 7r qqginci., in which ginci. refers to a gluon jet hemisphere recoiling 
against two almost coUinear quark jets q and q in the opposite hemisphere. 
For the results shown here, the ginci. hemisphere results in Ref. ||113|| have 
been multiplied by a factor of two both for the multiplicity and energy scales 
so that they correspond to gg "two-jet events" analogous to the T data. 
Although many studies of gluon jet multiplicity have been performed, only 



the data in Fig. |T6| are based on an unbiased definition of the jets. These 
data are summarized in Table 

The solid curve in Fig. Il6| shows the result of a fit of the 3NL0 expres- 
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Figure 16: The mean charged particle multiphcity of gg events from a color 
singlet point source versus energy scale Q. The solid curve shows a fit of 
the 3NL0 expression, eq. (Q), to the data, using n/=3, where Uf is the 
number of active quark flavors. The dashed curve shows the prediction of 
the HERWIG Monte Carlo for the mean charged particle multiplicity in gg 
events. 
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sion for gluon jet multiplicity, eq. (pi]), to the data, with ?/=ln(Q/A), where 
A is a fitted parameter related to the perturbative cutoff Qo by A=Qo/2 (see 
Section The other fitted parameter is the normalization constant K in 
eq. (|9^). For this fit, %=3. The T(1S) measurements |p.l4|| are not included 
in the fit because a systematic uncertainty was not provided for them: the 
smallness of the statistical uncertainty for these data would unduly bias the 
fit results. Nonetheless, these data are seen to lie near the fitted curve. The 
results for the fitted parameters, along with those found using n/=4 and 5, 
are given in the top portion of Table |^. The results of the fits with ^^/=4 and 



5 are virtually indistinguishable from that shown by the solid curve in Fig. 116 



The uncertainties given in Table |^ for the gluon jet parameters are defined 
by the maximum deviations observed when the gluon jet measurements are 
varied by their one standard deviation uncertainties. 

From Fig. and the results in Table ^ it is seen that the 3NL0 expression 
provides a good description of the growth of gluon jet multiplicity with energy 
using a physically sensible value of the coupling strength, corresponding to a 
value of A in the range from about 100 MeV to 1 GeV. For example the fitted 
resuh A = 0.64 ±0.17 GeV obtained for n/ = 5 yields ^^(Mz) = 0.142 ±0.008, 
compared to the world average value based on Af^ of ^^(Mz) ^0.120 ± 
0.002 ||80[. We remind the reader that A is not in general the same as Aj^. 

Experimental measurements of unbiased quark jet multiplicity are shown 
in Fig. |1^. These data are the inclusive charged particle multiplicity values 
of e^e~ hadronic annihilation events, corresponding to "two-jet events" as 
for the data of Fig. 16. The scale is Q=-E'c.m.- The results shown for the LEP 
experiments are combined values of ALEPH, DELPHI, L3 and OPAL. The 
combined values are obtained using the unconstrained averaging procedure 
described in , for which a common systematic uncertainty is defined by the 
unweighted mean of the systematic uncertainties quoted by the experiments. 
LEP-1 refers to data collected at the Z° peak, LEP-1.5 to data collected at 
-E-c.m.^lSS GeV, and LEP-2 to data collected at or above the threshold for 
W+W~ production. The quark jet data are summarized in Table 

The vertically striped band in Fig. O shows the 3NL0 prediction for 
quark jet multiplicity, eq. ( |104|) , for n/=3, using the values of A and K found 
from the fit to the gluon jet measurements (top portion of Table The 
width of the band corresponds to the uncertainty in the value of A presented 
in the top portion of Table ^ Almost identical results to those shown by 
the band are obtained using nf=A or 5. The band lies 15-20% below the 
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A (GeV) 


K 


X^/bins of data 


(a) Gluon jets, eq. (p4|) 


nf=3 


1.03 ±0.24 


0.288 ±0.037 


0.01/2 




0.84 ±0.21 


0.244 ± 0.034 


0.01/2 


nf=5 


0.64 ±0.17 


0.205 ±0.031 


0.01/2 


(b) Quark jets, eq. ([TO^ 


1), K=&xed 






0.670 ±0.036 


0.288 


26.3/18 




0.579 ± 0.034 


0.244 


31.3/18 


nf=5 


0.469 ±0.031 


0.205 


38.6/18 


(c) Quark jets, eq. (|T04|) 




0.262 ±0.072 


0.198 ±0.022 


6.6/18 




0.188 ±0.063 


0.156 ±0.020 


6.8/18 


n/=5 


0.119 ±0.048 


0.118 ±0.018 


7.1/18 



Table 5: Results of a fit of QCD expressions for the scale evolution of event 
multiplicity to the measured mean charged particle multiplicities of unbiased 
(a) gg (top), and (b) and (c) qq (center and bottom) events. For the gluon 
jets, the fitted expression is the 3NL0 result, eq. (^4]). For the quark jets, 
the expression is the 3NL0 result, eq. ( |104| ). In (b) the normalization K is 
set equal to the value found from the gluon jet fit; a one parameter fit of 
the scale parameter A is performed. In (c) a two parameter fit of A and K 
is performed. For the quark jets, the uncertainties are evaluated by varying 
the fit range as described in the text. 
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50 100 150 200 

Scale Q = ^ (GeV) 

Figure 17: The mean charged particle multiphcity of e"''e~ hadronic an- 
nihilation events versus energy scale Q=£^c.m.- The solid curve shows a fit 
of the 3NL0 expression, eq. (|104|) , to the data, using a one parameter fit 
of A (see text). The dotted curve shows the corresponding result from a two 
parameter fit of A and K. The hatched band shows the 3NL0 result if the 
parameter values used in the fit of gluon jet multiplicity are used. The width 
of the band corresponds to the uncertainty in the value of A. For the ana- 
lytic curves, ^/=3. The dashed curve shows the prediction of the HERWIG 
Monte Carlo for the mean charged particle multiplicity in e"'"e~ qq events. 
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Experiment 



ARGUS 1115 if 
JADE 



152 



TASSO [m 



TASSO 



HRS |T54| 
TASSO 

TASSO ig 

AMY 

LEP-1 ra 



LEP-1.5 [153, [15 
LEP-2 



m pm 1160] 



157, 159, 161 



T57I, IT621 



LEP-2 
LEP-2 
LEP-2 111, |164|, |16T^ 
LEP-2 
LEP-2 
LEP-2 |TB3 
LEP-2 



T65|, ITM Wm 



1651 , [161 H 



T65| , [164i |166^ 



g=Ee.m. (GeV) 



10.5 8.35 ±0.02 ±0.20 
12.0 8.4 ±0.3 ±0.6 
14.0 9.30 ±0.06 ±0.41 
22.0 11.3 ±0.08 ±0.46 
29.0 12.87 ±0.03 ±0.30 
34.8 13.59 ±0.02 ±0.46 

43.6 15.08 ±0.06 ±0.47 
57.0 17.19 ±0.07 ±0.48 
91.2 21.15 ±0.01 ±0.29 
133 24.03 ±0.20 ±0.53 
161 25.39 ±0.23 ±0.41 
172 26.38 ±0.26 ±0.56 
183 26.84 ±0.16 ±0.48 
189 27.04 ±0.11 ±0.34 
192 27.21 ±0.25 ±0.32 
196 27.47 ±0.18 ±0.33 
200 27.63 ±0.17 ±0.32 
202 27.70 ±0.25 ±0.38 



Table 6: The mean charged particle multiphcity, {rich.), measured in e"'"e~ 
annihilations at various cm. energies. The first uncertainty is statistical and 
the second is systematic. A more complete compilation for energies below 
91 GeV is given in |p.67|| . The results shown for LEP are combined values 
from ALEPH, DELPHI, L3 and OPAL (see text). 



data, representing an inadequacy of the 3NL0 calculations to simultaneously 
describe the gluon and quark jet measurements using precisely the same 
values of A and K. This problem can mostly be attributed to the analytic 
prediction for the multiplicity ratio r, discussed below in Section [11. 2| . 

It is also of interest to fit expression (|104|) to the quark jet data. For 
this, we use the normalization K found in the fit of the gluon jet multiplicity 
(compare eqs. (p^) and ( |104| )) and perform a one parameter fit of A to the 
quark jet data using ^/=3. The result is shown by the solid curve in Fig. 
The results found using ^/=4 or 5 are almost identical to that shown by 
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this curve. The fit is seen to provide a reasonable overall description of the 
measurements. The fitted values of A for nf=3, 4 and 5 are given in the 
central portion of Table ^. The uncertainties of the quark jet parameters in 
Table |^ are defined by the maximum difference between the results of the 
standard fit and those found by fitting only data between 29 and 202 GeV, 
between 10.5 and 161 GeV, or by excluding the LEP-1 data point. 

Comparing the parameter values obtained from the fits to the unbiased 
gluon and quark jet measurements (top and central portions of Table it 
is seen that the value of A from gluon jets is about 40% larger than that 
from quark jets. For nf=5, the quark jet result A=0.469 ± 0.031 GeV yields 
as{Mz) = 0.135 ± 0.002, to be compared to as ~ 0.14 ± 0.01 found from the 
gluon jet data as mentioned above. Thus although the gluon and quark jet 
data are not well described using exactly the same values of A and K, they 
both yield physically acceptable results for A, leading to similar values of the 
coupling strength. In this sense, the analytic description of multiplicity in 
single gluon and quark jets is quite consistent. 



For completeness, we include in Fig. O the result of a two parameter fit of 



A and K to the quark jet data, using nf=3. The result of this fit is shown by 



the dotted curve in Fig. 17. The corresponding parameter values are given in 
the bottom portion of Table |^, with systematic uncertainties defined as for the 
quark jet parameters in the central portion of that table. The results for A are 
seen to be about 3 times smaller than if the normalization factor K is required 
to be the same for quark and gluon jets as described above. For nf=5, the 
quark jet result A=0.119 ± 0.048 GeV yields as(Mz) =0.109 ± 0.005. 



The dashed curves in Figs. and [T^ show the prediction of the HERWIG 
Monte Carlo |j70| for the mean inclusive charged particle multiplicity of gg 
and e^e~ annihilation events, respectively. The prediction of HERWIG is 
generally similar to the fitted analytic results. 

11.2 Multiplicity ratio 

A test of the QCD prediction for the ratio r of mean multiplicities between 
gluon and quark jets is possible using the data in Table ^. Such a test is in 
principle limited to these data since they are the only unbiased measurements 



of gluon jet multiplicity currently available. See also Refs. |168|| , however, 
where results for r based on gluon jets from e^e~ three-jet events are also 
discussed. 
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Experiment 


Q (GeV) ^jet(GeV) r 


CLEO ||114|| 

CLEO 11150 , ARGUS |151| 
OPAL [|113|| 


5.5 2.75 1.04 ±0.02 ±0.04 
10.3 5.15 1.118 ±0.011 ±0.032 
80.2 40.1 1.514 ±0.019 ±0.034 



Table 7: Measurements of the multiplicity ratio r between unbiased gluon and 
quark jets. The first uncertainty is statistical and the second is systematic. 



To obtain a result for r from their measurements of gluon jet multiplicity 
in T(1S) decays, CLEO ||114|] divided the gluon jet multiplicity given in the 
top three rows of Table ^ by the charged particle multiplicity in radiative qq7 
events with similar hadronic recoil mass values. The three T(lS)-based data 
points are shown with their statistical uncertainties at ?/c~0.01 in Fig. ^ 
CLEO combined the three data points into a single value at an effective 
scale of about 5.5 GeV ||114||. The resulting value for r is given in the top 



row of Table |^. 

OPAL ||113|1 divided their result for the multiplicity in gluon jet hemi- 
spheres at 40 GeV by the multiplicity in hemispheres of light quark (uds 
flavored) events with about the same energy. The reason for selecting light 
quark jets is that it provides closer correspondence with the massless quark 
assumption of the calculations. The OPAL result for r is shown at Uc ~ 
3 X 10~^ in Fig. |^ and is listed in the bottom row of Table |^. 



CLEO did not report a result for r based on their measurement ||150 
of gluon jet multiplicity in T(3S) decays. The effective scale of this result, 
10.3 GeV, is not too different from the scale of the lowest energy quark jet 
result in Table |^, however, viz. 10.5 GeV. Dividing the 10.3 GeV gluon jet 
result from CLEO by the 10.5 GeV quark jet result from ARGUS yields the 
result for r given in the central row of Table |^. 



The experimental results for r are shown in Fig. |T8|. The results are pre- 
sented as a function of the jet energy E-^^t rather than the "two-jet event" 
scale Q to correspond to the analytic results. Fig. [T8| includes the QCD ana- 
lytic predictions for r at various orders of perturbation theory. The data are 
seen to lie substantially below the predictions (Section |6.1| ). The theoreti- 
cal predictions approach the experimental results more closely, however, as 
higher order terms are included. At the scale of the Z°, the 3NL0 expression. 
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Figure 18: Experimental results for the multiplicity ratio of unbiased gluon 
and quark jets, in comparison to QCD analytic predictions. The curve 
marked 3NL0(e) shows the analytic prediction at order 3NL0 with an ap- 
proximate accounting for truncation of the limits of integration at and 
1 — e~^, see Ref. 
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The dashed and dash-dotted curves show the pre- 
diction of the HERWIG Monte Carlo at the hadron and parton levels. The 
experimental and HERWIG results at the hadron level are based on charged 
particles only. 
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eq. (0), predicts r ^ 1.7, about 13% larger than the OPAL result Vexp. ~ 1.51. 
At the scale of the T, the difference between theory and data is larger. It 
is likely that uncalculated non-perturbative terms play an important role at 
this low scale, however. By itself, it is perhaps surprising to obtain any level 
of agreement between theory and data at all given the large value of the 
expansion parameter 70 ~ 0.5. 

The dashed and dash-dotted curves in Fig. 18 show the predictions of 
the HERWIG Monte Carlo for the ratio r at the hadron and parton levels. 
The HERWIG prediction is obtained by dividing the results for hemispheres 
of gg events by those for light quark qq events. The HERWIG results at 
the hadron level are based on charged particles only. HERWIG is seen to 
represent the OPAL measurement quite well at both the hadron and parton 
levels, implying a modest correction for hadronization at the TP energy. 

The computer solution of the QCD equations, discussed in Section [7], 
leads to near perfect agreement with the OPAL result for r, as shown in 
Fig. p. The computer solution differs from the analytic estimates by use of 
the explicit multiplicities in the integrals, with no Taylor series expansion, 
leading to a more accurate accounting of energy conservation, as mentioned 
in Section |^. This fact becomes especially important for higher order terms 
because of the large value of 70. The computer solution also implements 
the exact preasymptotic limits of integration, in contrast to the analytic 
solutions, as was also mentioned in Section 0. This agreement between the 
experimental and theoretical (computer) results for r at the scale of the TP 
is an impressive success for the QCD approach to multiplicity in jets. 

In contrast to the result at the scale of the Z°, the computer prediction 
for r at the scale of the T exceeds the data by about 25% (see Fig. |^). For 
such low scales, neglected non-perturbative effects are likely to be important, 
as noted above. 

Recently, results on the charged particle multiplicity of jets in dijet events 
produced in pp collisions have been reported by CDF Collaboration ||169| . 
The dijet masses range from 80 to 630 GeV. This substantially increases the 
range of accessible jet energies compared with e"'"e~ collisions. The multi- 
plicity ratio r has been estimated to be 1.7 ± 0.3 from these data. It should 
be noted that this result is obtained using a fit of an MLLA expression to 
the data. Since in MLLA quark and gluon jets differ only by the constant 
factor r=9/4, it is unclear how the experimental result r~1.7 should be 
interpreted. Also, the large uncertainty precludes a precise differentiation 
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between the different theoretical predictions in this region. 



11.3 Slopes 

According to eq. (p3D, the ratio r becomes smaller as the scale decreases 
since 70 increases. The same trend is observed experimentally (see Fig. ^ 
or |T8|). Thus the data and theory are in qualitative agreement. From Figs. ^ 
and |18| it is seen that the theoretical predictions decrease with scale more 



slowly than the data, however. For example, the computer solution for r [|78 
agrees perfectly with the data at the scale of the Z°, but it exceeds the 
experimental result by about 25% at the scale of the T (Fig. ^, as discussed 
in Section |11.2| . Thus there is a significant difference between the theoretical 
and experimental results for the slope r (see eq. (^5])), at least for energies 
below the Z°. From Fig. 2 in the slope predicted by the computer 
solution can be estimated to be about 0.096 at the Z". Using the data in 
Fig. ^ the corresponding result is about 0.174. Both these values differ from 
the analytic prediction of 0.06 given in Section |6.1| (cf. eq. (P^D). This implies 
a strong influence of higher order perturbative corrections, even at the TP, 
and of non-perturbative terms at lower energies. 

The ratio of slopes r«,eq. (||),is less sensitive to higher order corrections 
than r and is better approximated by the MLLA expression than r, as 
discussed in Section Experimental results for r^^-* are available from the 



DELPHI ||168|| and OPAL ||170|| Collaborations. These data are based on 
three-jet events selected using a jet finder and thus - unlike the other e^e~ 
data presented in Section |Tl] - do not employ the unbiased (hemisphere) 
definition of a jet utilized by the calculations. This makes a quantitative test 
of the QCD prediction uncertain. Nonetheless we proceed to a comparison 
of experiment with theory with an eye towards qualitative agreement. 

The energy scale of the experimental results ||168| , |1 70|| for r^^^ is presented 
in terms of the so-called jet hardness K=2i? sin(0/2) (see e.g. 0), where 
is the opening angle between the two lowest energy jets. For small B, 



y ^ hi{K/Qo). The jet hardness has been shown [16S] to be a more appropri- 
ate scale than the jet energy when comparing jets embedded in a three-jet 
environment. 

Both DELPHI and OPAL observe the experimental value of r^^p to de- 



pend very little on scale. The OPAL measurement r^^|j=2.27 ± 0.09 (stat.) 



±0.27(syst.) is about one standard deviation higher than the DELPHI result 
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1.97±0.10 (stat.). Note that the systematic uncertainty quoted by OPAL is 
mostly due to the jet selection bias. The experimental results agree with 
the MLLA prediction r^^^^^ ^ 2.01-2.03. The expressions and 



( |101D yield r^^^Q ^1.86-1.92 in the energy range from the T to the Z^, in 
agreement with DELPHI but about 1.3 standard deviations of the total un- 
certainty below OPAL. Fig. O displays the values of r^^^ in the MLLA and 



3NL0 approximations for nf=4: in comparison to the experimental limits 
from DELPHI. 

A recent study ||171|| attempts to correct for the jet bias introduced by 



the three-jet event selection by subtracting the biased quark jet component 
from the multiplicity of the qqg events. The result, r^^'' = 1.77±0.03 (stat.) is 
substantially smaller than the DELPHI and OPAL results presented above. 
There are a number theoretical and phenomenological assumptions which 
enter this analysis, however. 



From the DELPHI result | |168 | for r^^J , the measured values of pi^, ~ 



'"exp/S can be determined to range from 0.53 at the T to 0.75 at the Z". 
Its perturbative values according to eq. ( |101D vary from 0.85 to 0.9. Either 



the next order perturbative or else non-perturbative terms are needed ||172 
to reproduce the large difference di=r^^Jp — rea;p~ 2(1 — pi^exp) observed in 
experiment. The simplest Pade expression ( |102| ) improves the situation only 
slightly. 

11.4 Higher moments 

Measurements of higher moments of the multiplicity distribution in unbiased 
gluon and quark jets are presented in Ref. [|^. The data are collected at 
the Z°. The experimental results are given for factorial and cumulant mo- 
ments with ranks 2-5. These data are shown in Figs, ^(a) and ^(b) in 
comparison to the predictions of QCD Monte Carlo programs, and are listed 
in Table ^. The quark jet moments are observed to be larger than those of 
gluon jets in accordance with theoretical expectation (Section |^), i.e. event- 
to-event fluctuations are larger for quark jets than for gluon jets. Note that 
factorial moments of different rank are highly correlated with each other sta- 
tistically. In contrast, cumulant moments of different rank are not (see e.g. 
Ref. |15[)- Thus there is a high degree of bin-to-bin correlation in Fig. |20|(a) 



but not in Fig. |2^(b). It is worth noting that cumulants of ranks 2, 3 and 4 
are directly correlated to dispersion, skew and kurtosis, respectively, i.e. to 
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Figure 19: The ratio of slopes of the average multiphcities in gluon and quark 
jets, r^^\ in the MLLA and 3NL0 (with its Fade expression) approximations. 
The theoretical results are obtained using nf=4. Experimental limits (statis- 



tical uncertainties only) from Ref. ||168| are shown by the dash-dotted lines. 
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the ordinary moments, cf. Ref. |TH . 
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Figure 20: The (a) factorial (Fg) and (b) cumulant (Kg) moments of the 
charged particle multiplicity distributions of unbiased single gluon and quark 
jets, for 41.8 GeV ginci. gluon jets and 45.6 GeV uds quark jets |jT5[, as 
a function of the rank q. The predictions of the HERWIG and JETSET 
parton shower Monte Carlo event generators are also shown. 

The experimental results for 41.8 GeV gluon jets, =1.023, and for 
45.6 GeV uds quark jets, =1.082, are much smaller than the DLA predic- 
tions, viz. 1.33 and 1.75, respectively (see eq. ( |11U| )). The analytic results in 
MLLA, F^{MLLA) ^ 1.11 and Fi{MLLA) ^ 1.22 (for n/=4), agree much 
better with the data. This is analogous to the situation described in Sec- 



tion |TL2| for the multiplicity ratio r, i.e. the inclusion of higher orders results 
in a marked improvement in the description of the data compared to lead- 
ing order. If one accepts the effective value of as averaged over all the 
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Q 


Fg, gluon jets 


Fg, quark jets 


2 


1.023 ±0.008 ±0.011 


1.0820 ± 0.0006 ± 0.0046 


3 


1.071 ±0.026 ±0.034 


1.275 ±0.002 ±0.017 


4 


1.146 ±0.059 ±0.074 


1.637 ±0.005 ±0.042 


5 


1.25 ±0.11 ±0.13 


2.274 ±0.014 ±0.093 


Q 


Kg, gluon jets 


Kg, quark jets 


2 


0.0233 ±0.0083 ±0.0109 


0.0820 ±0.0006 ±0.0048 


3 


0.0010 ±0.0039 ±0.0048 


0.0291 ±0.0006 ±0.0035 


4 


0.0000 ±0.0023 ±0.0015 


0.0081 ±0.0007 ±0.0015 


5 


-0.0005 ±0.0018 ±0.0014 


-0.00300 ± 0.00096 ± 0.00095 



Table 8: The factorial (Fg) and cumulant {Kg) moments of the charged 
particle multiplicity distribution of unbiased 41.8 GeV gmci. gluon jets and 
45.6 GeV uds quark jets as a function of the rank q. The first uncertainty 
is statistical and the second is systematic. 

energies of the partons during the jet evolution to be as ~0.2, one obtains 
F^{MLLA)^ 1.039 and F^iMLLA) ^1.068, which are quite close to the 
experimental results. In this sense the MLLA prediction can be said to de- 
scribe the widths of the gluon and quark jet multiplicity distributions at the 
Z" energy to within 10% accuracy. 

Unfortunately the NNLO and 3NL0 terms worsen the agreement with 
data compared to MLLA (but not compared to DLA), as was already men- 
tioned in Section |6.2| . The analogy with r, for which each subsequent correc- 
tion brings the theory closer to data, does not hold here. The NNLO calcu- 
lation yields slightly larger values for F^ and F^ than in MLLA, resulting in 
a slight worsening of the description of experiment. The 3NL0 corrections 
are large and negative, yielding F2^(3A^LO) ^ 1.01 and F^{3NLO)^0M 
for nf=4:. These latter results are in qualitative conflict with the data in 
that they imply F2 < F^. The problem in the analytic description of higher 
moments can be traced to an inappropriate treatment of soft particles, as 
discussed in Section 

Given this problem with the analytic approach, it is remarkable that the 
computer solution of the QCD equations provides a near perfect description 
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of the higher moments. This is shown in Fig. which displays the com- 
puter solutions for the factorial moments of single gluon and quark jets [ 104 
in comparison to the corresponding data from Table |^. Overall, excellent 
agreement is observed between experiment and theory for ranks up to g ~ 5. 
This suggests that the failure of the analytic approach to describe the widths 
of single quark and gluon jets is mainly a technical issue. The success of the 
computer solution again emphasizes the importance of an exact treatment of 
energy conservation and the limits of integration. 

Further insight into this issue can be gained by examining the predictions 
of the ARIADNE Monte Carlo |7l| at the parton level, as was done in ||173| 
using Monte Carlo parameter values consistent with LPHD.[] At the parton 
level, ARIADNE is roughly equivalent to a DLA analytic calculation except 
that it includes exact conservation of both energy and momentum. The par- 
ton level predictions of ARIADNE for the factorial moments of uds quark 
jet hemispheres ||1 73|| are listed in the rightmost column of Table ^ These 
values are remarkably similar to the corresponding measurements in Table |[ 
This suggests that an exact treatment of energy-momentum conservation is 
essential for an accurate prediction of the higher moments. The correspond- 
ing ARIADNE results for hemispheres of gg events are also listed in Table ^ 
The agreement with the corresponding data in Table |^ is not as good as for 
quark jets. Nonetheless, the ARIADNE results are much more similar to the 
data than either the LO or 3NL0 analytic results discussed above, for both 
quark and gluon jets. 



Q 


Fq, gluon jets 


Fq, quark jets 


2 


1.044 


1.092 


3 


1.132 


1.310 


4 


1.271 


1.720 


5 


1.474 


2.444 



Table 9: The factorial {Fq) moments of the parton level multiplicity distribu- 
tions of hemispheres of gg and uds qq events generated using the ARIADNE 
multihadronic Monte Carlo event generator, as a function of the rank q. 



Qualitatively, the results shown in Fig. p| for the energy dependence of 



''Essentially this means using a small value of the cutoff parameter, Qo ~ 0.2 GeV |173| 
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Figure 21: The measured, normalized factorial moments of the charged parti 
cle multiplicity distributions of unbiased, separated gluon and quark jets [|l^ 



in comparison to the corresponding QCD predictions obtained from a com- 



puter solution of the equations for the generating functions ||104| . 
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the second factorial moments are similar to experimental trends reported at 
LEP [ 1471| . The NNLO terms tend to slightly violate the level of agreement 
with experiment while the agreement is even worse if the 3NL0 terms are 
included, analogous to the situation discussed above in connection with the 
data of Table 0. 



11.5 Oscillating Hq moments 



The Hq moments of quark jets, eq. (|T^), have been experimentally studied 
up to rank qK,11 in e"'"e~ annihilations | p.74| , |175| |. Results from the SLD 
Collaboration based on the charged particle multiplicity distribution from 
hadronic TP decays are shown in Fig. [ 174|| . Note that these results em- 
ploy full inclusive e"'"e~ hadronic annihilation events and thus correspond to 
unbiased "two-jet" qq configurations, unhke the data in Table || and Fig. |20| 
which are based on event hemispheres, i.e. single jets. 

The results in Fig. are strikingly similar to the QCD analytic predic- 
tions discussed in Sections O 



and 8.2 



VIZ. 



the Hq moments oscillate with 



increasing rank, with the occurrence of the first minimum at rank g^5. We 
again note the essential role that energy conservation plays for the theoret- 
ical prediction. Similar oscillations also have been observed in ep, pp, and 
11781, 



heavy ion collisions 



see Fig. |2^. From Fig. ^ it is seen that the 
amplitude of the oscillations increases as the structure of the colliding parti- 
cle becomes more complicated. The large positive values for ranks less than 
about 3 are due to strong correlations between particles from resonance de- 
cays. The positive values for q 7 or 8 imply the production of clusters (or 
mini-jets) at this multiplicity scale. The negative cumulants at q~ 5 suggest 
anti-correlations between resonances (and particles) inside these mini-jets. 
The behavior at larger ranks arises from a complicated mixture of attractive 
and repulsive forces inside higher multiplicity groups (jets). It is interesting 
that these features are rather universal in nature. 

In a study by the L3 Collaboration ||175|| , it is noted that Monte Carlo 
simulations provide a good representation of the oscillations in the data. 
This is found to be true even for simulations based on the MLLA approx- 
imation, despite the fact that analytic calculations predict oscillations only 
at order NNLO and above (Section p.2| ). This situation can easily be un- 
derstood to be a consequence of energy conservation. Unlike the MLLA 
analytic results, Monte Carlo results based on MLLA formulas incorporate 
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Figure 22: Measured ratio of cumulant to factorial moments, Hq, as a func- 
tion of the rank q, for the charged particle multiplicity distribution in e+e" 
hadronic Z° decays [|1 74\\ . The uncertainties are statistical only. 
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exact energy-momentum conservation. However, one of the principal differ- 
ences between analytic results at NNLO compared to lower orders is that 
energy conservation is included. Thus Monte Carlo results based on the 
MLLA effectively incorporate terms at NNLO and above. This provides 
an explanation of why MLLA based Monte Carlo event generators exhibit 
oscillations of Hq in agreement with the data. 



11.6 Soft particles 

So far, our discussion has focused exclusively on multiplicity in full phase 
space, i.e. within an entire jet. It is of considerable interest to also consider 
multiplicity within a limited region of phase space. In this section, we con- 
sider the hmited phase space volume defined by soft particles, i.e. particles 
with an energy or momentum below some cutoff. See also Section 14 where 
multiplicity in limited phase space is discussed in the context of intermit- 
tency. 

Particle multiplicity in limited phase space is not described by the ana- 
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Figure 24: Soft gluons radiated from a gluon (top) or quark (bottom) jet can- 
not resolve the color substructure of the jet because of their long wavelength. 
They therefore couple to gluon and quark jets with relative strengths given 
by the respective color factors Nc=3 and Ci?=4/3, see Refs. |113| , |176 . 
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lytic approach based on generating functions presented above in Sections 
Nonetheless soft particle multiplicity possesses a simple theoretical interpre- 
tation when considered in terms of the ratio between unbiased gluon and 
quark jets (see below). Recent theoretical and experimental progress in this 
area merits discussion in our review. 

In QCD, the ratio of mean multiplicities between gluon and quark jets, 
r, equals rQ=Nc/Cp=2.25 in the asymptotic limit Ej^t oo, as discussed in 
Section |6.1| . The approach to this limit is very slow, however. For finite en- 
ergies such as Ejet ~Mz, preasymptotic corrections from higher order terms 
and conservation laws limit r to values of approximately 1.5-1.7 (Sections ^]1| 
and 1^). Energy and momentum conservation only apply to full phase space, 
however. In limited phase space, conservation laws are not a constraint. Fur- 
thermore, the asymptotic condition E^^^t oo may effectively be replaced by 
the condition that the energy of the particles being considered, E, satisfy 
Qo < E « Ejct- By selecting soft particles, one therefore fulfils the asymp- 
totic condition in at least an approximate manner, opening the possibility of 
observing r— *• ro=2.25 even for finite energies. 



These ideas are developed in Refs. [^, p.76|] . Because of color coherence 
in the radiation of soft gluons, the low energy end of the particle spectrum 
is found to be nearly independent of the energy scale Ejet- Qualitatively, 
this arises because soft gluons have long wavelengths and cannot resolve 
individual partons within a jet. Instead, they couple to the jet with a strength 
proportional to the color charge of the parton which initiated it, given by A^c 
for gluon jets and Cp for quark jets, see Fig. ^ Thus for soft particles, the 
ratio of mean multiplicities between gluon and quark jets, r^o/t, is predicted 
to approach Vsoft Nc/Cp as E ^ Qq irrespective of the energy scale. 

Using their samples of unbiased gluon ("ginci.") and quark jets (see Sec- 
tion |1L1]), the OPAL Collaboration tested this prediction by measuring r as 



a function of particle momenta p |113|. The results are shown in Fig. 25. As 



p becomes small, r is seen to saturate at a value of approximately 1.8. This 
result is larger than the result r ?s 1.5 obtained in full phase space (Table |^). 
However it is considerably smaller than the theoretical prediction Vsoft ~ 2.25 
discussed in the previous paragraph. 

Insight into this difference between the measured and theoretical values 
of Tsoft can be obtained by examining the predictions of the HERWIG Monte 
Carlo. The predictions of HERWIG for r versus p at the hadron and parton 
levels are shown in Fig. |20. For small momenta, the hadron level result 
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Figure 25: The ratio of mean charged particle multiphcities between un- 
biased gluon and quark jets as a function of \np, with p particle momen- 



tum [|113|| . The results are shown in comparison to the predictions of QCD 
Monte Carlo programs. 
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Figure 26: Predictions of the HERWIG Monte Carlo fTg] for the ratio 
of mean particle multiplicities between unbiased gluon and quark jets as a 
function of particle momentum p, for partons and charged hadrons. 

(dashed curve) converges to Vgoft ~ 1-8, in agreement with OPAL (see also the 
Monte Carlo curves in Fig. ^5|) . The parton level prediction (solid histogram) 
converges to r^o/t^ 2.25, in agreement with the theoretical expectation [fT^j . 
Thus the difference between the measured and theoretical values of Vsoft can 
be attributed to the effects of hadronization. 

For soft particles, hadronization is most important within the core of a 
jet, corresponding to small transverse momenta p± with respect to the jet 
axis. For example. Fig. |2^ shows the predictions of HERWIG for the p± 
distributions of particles in quark jets, for partons (solid histogram) and 
charged hadrons (dashed curve). Partons are overwhelmingly produced with 
p± values larger than 0.8 GeV/c (because of the cutoff Qq in the Monte 
Carlo) in contrast to hadrons which predominantly appear at smaller p^. 
On the basis of this information, OPAL repeated their analysis selecting 
particles with p±_ > 0.8 GeV/c only, i.e. they restricted their analysis to the 
perturbative region of phase space as predicted by HERWIG. 



The results are shown in Fig. The charged particle multiplicity ratio 
r is shown function of the softness of the particles. The softness 

of particles is defined by the maximum particle momentum Pmax. considered 



103 



— 1 1 1 1 1 1 1 r- 



— 1 1 r- 



— 1 1 r- 



^ 0.8 

'S 

3 0.6 
U 

p 0.4 

u 

0.2 - 



. ' \ 
_/ \ 
-I \ 
-I \ 
-I \ 
-^ \ 
-I \ 



HERWIG 5.9 partons 
HERWIG 5.9 hadrons 



_l I I I I I I L 



- .- -j - I 



__, L_ 



J I L_ 



1 . 



2 3 

p , [GeV/c] 



Figure 27: Predictions of the HERWIG Monte Carlo |70| for the transverse 
momentum distributions of partons and charged hadrons in unbiased uds 
quark jets. The transverse momentum p± is defined with respect to the jet 
axis. The ordinate units are arbitrary. 
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Figure 28: Ratio of charged particle multiplicities between unbiased gluon 
and quark jets for particles with large transverse momenta to the jet axis 
defined by pj_ > 0.8 GeV/c. The results are shown as a function of the softness 
of the the particles, defined by the maximum particle momentum Pmax. used 
to determine r. 
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when determining r. Unlike the results in Fig. particles are required to 
have > 0.8 GeV/c as explained in the previous paragraph. The solid curve 
shows the prediction of HERWIG. 

With no explicit cut on Pmax. ("All momenta"), the multiplicity ratio is 
predicted to be about 1.8. As softer and softer particles are selected (pmax. is 
decreased), the HERWIG prediction approaches the QCD result rgojjRi 2.25 
in conformity with theory ||7^. Thus a measurement of r^o/t for particles 
with p±>0.8 GeV/c effectively yields a measurement of the color factor 
ratio Nc/Cf- OPAL results are shown for Pmax.=2 GeV/c and 4 GeV/c. 
The result using Pmax.=4 GeV/c is r=2.29 ± 0.017 (stat. + syst.) jn|] which 
provides one of the most accurate experimental determinations of Nc/Cp 
currently available. Note that unlike all other measurements of Nc/Cp, this 
result is not based on a fit of a QCD motivated expression - in which N^/Cf 
is extracted as a fitted parameter - but is the ratio of directly measured 
quantities. This OPAL result confirms one of the oldest, previously unver- 
ified predictions for multiplicities in gluon and quark jets, that the parti- 
cle multiplicity ratio between gluon and quark jets should equal the color 
factor difference (in an appropriate region of phase space), dating back to 



Refs. [|T05|, [TTTf 



12 Jet and sub jet multiplicities 

To this point, we have considered only the final result of well developed 
cascades, namely hadron multiplicities. It is also of interest to study the 
intermediate stages of the cascade evolution, i.e. the jet substructure. Re- 
calling the phrase about "whorls inside whorls inside whorls" from the field of 
turbulence, one can ask what is predicted by pQCD for the structure of "jets 
inside jets inside jets" in e^e~ events. The ordering of emissions by trans- 
verse momentum (kt) is at the heart of this issue. Thus one can ask about the 
number of jets or subjets (the "jets inside jets" ) as the angular (fc^) resolution 
scale is varied. With very low resolution one obtains two jets, corresponding 
to the condition eq. (^) imposed on the equations for the generating func- 
tions, eqs. (|57|) and (^8|), which requires that a single quark-antiquark pair 
be initially created in an e"'"e~ annihilation. A three-jet structure appears 
when a gluon with large transverse momentum is emitted by the quark or 
antiquark. Three-jet event production is suppressed by an additional factor 
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of as, which is small for large transferred momenta. Therefore, this process 
can be calculated perturbatively. The well known exact three-jet matrix ele- 
ment ||180[| can then be directly compared with the results of the generating 
function approach described above. At relatively low transferred momenta, 
the jet evolves to angular (or kt) ordered subjets. By increasing the resolu- 
tion, more and more subjets are observed. In the ideal case, each subjet is 
an individual final-state hadron in the limit of very small resolution scales. 
The resolution criteria must be chosen in a manner to provide infrared safe 
results. Then the perturbative expansion gives rise to finite answers which 
can be confronted to experimental data, so long as the hadronization stage 
does not significantly alter the situation as is assumed by local parton-hadron 
duality. 

Details of the resolution criteria are not too important if ratios are con- 
sidered. In particular, one can form the ratio of subjet multiplicities in three- 
and two-jet events. Naively, this ratio should equal the corresponding ratio 
of Casimir factors: 

^ 2Cf + Ca ^ 17 
(n^^) 2Cf 8' ^ ' 

However soft gluon coherence suppresses this ratio to a value below 1.5 (see 
Ref. [9C]), demonstrating the importance of correlations when dealing with 
interfering jets. 

The first proposal to eliminate collinear and mass singularities when re- 



solving subjets was formulated in Ref. [p.81|| . According to this proposal, one 



should consider the relative angles between any two partons 9ki and their 
fractional energies and xi which exceed some finite values: 6ki > S and 
Xk,i > e. The algorithm can be simplified by replacing these two conditions 
by a single one, as is accomplished by the JADE and Durham proposals. In 
the JADE algorithm ||138|| , the relative invariant mass of two partons M^i/s^^'^ 
is required to exceed some value {MlJ s> ^c)- In the Durham (or kt) algo- 



rithm ||136| , |137|| , a similar inequality is imposed on the relative transverse 



momentum, or, more precisely, 

iki = 2(1 - cos 9ki) min ™^ > (145) 



where Ek,i are the energies of partons k and /. For small relative angles 6ki, 
the parameter C,c approximates the rescaled transverse momentum k^j^JQ'^ 
of the jet of lower momentum with respect to that of higher momentum. 
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Figure 29: The measured ratio of subjet multiplicities between three- and 
two-jet events as a function of the subjet resolution scale ||184|| , in com- 
parison to analytic calculations [^]. The results are obtained from e+e~ Z° 
data using the Durham jet finder with an initial clustering scale ?/i=0.010. 
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The Durham algorithm is well suited for analytic calculations and gives 
rise to small hadronization corrections according to Monte Carlo schemes. 
Therefore it is widely used for comparison of theory with experiment. The 
two limiting cases of ^ 1 and correspond to low resolution, for 
which the minimum number of jets are resolved|, and to the final parton (or 
particle) limit of the highest available resolution, discussed in detail in the 
previous sections. In terms of virtualities, ic=Ql/Q'^- Thus for Qc=Q one 
deals with the initial stage of the process, while ever smaller values of Qc 
correspond to its later subjet stages, and final hadrons appear at values of 
Qc approaching the non-perturbative scale Qo of a few hundred MeV. The 



number of jets increases as decreases, as calculated in Ref. 



Jet multiplicities in e+e processes are usually fitted by the formula 

«e-(Qc,Q)) = 2(np(^,^)), (146) 

whereas the multiplicity of charged hadrons has been approximated by a 
formula with two additional constants: 

{nX-{Qo,Q)) = 2Kch{nF ( ^^,^)) + constant (147) 



It is interesting to determine if eqs. ( [57| ) and ( pH]) can describe jet, subjet 
and hadron multiplicities in e"'"e~ annihilations in a unified manner. The 
answer is yes ||104|, [f8[| , so long as the equations are solved numerically using 



the proper limits of integration. By proper limits of integration, we mean a 
lower limit Xc and an upper limit 1 — Xc instead of and 1, with Xc=QcV^/Q 
as has been suggested in Refs. |p.04] , [f^ . The numerical solutions presented in 



Refs. ||104| , [78[] describe the average multiplicities of jets, subjets and hadrons 
using a common normalization corresponding to "one parton=one hadron" 
without additional constants, over a wide energy interval from the threshold 
region to the Z° energy. This result is especially remarkable given that only 
two hard partons are present at the initial stage, whereas up to 60 charged 
particles appear in the final stage at the Z° energy, some of which originate 
from the decays of intermediate resonances. The most common approach 
consists of inserting a so-called K-iactor, Kch (see eq. ( |147[ )), i.e. some con- 
stant relating the parton and hadron distributions which usually differs from 



,(iet) 



(^c=l)) corresponds to two initial quark jets in e^'^e collisions. 
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unity (typically, i^c/i^ 1-2). This difference from unity is ascribed to the 
confinement process. In Monte Carlo models, this stage is fitted by some 
phenomenological formulas which are not yet explained in the framework of 
the more theoretical approaches discussed above. 



The ratio of hadron multiplicities in gluon and quark jets is found ^8 



to be even smaller than in the analytic solutions of the evolution equations, 
bringing it into closer agreement with experiment, see Fig. ^ and the dis- 
cussion in Section |11.2| . The running property of the coupling constant is 
a crucial consideration for this result. Even more astonishing is that high 
moments of the multiplicity distributions and their oscillations fit the exper- 



imental data using the same normalization as used for r ||104|| , see Figs. |T0 
and This implies that power corrections due to energy conservation, ob- 
tained by using the correct integration limits, are essential, raising questions 
about the potential importance of other possible power corrections, due to 
instantons, etc., as discussed above. 

We have briefly described the latest theoretical developments in the prob- 
lem, leaving aside the corresponding formulas which can be found in Refs. Jl^ 

s n, m ^ m, 11, m. 

Subjet multiplicities have been studied experimentally by a number of 
groups [p.83|| - [p!85[] . Measured results for the ratio (|144|) as a function of the 
subjet resolution scale yo=^c are shown in Fig. ^ These results are obtained 
using the Durham jet finder with an initial clustering scale ^c=?/i=0.010 to 
select the two- and three-jet events. The ratio has a value of 1.5 for yo=yi, 
as required from its construction, and decreases rapidly at smaller resolution 
scales yo- The results are far below the naive prediction of 17/8 for all values 
of yo- Fig. ^ also shows the perturbative predictions of Ref. |Q, both for 
fixed order 0{as) calculations and including resummed terms in the MLLA 
approximation. The MLLA results are shown for different choices of the 
scale parameter A. Both the fixed order and MLLA results reproduce the 
qualitative behavior of the data for yo ^ yi, i-e. they exhibit a falloff from 
the value of 1.5 with about the same slope as the data as yo decreases. 

Subjet multiplicities have also been studied in separated gluon and quark 
jets, both experimentally ||185|| and theoretically ||186|| . An example is shown 
in Fig. ^ The jets in this example are defined using the Durham algorithm. 
The analytic results [|186|| are seen to represent the data fairly well for large 
values of the subjet resolution scale yo- 
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Figure 30: The subjet multiplicities of separated (a) gluon and (b) quark 
jets ||185|| in comparison to analytic results from Ref. ||186|| and to Monte 



Carlo predictions. 

13 Multiplicity of three-jet events 

The multiplicity of two-jet events was defined from the product of the gen- 
erating functions for single jets, eq. (^), implying their independence. The 
assumption of independence is valid because the angular separation of the jets 
is large, i.e. about 180°. For three-jet events, various angular combinations 
are possible. Mutual interference between the jets cannot be neglected. Thus 
the particle multiplicity in three-jet events depends on the angular topology 
of the events. 

To describe this effect it is necessary to go beyond the equations presented 
in this work, to consider angular correlations between jets in the framework of 
the equations for the generating functionals or directly in terms of Feynman 
graphs. Therefore we refer to Refs. |T^, |139|| where this problem is treated 



in detail. Up to further correction terms of order 0{as), the multiplicity of 
three-jet events can be approximately written as the sum of the multiplicities 
of independent jets at properly defined energy scales: 

{ripp^) ^ 2{nF{yF)) + (ncd/c)), (148) 

where yF=^'^{E* / h) and yG=ln(pt/2A), with E* the quark energy in the 
quark-antiquark cms and pt the gluon transverse momentum in this system. 



Ill 



T I I I I 1 1 I I I 1 I 1 1 I I I 1 1 I I 1 r 




15 



20 



40 



60 



80 



100 



120 



Figure 31: Measurements ||115|| of the mean charged particle multiphcity of 
three-jet Y events |[1871| as a function of the opening angle 6i between the 
two lower energy jets, for i?c.m.=91 GeV. The events are selected using the 
k± jet finder. The solid curve shows the result of a one parameter fit 
eq. (|149|) to the data within the fit region shown. 
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This expression has the correct hmit when pt^O and the two jet configura- 
tion is restored. 

Eq. ( |148| ) can be re-expressed using directly observable quantities such 
as the multiplicity of e^e~ events and the ratio r of multiplicities between 
gluon and quark jets: 



(^ffg) ~ K+e-(2^*)) + ^r{pt){n^+e-{Pt))- 



(149) 



These formulas also can be written in forms suitable for studies of gluon jets 
or to relate multiplicities at different scales (see Ref. [1T7|). 

Expression (|149|) has been tested in several experimental studies ||115| , 
|188|| . As an example, the data in Fig. ^ show the charged particle multi- 
plicity of so-called "Y events"^ namely three-jet events for which the angle 
between the highest energy jet and each of the two other jets is about the 
same, versus the opening angle 6i between the two lower energy jets. In 
Ref. ||188|| it is found that eq. ( |149|) describes this data accurately only if the 
angular separation between the jets is about 80° or larger, as shown by the 
solid curve in Fig. In Ref. |p.89|| , this effect is explained to arise as a conse- 
quence of the bias in quark jet multiplicity introduced by the three-jet event 
selection procedure. A modified version of eq. ( |149| ) is suggested [[L89|| based 
on the dipole model of QCD (see Ref. ||19CI|| and references therein), which 
implicitly includes non-perturbative string effects. The modified prediction 
is found to yield a better description of the experimental measurements, see 
Ref. im| . 



14 Evolution of distributions with decreasing 
phase space: intermittency and fractahty 

Multiplicity distributions can be studied not only in total phase space (as 
discussed in the previous sections for very large phase space volumes) but 
within any subset of it. For a homogeneous distribution of particles, the av- 
erage multiplicity decreases in proportion to the considered volume, whereas 
the fiuctuations increase. The most interesting problem here is the law gov- 
erning the growth of the fiuctuations and its possible deviation from a purely 



events were first studied in |187|| 
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statistical law related to the decrease of the average multiplicity. Such a de- 
viation necessarily would be a consequence of the dynamics of the underlying 
interaction. In particular, it has been proposed | |191|| to search for a power 
law behavior of the factorial moments in small rapidity intervals 5y: 

F,^{5yr^^'i^ {5y^Q\ (150) 

where 0(g) > 0. This proposal is motivated by an analogy to turbulence in 
hydrodynamics, where a similar behavior is known as intermittency and the 
(f){q) are called intermittency exponents. 

The power-like dependence in eq. (|150|) can be related to fractal properties 
of particle distributions in the available phase space, as was first discussed in 
Ref. ||192|| . Earlier, the fractal properties of branching processes and, in par- 



ticular, of QCD jets in e+e~ annihilation were considered in Refs. [|193|] - |195|| . 
From the point of view of multiplicity, intermittency implies a rather strong 
increase of the distribution's width, and a longer tail at high multiplicities. 

Experimental data from various processes over a wide energy range re- 
veal a power law dependence in agreement with eq. (|150D and thus suggest 
the relevance of intermittency to high energy multiparticle processes. Simple 
branching models have been proposed to explain this phenomenon ||196| - 



r^]. Calculations based on perturbative QCD provide more specific predic- 



tions than these models, however. The current state of affairs is reviewed in 



Ref. |T6|. In the following we show how Quantum Chromo dynamics produces 



intermittency |^4[, [P0(]|| - p0^ and then compare the theoretical predictions 



to recent experimental results. 

We again stress that QCD deals with partons (quarks and gluons), in con- 
trast to experiment which provides results based on hadrons. Local parton- 
hadron duality implies the proportionality of inclusive distributions. The 
validity of LPHD is not obvious for correlations, however, and is sometimes 
not fulfilled in Monte Carlo schemes. For phenomena such as intermittency, 
one can therefore anticipate a qualitative agreement between theory and data, 
at best. 

In contrast to the previous sections, here we rely on the diagrammatic 
approach rather than on the equations for the generating functions. This 
is because we must now deal with a small part of the parton content of a 
well developed jet, namely with those partons which are present in a small 
phase space volume. The pre-history of the jet as a whole is significant for 
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the subjet under consideration. This is most readily addressed using the 
diagrammatic technique, viz.: 

1. the primary quark emits a hard gluon with energy E in the direction 
of the angular interval 6, but not necessarily hitting the window; 

2. the emitted gluon produces a jet of partons with parton sphtting angles 
larger than the window size; 

3. among those partons there exists a parton with energy k which initiates 
a subjet which hits the window; 

4. all decay products of the subjet exactly cover the angle 6. 

This picture dictates the rules of calculation of the gth correlator of the whole 
jet. The gth correlator of the subjet AN^'^^kO) should be averaged over all 
possible ways it can be produced, i.e. convoluted with the inclusive spectra 
of partons in the whole jet and with the probability of creation of the jet 
(asKp). Analytically, this is represented by 

'j) « r f g-? r t-' -) 

(151) 

where AN^'^^ = Fq{nY is the unnormalized factorial moment (on the left hand 
side for the whole jet, and on the right hand side for the parton subjet with 
momentum k within the angle 9). Since the unnormalized moments increase 
with energy, whereas the parton energy spectrum decreases, the product 
D^AN^i\ke) has a maximum at some energy, and the integral over momenta 
can be calculated using the method of steepest descent. Leaving aside the 
details (see Ref. 0]), we describe the general structure of the correlator in 
the case of a fixed coupling constant, 7o=constant: 

if) \ ^ ^ Ff)\ 

AN^'^^ ocAQijj i—j , {fi = constant), (152) 

where the three factors represent the phase space, the energy spectrum factor, 
and the gth power of the average multiplicity. To obtain the normalized 
moment, expression (|152|) should be divided by the gth power of that part 
of the mean multiplicity of the whole jet which appears inside the window 6, 
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i.e. by the proportion of the total average multiphcity corresponding to the 
phase space volume Afl: 



AN{9) AnAN{eo). (153) 

If the analysis is performed in D-dimensional space, the phase space volume 
obeys 

An ~ 9^, (154) 

where 6 corresponds to the minimum linear size of the £)-dimensional window 
which stems from the singular behavior of parton propagators in Quantum 
Chromodynamics (see Ref. |0). This allows the factorial moments to be 
represented as the product of a purely kinematic factor depending on the 
dimension of the analyzed space, and of a dynamic factor not related to that 
dimension which depends on the coupling constant, i.e. 



~ r^(^-i)0^^°. (155) 

For small angular windows, 6 ~ 6y. The intermittency indices defined by 
eq. ( |15CI| ) are then 

0(g)=D(g-l)-^-^7o. (156) 

q 

Expression ( |156|) is valid for moderately small windows, for which the con- 
dition a^ln^o/^ < 1 is fulfilled. For extremely small windows it is necessary 
to account for the running property of the coupling constant. In this case 
the constant 70 should be replaced by an effective value (7), which depends 
logarithmically on the width of the window 6 and may be approximated 
by H 

+ 1 
4g2 



(7)^70(1 + ^6), (157) 



where 

ln(E6'o//i) 

Thus intermittency indices for very small windows are markedly smaller than 
for intermediate sized windows (fixed coupling regime), especially for low 
rank moments. Moreover, the simple power law, eq. (|150|) , becomes modified 



by logarithmic corrections, and the intermittency indices depend on the size 



116 



of the chosen interval. The resuhing curve of InFq versus — In 6' is charac- 
terized by two distinct regions with different slopes. A rather steep linear 
increase is predicted for moderately small windows 6, with a slope given by 
eq. ( |156| ). For smaller window sizes, this behavior is replaced by a much 
slower quasi-linear increase given by eqs. (|157|) and ( |158| ). It is easy to de- 
termine the location of the transition point between the two regimes and 
to show that this point shifts to smaller window sizes at larger values of q. 
In any case, the factorial moments of any rank increase as the size of the 
interval decreases. This demonstrates that the fluctuations of the multiplic- 
ity distribution become stronger for smaller intervals and, more importantly, 
that they significantly exceed Poisson fluctuations. 

Above, we have described the results of the double-logarithmic approxi- 
mation, with corrections from the running coupling constant, eq. ( |1571 ). Cor- 
rections from the modified-leading logarithm approximation (see eq. ( |164| ) 
below) are comparatively small, i.e. about 10%. The MLLA terms move 
the transition point between the power law and quasi-power law regimes to 
slightly smaller windows for all moments except the second one, which moves 
to somewhat larger windows. This tendency can be attributed to the mutual 
influence of the energy spectrum and the average multiplicity. 

Of greater interest are qualitative effects introduced in MLLA which yield 
a functional dependence on the rank q from terms proportional to q'j. Con- 
sider an analogy to statistical mechanics, presented in Ref. |Q, where the 
quantity 

$=1-M±1 (159) 

q 

is interpreted as the free energy and the rank q as the inverse temperature 
(3 = 1/kT. In the lowest order approximation, $ increases monotonically 
with q. Higher order corrections produce a maximum in $ at values of q 
where Hq exhibits its first minimum (see eq. (0)), i.e. at 

^ + ^ ~ 5. (160) 



hi'jo 2 

In statistical mechanics, this corresponds to zero entropy, i.e. to a phase 
transition. In Quantum Chromodynamics it just indicates the role of the new 
parameter q^y. This feature is related to the singularities of the generating 
function and to the pinching behavior of the zeros of the truncated generating 



function discussed in Section 10 
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The above results can be restated in terms of fractals. The power-like 
behavior of factorial moments suggests fractal properties of particle (parton) 
distributions in phase space. According to the general theory of fractals (see 
Ref. |16| and references therein), intermittency indices are related to fractal 
(Renyi) dimensions Dq by the formula 



<P{q) = {q-l){D-Dq). 



(161) 



In the double-logarithmic approximation taking into account eqs. (|156|) and ( p.571 ), 
one obtains 

-(7) = (70 + —)(! + ^-^e)- (162) 



q ' ' " q 4g2 

The first term in the first bracket corresponds to mono-fractal behavior in 
the case of fixed coupling and is due to the increase in average multiplicity. 
The second term in the first bracket corresponds to multi-fractal properties 
and is related to the falloff of the energy spectrum. One can easily obtain 
the multi-fractal spectral function in this case (see Ref. ||9^). It is found 
that the fractality in Quantum Chromodynamics has a purely dynamical 
origin {Dq ~ 70) related to the cascade nature of the process, and that the 
kinematic factor in relation ( p.55| ) has an integer dimension. The term with 
e in eq. ( |162|) accounts for the running property of the coupling constant. 
This term slightly violates pure multi-fractal behavior, giving rise to a slower 
increase of factorial moments at small angles. The fixed coupling regime 
corresponds to e=0, yielding 



7o- 



q + l 



(163) 



Further corrections to eq. ( |162|) appear in MLLA |^|. The general form of 
eq. ( [162| ) remains valid, but 70 (Q) is replaced by an effective 'Jo^^ {Q) which 
depends on q: 



Yo"{Q)=io{Q)+im 



j3o_ 
4iV. 



where 



— + 



/ 



+ 



1 



2{q+l) 2(g + l)(g2 + l) 4 



(164) 



(165) 
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For large ranks q one obtains a negative shift: 



This corresponds to the general tendency of negative shifts from MLLA cor- 
rections. However, this represents only a partial accounting for higher order 
effects. Other corrections stemming from energy conservation are more dif- 
ficult to include and are not incorporated into the above equations. 

An alternative approximation for the e dependence of Dg in DLA has 
been suggested in Ref. |203|: 



D, = 27o(g)^- ^"^ . (167) 



A somewhat different expression for Dg was derived in Refs. ||201| , |202|| start 



ing from the expression for cumulant moments, supposing that they converge 
to factorial moments at high energies up to a factor of as discussed above: 

Dg = 270(g) ^ ■ ^ (168) 

q — I e 

The variety of proposed forms for Dg suggests that there are neglected 
effects which are important for a proper analytic treatment. In particular, en- 
ergy conservation constraints are just as severe in small phase space windows 



as in total phase space (see [205 1). This problem has not yet been solved. 
Thus the above formulas should be compared to experiment on a qualitative 
level only. In addition, variables such as rapidity or pseudo-rapidity should 
be used. Other variables, such as the invariant mass Q^, are more suitable for 
studies of Bose-Einstein correlations than of intermittency (for more details, 
see M )- 

In general, the qualitative trends predicted by the QCD equations are 
observed in experiment. Intermittency is observed to be stronger in e"'"e~ an- 
nihilations than in reactions with initial hadrons, and weaker yet in nucleus- 
initiated processes [0. This is related to the increasing number of com- 
peting subprocesses (sources) in reactions with complicated structure, which 
smoothes out the intermittent behavior of individual sources. Due to the 
stronger effect and higher statistics, intermittency studies in e^e~ collisions 
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are more conclusive. Recent studies from LEP p07] , p08|| represent con- 



siderable improvements upon earlier results (summarized in Ref. [16|). In the 



new studies, the experimental data are compared with QCD predictions for 
local multiplicity fluctuations in 1-, 2- and 3-dimensional intervals at LEP-1 
and LEP-2 energies. The general predictions of QCD are confirmed by the 
data. The factorial moments rise in an approximately linear manner for large 
bins and level off for smaller ones. The dependence on the bin dimension, the 
order of the factorial moments, and on energy are qualitatively reproduced, 
i.e. the moments increase as the bin dimension increases, become larger for 
higher orders, and increase with energy. The analytic predictions depend 
sensitively on the QCD parameter Qq, as is also true for the predictions of 
the QCD Monte Carlo models. 

On a quantitative level, significant differences are observed between the- 
ory and data, however. The factorial moments increase faster in experiment 
than in the calculations as the phase space interval decreases, and gener- 
ally level off at larger bin sizes than the theory as the intervals continue to 
decrease. An example is shown in Fig. ^ ||207|| , for which z oc ln(0o/0) mea- 
sures the size of the phase space window < Oq. Thus z=0 corresponds to 
the maximum phase space region 0=0o (where 0o=25° in this example), 
while z increases as the window size decreases. Variation of the QCD param- 
eters, e.g. Qo, does not improve the level of agreement between theory and 
data. 

Given the quantitative discrepancies between theory and data, an eval- 
uation of QCD parameters using intermittency is not possible at present. 
The discrepancies are not surprising given the crucial role played by energy 
conservation in total phase space, see Sections and the absence of en- 
ergy conservation from the analytic description here. Energy conservation 
is included in Monte Carlo studies, and they provide a reasonably good de- 



scription of the data as is shown in Refs. |66| , |207| , |208 |. 

The increase of correlations for smaller bin sizes, as described by the 
behavior of the cumulants, reveals the important role of genuine four- and 



five-particle correlations in e"'"e~ annihilations [Q. This is in contrast to 
nucleus collisions where it has not yet proved possible to observe correlations 
for groups of more than two particles. The fractal nature of the particle 
distributions in e"'"e~ events has been also demonstrated ||209|| using the so- 



called bunching parameters [plOj , |211 



The fractality of particle distributions in limited phase space volumes 
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Figure 32: Measurements of the ratio Fq{z)/Fq{{)) from e'^e" Z° data for 
g=2-5, as a function of the scahng variable z = t ||207|| . The data are 
compared to QCD analytic calculations with fixed coupling strength as 
(eq. (IID): DLA (a) (eq. (|T62|)), DLA (b) (eq. {^)), DLA (c) (eq. (|168D), 
and MLLA (eq. ([164D ). Note that DLA is referred to as "DLLA" in the figure 
legend. 
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suggests fractal properties for colliding objects in ordinary space. Surely, 
owing to its dynamical origin, such a structure would itself be dynamical, 
i.e. rapidly evolving in space and time. There are two reasons to favor this 
interpretation. First, the evolution of the parton shower in ordinary space 
must give rise to a tree-like structure of a fractal type which should evolve 
due to the cascade evolution. Second, lattice computations in SU(2) gluo- 
dynamics ||212 have shown that the system of gluons in the vicinity of the 
phase transition is fractal in the sense that it fills a volume V bounded by 
a surface S which are related by the formula V ~ S^-^"^. This is typical of 
fractal objects: the exponent would equal 1.5 for ordinary three-dimensional 
objects. 

The geometric fractality of macroscopic bodies has been revealed by mea- 
surements of the power-like shape of their structure functions when point-like 
particles (photons, electrons, neutrons, etc.) are scattered from them. Using 
this example, one might try to measure ||213|| structure functions in deep in- 
elastic processes, in order to determine the fractal dimensions of the scattered 
partners. Theoretically, this has been considered only at the level of models. 
Experimental difficulties have prevented direct tests of these models. The 
question of the fractal geometry of particles in ordinary space-time therefore 
remains open. 



15 Heavy quark jets 

The multiplicity of heavy quark jets is of special interest because it exhibits 
two peculiar effects with analogues in QED. It is well known from QED 
that photon emission from muons is strongly suppressed compared to emis- 
sion from electrons (as the ratio of corresponding masses squared). This 
suppression is related to the different masses in the propagators and to the 
vector nature of the emitted photons. The second effect, known from cos- 
mic ray studies and named after Chudakov, is that two opposite charges in 
an electron-positron pair screen each other if they are not well separated, 
i.e. until the distance between them becomes larger than the inverse trans- 
verse momentum of the resolving quanta. 

In QCD, the emission of gluons by heavy quarks is suppressed in a similar 
manner. At relatively small angles (^ -C 1), gluons are emitted according to 
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the formula 

e 



{e^ + elY' ^^^^^ 

where 6m = M/Eq, with M and Eq the mass and energy of a heavy quark. 
The form of the numerator in eq. ( [169| ) is due to the vector nature of the 
radiated quanta. The denominator equals the quark propagator squared. 
The larger the value of M, the stronger the suppression of forward {6 < 6m) 
emissions. This gives rise to "ring-like" events as advocated in Refs. ||214|| - 
[pl6|| using somewhat different but similar arguments, or to the "dead cone" 
phenomenon ||217| , |218|| . Note that the suppression of gluon emission at small 
angles does not necessarily imply a suppression of emissions with low trans- 
verse momenta. The origin of the ring-like events is the large mass of heavy 
quarks, which is not the same as kt ordering. 

There is a difference between QED and QCD with respect to the prop- 
erties of the emitted quanta, however. While photons do not have electric 
charge, gluons carry color. Radiated photons do not multiplicate, in contrast 
to gluons which produce jets. Gluons emitted at rather large angles {6 > 6m) 
multiplicate in the usual manner. For gluons emitted at smaller angles, it is 
necessary to account for interference between the secondary emissions and 
possible emission by the quark. This situation is analogous to the Chudakov 
effect. In QCD, the analogous effect for partons moving close to one another 
is called color transparency. 

However, there is a substantial difference for heavy quarks (for a more 
detailed discussion, see Ref. Because its mass is large, the heavy quark 

moves more slowly than the emitted gluon, and the two become separated 
more quickly than is the case when the quark is light. Therefore color charge 
screening is not as effective for heavy quarks as it is for light quarks and 
interference effects are less important. It can be shown that the primary 
gluon can be treated as an incoherent state with respect to the quark, and 
that, as a consequence, gluon emission by massive quarks is suppressed for 
angles 6 < 6m compared to light quarks. 

To obtain an expression for the multiplicity of partons accompanying a 
heavy quark pair in e"'"e~ annihilations, the distribution of gluons emitted by 
the quark, eq. ( |169| ), should be convoluted with the multiplicity of a gluon 
subjet. According to the above discussion, this expression can be integrated 
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over all transverse momenta exceeding a fixed value Qq to obtain 



<"«5(^'')) = 2 /„, ^ (l — j i ^A-F(x)^(„a(fe)>. (170) 

with K^{x) given by eq. fg), h = xk, and k'^ = W^[sm^{e/2) + ^^/4]. 
Eq. ( |170| ) yields the relation between the multiplicities accompanying quark- 
antiquark pairs with heavy and light components produced in e"'"e~ pro- 
cesses |pi7|| : 

Kgm) = MW)) - MV~eM)), (171) 

which is valid up to terms of order 0{as). Therefore the difference between 
the total multiplicities of events with heavy and light quarks is 

KU^)) - Kf(W^)) = - «(^/^M)), (172) 

where (n*^^) denotes the decay multiplicity of a heavy quark with mass M. 
The difference does not depend on the primary energy W but only on the 
heavy quark mass. 

In contrast to eq. ( |172| ), the so-called naive model ||219|| - p2T| considers 



the difference between heavy and light quark event multiplicities to be the 
consequence of energy rescaling only: 

{nQQiW)) = Mil-{xQ))W]), (173) 

where (xq) denotes the mean energy fraction of the heavy quark, (xq) = 
2{Eq)/W. In this case the difference in total multiplicities between heavy 
and light quark events is energy dependent, decreasing towards zero as W 
becomes larger. 

Measurements of the difference Sqe = {ng^iW)) — {n';^{W)) between 
heavy and light quark event multiplicities in e"'"e~ annihilations are shown 



as a function of the center-of-mass energy in Fig. ^ ||222|| . The heavy and 
light quark jet samples consist of b and uds flavored events, respectively. 
The data support the QCD hypothesis that Sqe is independent of energy, 
eq. ( |172| ), over the energy dependent hypothesis of the naive model. Note 
the essential contribution that the highest energy points (from LEP-2) pro- 



vide for this conclusion. The initial QCD estimates ||217|| predicted a rather 



large value 5q£~5.5. Later this prediction was reduced to (5q^~3.7 |p21 
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Figure 33: Measurements of the difference 5qi between the total event mul- 
tiphcities of heavy and light quark events in e^e~ annihilations, versus the 
center-of-mass energy, compared to a prediction motivated by MLLA QCD 
and the prediction of the so-called naive model of energy rescaling (the naive 
model is referred to as "Flavour independent" in the figure legend). The 



figure is taken from Ref. p22 . 
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The horizontal band in Fig. ^ at 5q^~3 [p22|| is not, strictly speaking, a 
theoretical estimate, but an admixture of the the MLLA statement that 6qi 
is constant together with the experimental measurement from the Z° pole. 

We will also briefly mention an interesting feature of heavy quark frag- 
mentation. Because of their large masses, heavy quark hadrons tend to retain 
a substantial portion of the primary quark energy, which can be described 
in the framework of perturbative QCD (for more details, see Ref. [1^1). The 
comparatively long lifetime of heavy quark hadrons plays an important role 
at extremely high energies in the spatial multiplicity distribution of cosmic 
ray showers |p23|| . Because of relativistic effects and smaller interaction cross 
sections, high energy heavy hadrons created in cosmic ray showers penetrate 
more deeply into detectors on average than light quark hadrons, giving rise to 
so-called long flying cascades. The multiplicity of the main cascade in these 
events is somewhat reduced, while particles appear at longer distances from 
the primary interaction point, compared to cascades without heavy quarks. 
Thus the mean length of the cascade is larger than normal, from which term 
"long flying" is derived. 



16 Conclusions 

The distribution of particle multiplicity in high energy inelastic processes 
provides essential information on the dynamics of strong interactions. Its 
importance stems from the fact that it contains all the correlations between 
final-state hadrons in an integrated form. It is comparatively easy to mea- 
sure, in contrast to differential correlation functions which require much more 
detailed analysis. 

Higher order perturbative solutions of the equations for the generating 
functions based on the running coupling constant, or the exact solution based 
on a fixed coupling, exhibit qualitative features which are absent in lowest 
order. Features predicted for partons, such as the oscillations of the Hq ra- 
tios or the larger factorial moments of quark jets compared to gluon jets, are 
experimentally observed for hadrons. The shape of the multiplicity distri- 
bution is mostly defined by soft particles created during the final stages of 
development of the cascade. The success of the theoretical methods implies 
the applicability of higher order perturbative techniques to the soft stages 
of evolution, and the suitability of local parton-hadron duality to multiplic- 
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ity based measurements. The principal features of multiplicity predicted by 
QCD are observed qualitatively both in e+e^ and hadron-initiated processes. 
This prompts speculation about the similarity of the production mechanisms 
in the two cases. This similarity probably occurs because the parton wave 
functions of the colliding particles in hadron collisions are prepared long be- 
fore the scattering, and are described by a space-like bremsstrahlung cascade 
with angular ordering similar to e^e^ case. 

The historical development of the perturbative approach to multiplicity is 
summarized in the introduction. The initial excitement created by the first 
applications of QCD to multiplicity gave way to disappointment over the 
extremely wide multiplicity distribution predicted by the theory at lowest 
order, although it was soon realized that the corrections were rather large. 
Now the importance of higher orders, up to NNLO and beyond, is well rec- 
ognized, and it is possible to account for them in a consistent manner. The 
results reveal approximate F (or KNO) scaling, with some dependence on 
the coupling constant in the preasymptotic regime. 

The long standing discrepancy between theory and experiment for the 
value of the ratio of the average multiplicities in gluon and quark jets has 
now been effectively resolved. Measurements of the higher moments of the 
multiplicity distribution are found to be in excellent agreement with the com- 
puter solutions of the QCD equations. The rise in average multiplicity with 
energy is described in a satisfactory manner, using a sensible value of the 
strong coupling strength as for both quark and gluon jets. The evolution 
of the distributions in smaller phase space regions has also been successfully 
described, at least at a qualitative level. Energy conservation in the calcula- 
tions and a proper unbiased definition of jets in the experiments are crucial 
for the quantitative tests of the QCD results, and are the principal develop- 
ments which allowed recent progress in the field. The astonishing success of 
the computer solutions also demonstrates the importance of using the cor- 
rect limits of integration over the parton energy splitting variables, i.e. the 
importance of the treatment of the boundary between the perturbative and 
non-perturbative regions. 

In combination with predictions for inclusive spectra and various corre- 
lation functions, the above results suggest that Quantum Chromodynamics 
can be successfully applied to predict qualitative features of soft processes. 
Angular ordering inside each jet and the collective nature of the inter-jet 
fiows are essential considerations for multiplicity and its energy evolution in 
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individual events. Analytic approaches provide a means to predict qualitative 
new effects and to gain insight into the behavior of the solutions (dependence 
on energy scale, the number of active quark flavors, etc.). Recent progress in 
that direction, described at some length above, gives hope for further success. 

17 Acknowledgments 

We are very much indebted to A. Capclla, Yu.L. Dokshitzer, G. Gianini, 
R.C. Hwa, C.S. Lam, B.B. Levtchenko, V.A. Nechitailo and J. Tran Thanh 
Van for collaboration and discussions. 

This work was supported by the Russian Fund for Basic Research, by IN- 
TAS, and by the US Department of Energy under grant DE-FG03-94ER40837. 



128 



References 

[1] I.M. Dremin, Physics-Uspekhi, 37 (1994) 715. 

[2] A Giovannini and L Van Hove, Z. Phys. C30 (1986) 391; Acta Phys. 
Pol. B19 (1988) 495; 917; 931. 

[3] I.V. Andreev, Chromodynamics and hard processes at high energies 
(Moscow, Nauka, 1981) (in Russian). 

[4] B.L. loffe, L.N. Lipatov and V.A. Khoze, Deep inelastic processes 
(Moscow, Energoatomizdat, 1983) (in Russian). 

[5] F.J. Yndurain, Quantum Chromodynamics (N.-Y.-Berlin- Heidelberg- 
Tokyo, Springer Verlag, 1983). 

[6] M.B. Voloshin and K.A. Ter-Martirosyan, Theory of gauge interactions 
of elementary particles (Moscow, Energoatomizdat, 1984) (in Russian). 

[7] Yu.L. Dokshitzer, V.A. Khoze, A.H. Mueller and S.I. Troyan, Basics of 
perturbative QCD (Gif-sur-Yvette, Editions Prontieres, 1991). 

[8] Ya.I. Azimov, Yu.L. Dokshitzer, V.A. Khoze and S.I. Troyan, Z. Phys. 
C27 (1985) 65. 

[9] D. Amati and G. Veneziano, Phys. Lett. B83 (1979) 87. 

[10] A. Bassctto, M. Ciafaloni and G. Marchesini, Phys. Rep. ClOO 
(1983) 201. 

[11] A. Bassctto, M. Ciafaloni and G. Marchesini, Nucl. Phys. B163 
(1980) 477. 

[12] Z. Koba, H.B. Nielsen and P. Olesen, Nucl. Phys. B40 (1972) 317. 
[13] Yu.L. Dokshitzer, Phys. Lett. B305 (1993) 295. 
[14] I.M. Dremin, Phys. Lett. B313 (1993) 209. 

[15] OPAL Collaboration, K. Ackerstaff et al., Eur. Phys. J. CI (1998) 479. 
[16] E.A. DeWolf, I.M. Dremin and W. Kittel, Phys. Rep. 270 (1996) 1. 

129 



[17] V.A. Khoze and W. Ochs, Int. J. Mod. Phys. A12 (1997) 2949. 

[18] S.K. Kaufmann and M. Gyulassi, J. Phys. All (1978) 1715. 

[19] S. Hegyi, Phys. Lett. B309 (1993) 443; B318 (1993) 642; B463 (1999) 
126. 

[20] E.M. Friedlander and I. Stern, LBL-31354, 1991. 

[21] K. Oldham, The fractional calculus (Orlando, Academic Press, 1974) 
p.60. 

[22] B. Ross, Fractional calculus and its applications, in Lecture Notes in 
Mathematics (Berlin, Springer Verlag, 1975) v.457, p.l. 

[23] I.M. Dremin, JETP Lett. 59 (1994) 585. 

[24] A.M. Polyakov, Sov. Phys. JETP 32 (1971) 296; 33 (1971) 850; 34 
(1972) 1177. 

[25] S. Hegyi, Proc. 8th Int. Conf. on Correlations and Fluctuations, Ma- 
trahaza, Hungary, 1998 (World Scientific, Singapore, to be published). 

[26] M. Ploszajczak, Proc. 8th Int. Conf. on Correlations and Fluctuations, 
Matrahaza, Hungary, 1998 (World Scientific, Singapore, to be pub- 
lished) . 

[27] R. Botet and M. Ploszajczak, Z. Phys. C76 (1997) 257. 
[28] S. Hegyi, Phys. Lett. B467 (1999) 126. 

[29] A.I. Golokhvastov, Sov. J. Nucl. Phys. 27 (1978) 430; 30 (1979) 253. 
[30] A. Wehrl, Rev. Mod. Phys. 50 (1978) 221. 

[31] V. Simak, M. Sumbera and I. Zborovsky, Proc. XVIII Int. Symp. on 
Multiparticle Dynamics, Tashkent, USSR, 1987, cds. I. Dremin and K. 
Gulamov (World Scientific, Singapore, 1988) p. 205. 

[32] V. Simak, M. Sumbera and I. Zborovsky, Phys. Lett. B206 (1988) 159. 



130 



[33] M. Pachr, V. Simak, M. Sumbera and I. Zborovsky, Mod. Phys. Lett. 
A7 (1992) 2333. 

[34] F. Takagi, Proc. XXI Int. Symp. on Multiparticle Dynamics, Wuhan, 
China, 1991, ed. C.S. Liu et al. (World Scientific, Singapore, 1992) 
p. 346. 

[35] I.M. Dremin, Mod. Phys. Lett. A8 (1993) 2747. 

[36] I.M. Dremin and R.C. Hwa, Phys. Rev. D49 (1994) 5805. 

[37] A. Giovannini, Nuovo Cim. A15 (1973) 543. 

[38] S. Hegyi, Phys. Lett. B387 (1996) 642; B417 (1998) 186. 

[39] M. Biyajima, Phys. Lett. B137 (1984) 225. 

[40] P.V. Chhapnikov, O.G. Tchikilev and V.A. Uvarov, Phys. Lett. B352 
(1995) 461. 

[41] V.D. Rusov, T.N. Zelenzova and S.J. Kosenko, OPU-98-1-ENIN, 1998. 

[42] S. Carius and G. Ingelmann, Phys. Lett. B252 (1990) 647. 

[43] R. Szwed and G. Wrochna, Z. Phys. C47 (1990) 447; 

R. Szwed, G. Wrochna and A.K. Wroblewski, Mod. Phys. Lett. A5 
(1990) 981. 

[44] I.M. Dremin and A.M. Dunaevskii, Phys. Rep. 18 (1975) 159. 

[45] B. Andersson, G. Gustafson and T. Sjostrand, Phys. Lett. B94 (1980) 
211. 

[46] B. Andersson, G. Gustafson, G. Ingelman and T. Sjostrand, Phys. Rep. 
97 (1983) 33. 

[47] A. Capclla, U. Sukhatme, C.I. Tan and J. Tran Thanh Van, Phys. Lett. 
B81 (1979) 68. 

[48] A. Capella and J. Tran Thanh Van, Phys. Lett. B93 (1980) 146; Z. 
Phys. CIO (1981) 249. 



131 



[49] A.B. Kaidalov, Phys. Lett. B116 (1982) 459. 

[50] A.B. Kaidalov and K.A. Ter-Martirosyan, Phys. Lett. B117 (1982) 247; 
Sov. J. Nucl. Phys. 39 (1984) 979; 40 (1984) 135. 

[51] R. Ugoccioni, A. Giovannini and S. Lupia, Proc. XXIII Int. Symp. on 
Muhiparticle Dynamics, Aspen, USA, 1993, eds. M.M. Block and A.R. 
White (World Scientific, Singapore, 1994) p. 297. 

[52] B.B. Levtchenko, Proc. VIII Workshop on High Energy Physics and 
Quantum Field Theory and III Workshop on Physics at VLEPP, 
Zvenigorod, Russia, 1993, ed. B.B. Levchenko (Moscow State Univ., 
Moscow, 1994) p. 68. 

[53] R. Ugoccioni and A. Giovannini, Nucl. Phys. Proc. Suppl. 71 (1999) 
201. 

[54] F. Bianchi, A. Giovannini, S. Lupia and R. Ugoccioni, Z. Phys. C58 

(1993) 71. 

[55] I.M. Dremin and D.S. Chernavsky, Zh. Exper. Teor. Fiz. 38 (1960) 229. 

[56] L. McLerran and R. Venugopalan, Phys. Rev. D49 (1994) 2233; D49 

(1994) 3352; D50 (1994) 2225. 

[57] Yu.V. Kovchegov, E.M. Levin and L. McLerran, BNL-NT-99/7, TAUP- 
2613-99. 

[58] E. Fermi, Progr. Theor. Phys. 5 (1950) 570. 

[59] I.Ya. Pomeranchuk, Doklady AN SSSR 78 (1951) 78. 

[60] L.D. Landau, Izvestia AN SSSR 17 (1953) 51. 

[61] E.A. Kuraev, J. Manjavidze, A.N. Sisakian, |hep-ph/9909287| ; |iep^ 



ph/0003094 



[62] UA5 Collaboration, G.J. Alner et al., Phys. Lett. B160 (1985) 193; 199. 
[63] A. Zichichi, M. Basile, G. Bonvicini et al, Proc. Erice 1982, p.701. 
[64] M. Basile et al, Lett. Nuovo Cim. 41 (1984) 293; 298. 



132 



[65] ALEPH Collaboration, D. Buskulic et al., Z. Phys. C69 (1995) 15. 
[66] OPAL Collaboration, G. Abbiendi et al, Eur. Phys. J. Cll (1999) 239. 
[67] E.K.G. Sarkisyan, Phys. Lett. B477 (2000) 1. 

[68] M. Biyajima, T. Kawabe and N. Suzuki, Phys. Lett. B189 (1987) 466. 

[69] T. Sjostrand, Comp. Phys. Comm. 82 (1994) 74; CERN-TH.7112/93 
(revised August 1995). 

[70] G. Marchesini, B.R. Webber et al., Comp. Phys. Comm. 67 (1992) 465. 

[71] L. L5nnblad, Comp. Phys. Comm. 71 (1992) 15. 

[72] B. Bambah et al., Z Physics at LEP, CERN 89-08, eds. G. AltareUi, R. 
Kleiss and C. Verzegnassi, (Geneva, 1989) vol.3, p. 143. 

[73] K. Goulianos, Phys. Lett. B193 (1987) 151. 

[74] X.N. Wang, M. Gyulassy and M. Pliimer, Phys. Rev. D51 (1995) 3436. 

[75] B.I. Ermolacv and V.S. Fadin, JETP Lett. 33 (1981) 269; 
A.H. Mueller, Phys. Lett. B104 (1981) 161. 

[76] K. Konishi, A. Ukawa and G. Veneziano, Nucl. Phys. B157 (1979) 45. 

[77] V.A. Khoze, S. Lupia and W. Ochs, Eur. Phys. J. C5 (1998) 77. 

[78] S. Lupia and W. Ochs, Phys. Lett. B418 (1998) 214; Nucl. Phys. Proc. 
Suppl. 64 (1998) 74. 

[79] W. A. Bardeen, A.J. Buras, D.W. Duke and T. Muta, Phys. Rev. D18 
(1978) 3998. 

[80] PDG, R.M. Barnett et al, Phys. Rev. D 54 (1996) 1. 

[81] Yu.L. Dokshitzer, V.A. Khoze and S.I. Troyan, Phys. Rev. D53 
(1996) 89. 



133 



[82] V.N.Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438; 
L.N. Lipatov, Sov. J. NucL Phys 20 (1975) 95; 
G. Altarelli and G. Parisi, NucL Phys. B126 (1977) 298; 
Yu.L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 

[83] Yu.L. Dokshitzer and M. Olsson, NucL Phys. B396 (1993) 137. 

[84] F. Cuypers and K. Tesima, Z. Phys. C54 (1992) 87. 

[85] I.M. Dremin, C.S. Lam, V.A. Nechitailo, Phys. Rev. D61 
(2000) 074020. 

[86] Yu.L. Dokshitzer, V.A. Khoze and S.I. Troyan, Perturbative QCD, ed. 
A.H. Mueller (World Scientific, Singapore, 1989). 

[87] E.D. Malaza and B.R. Webber, Nucl. Phys. B267 (1986) 702. 

[88] J.B. Gaffney and A.H. Mueller, Nucl. Phys. B250 (1985) 109. 

[89] F. Cuypers and K. Tesima, Z. Phys. C52 (1991) 69. 

[90] S. Catani, Yu.L. Dokshitzer, F. Fiorani and B.R. Webber, Nucl. Phys. 
B383 (1992) 419. 

[91] M. Olsson and G. Gustafson, Nucl. Phys. B406 (1993) 293. 

[92] S. Catani, Yu.L. Dokshitzer, F. Fiorani and B.R. Webber, Nucl. Phys. 
B377 (1992) 445. 

[93] I.M. Dremin and V.A. Nechitailo, Mod. Phys. Lett. A9 (1994) 1471. 

[94] Yu.L. Dokshitzer and I.M. Dremin, Nucl. Phys. B402 (1993) 139. 

[95] I.M. Dremin, Phys. Lett. B341 (1994) 95. 

[96] I.M. Dremin and V.A. Nechitailo, JETP Lett. 58 (1993) 881. 

[97] R. Ugoccioni, A. Giovannini and S. Lupia, Phys. Lett. B342 (1995) 387. 

[98] A. Capella, I.M. Dremin, J.W. Gary, V.A. Nechitailo and J. Tran 
Thanh Van, Phys. Rev. D61 (2000) 074009. 



134 



[99] I.M. Drcmin, B.B. Levtchenko and V.A. Nechitailo, Phys. Atom. Nucl. 
57 (1994) 1029. 

[100] I.M. Dremin and J.W. Gary, Phys. Lett. B459 (1999) 341. 

[101] K. Shizuya and S.-H. H. Tye, Phys. Rev. Lett. 41 (1978) 787. 

[102] M.B. Einhorn and B.G. Weeks, Nucl. Phys. B146 (1978) 445. 

[103] Ya.I. Azimov, Yu.L. Dokshitzer and V.A. Khoze, Sov. Phys. Uspekhi 
23 (1980) 732. 

[104] S. Lupia, Phys. Lett. B439 (1998) 150. 

[105] S.J. Brodsky and J.F. Gunion, Phys. Rev. Lett. 37 (1976) 402. 
[106] A.H. Mueller, Nucl. Phys. B241 (1984) 141. 
[107] P. Eden and G. Gustafson, JHEP 09 (1998) 015. 
[108] B.R. Webber, Phys. Lett. B143 (1984) 501. 

[109] Yu.L. Dokshitzer, V.A. Khoze and S.I. Troyan, Int. J. Mod. Phys. A7 
(1992) 1875. 

[110] W.J. Marciano, Phys. Rev. D29 (1984) 580. 

[Ill] V.I. Zakharov, Progr. Theor. Phys. Suppl. 131 (1998) 107. 

[112] OPAL Collaboration, G. Alexander et al., Phys. Lett. B388 (1996) 659. 

[113] OPAL Collaboration, G. Abbiendi et al, Eur. Phys. J. Cll (1999) 217. 

[114] CLEO Collaboration, M.S. Alam et al., Phys. Rev. D56 (1997) 17. 

[115] DELPHI Collaboration, P. Abreu et al., Z. Phys. C56 (1992) 63. 

[116] A. Giovannini, S. Lupia and R. Ugoccioni, Phys. Lett. B374 (1996) 231. 

[117] I.M. Dremin and R.C. Hwa, Phys. Lett. B324 (1994) 477. 

[118] Yu.L. Dokshitzer, V.A. Khoze, A.H. Mueller and S.I. Troyan, Rev. 
Mod. Phys. 60 (1988) 373. 

135 



[119] A. Giovannini and R. Ugoccioni, Proc. XXVIII Int. Symp. on Multi- 
particle Dynamics, Delphi, Greece, 1998 (World Scientific, Singapore, 
to be published). 

[120] A. Giovannini, S. Lupia and R. Ugoccioni, Phys. Lett. B388 (1996) 639. 
[121] J.W. Gary, unpublished. 

[122] I.M. Dremin, V.A. Nechitailo and M. Biyajima, J. Atom. Nucl. 59 

(1996) 2264. 

[123] I.M. Dremin, JETP Lett. 31 (1980) 185. 

[124] I.M. Dremin and A.V. Leonidov, Sov. J. Nucl. Phys. 35 (1982) 247. 

[125] A.V. Leonidov and D.M. Ostrovskii, Phys. Atom. Nucl. 60 (1997) 119. 

[126] A.V. Leonidov and D.M. Ostrovskii, Lebedev Physical Institute re- 
ports, 1-2 (1997) 36. 

[127] J. Ellis and K. Geiger, Phys. Rev. D52 (1995) 1500. 

[128] J. EUis and K. Geiger, Nucl. Phys. A590 (1995) 609c. 

[129] I.M. Dremin, JETP Lett. 60 (1994) 764. 

[130] R. Blankenbecler, Phys. Rev. D8 (1973) 1611. 

[131] C.N. Yang and T.D. Lee, Phys. Rev. 87 (1952) 404. 

[132] T.D. Lee and C.N. Yang, Phys. Rev. 87 (1952) 410. 

[133] E.A. De Wolf, Proc. XXIV Int. Symp. on Multiparticle Dynamics, Vi- 
etri sul Mare, Italy, 1994, eds. A. Giovannini, S. Lupia and R. Ugoccioni 
(World Scientific, Singapore, 1995) p. 15. 

[134] A. Capella, I.M. Dremin, V.A. Nechitailo and J. Tran Thanh Van, 
Z. Phys. C75 (1996) 89. 

[135] T.C. Brooks, K.L. Kowalski and C.C. Taylor, Phys. Rev. D56 

(1997) 5857. 



136 



[136] N. Brown and W.J. Stirling, Phys. Lett. B252 (1990) 657; Z. Phys. C53 
(1992) 629. 

[137] S. Catani, Yu.L. Dokshitzer, M. Olsson, G. Turnock and B.R. Webber, 
Phys. Lett. B269 (1991) 432. 

[138] JADE Collaboration, W. Bartcl et al, Z. Phys. C33 (1986) 23; 
S. Bethke et al., Phys. Lett. B213 (1988) 235. 

[139] Yu.L. Dokshitzer, V.A. Khoze and S.I. Troyan, Sov. J. Nucl. Phys. 47 

(1988) 881. 

[140] J.W. Gary, Phys. Rev. D49 (1994) 4503. 

[141] V.A. Khoze, Proc. Workshop on Physics with Linear Colliders, 
Saariselka, Finland, 1991, eds. R. Orava, E. Eerola and M. Nordberg 
(World Scientific, Singapore, 1992) p. 547. 

[142] JADE Collaboration, W. Bartel et al., Z. Phys. C21 (1983) 37. 

[143] UA2 Collaboration, P. Bagnaia et al, Phys. Lett. B144 (1984) 291. 

[144] TASSO Collaboration, M. Althoff et al., Z. Phys. C22 (1984) 307. 

[145] HRS CoUaboration, M. Derrick et al., Phys. Lett. B165 (1985) 449. 

[146] OPAL Collaboration, P.D. Acton et al, Z. Phys. C58 (1993) 387. 

[147] DELPHI Collaboration, P. Abreu et al, Z. Phys. C70 (1996) 179. 

[148] ALEPH CoUaboration, D. Buskulic et al., Phys. Lett. B384 (1996) 353. 

[149] OPAL Collaboration, R. Akers et al., Z. Phys. C68 (1995) 179. 

[150] CLEG Collaboration, M.S. Alam et al, Phys. Rev. D46 (1992) 4822. 

[151] ARGUS Collaboration, H. Albrecht et al., Z. Phys. C54 (1992) 13. 

[152] JADE Collaboration, W. Bartel et al., Z. Phys. C20 (1983) 187. 

[153] TASSO CoUaboration, W. Braunschweig et al., Z. Phys. C45 

(1989) 193. 



137 



HRS Collaboration, M. Derrick et al., Pliys. Rev. D34 (1986) 3304. 

AMY Collaboration, H.W. Zheng et al., Phys. Rev. D42 (1990) 737. 

ALEPH Collaboration, R. Barate et al., Phys. Rep. 294 (1998) 1; 
DELPHI Collaboration, P. Abreu et al., Eur. Phys. J. C6 (1999) 19; 
L3 Collaboration, B. Adeva et al., Z. Phys. C55 (1992) 39; 
OPAL Collaboration, K. Ackerstaff et al., Eur. Phys. J. C7 (1999) 369. 

ALEPH Collaboration Note 98-025, CONE 98-014. 

DELPHI Collaboration, P. Abreu et al., Phys. Lett. B372 (1996) 172; 
L3 Collaboration, M. Acciarri et al., Phys. Lett. B371 (1996) 137; 
OPAL Collaboration, G. Alexander et al., Z. Phys. C72 (1996) 191. 

DELPHI Collaboration, P. Abreu et al., Phys. Lett. B416 (1998) 233; 
L3 Collaboration, M. Acciarri et al., Phys. Lett. B404 (1997) 390. 

OPAL Collaboration, K. Ackerstaff et al., Z. Phys. C75 (1997) 193. 

OPAL Collaboration, G. Abbiendi et al., Eur. Phys. J. C16 (2000) 185. 

L3 Collaboration, M. Acciarri et al., Phys. Lett. B444 (1998) 569. 

ALEPH Collaboration Note 99-023, CONE 99-018. 

L3 Collaboration Note 2505. 

ALEPH Collaboration Note 2000-017, CONE 2000-014. 
OPAL Collaboration, Physics Note PN420. 

G. D. Lafferty, P. I. Reeves, and M. R. Whalley, Journal of Physics G, 
Nucl. Part. Phys. 21 (1995) A1-A151; 

The Durham RAL Databases, [http: / / durpdg.dur.ac.uk/HEPDATA . 



DELPHI Collaboration, P. Abreu et al., Phys. Lett. B449 (1999) 383. 

A.N. Safonov, Nucl. Phys. Proc. Suppl. 86 (2000) 55. 

OPAL Collaboration, G. Abbiendi et al., CERN-EP-2000-070. 



138 



[171] DELPHI Collaboration, CERN-OPEN-2000-134. 
[172] I.M. Dremin, JETP Lett. 68 (1998) 559. 

[173] S. Lupia, W. Ochs and J. Wosiek, Nucl. Phys. B540 (1999) 405. 

[174] SLD Collaboration, K. Abe et al., Phys. Lett. B371 (1996) 149. 

[175] L3 Collaboration, L3 Note 2407; W. Metzger, Proc. XXIX Int. Symp. 
on Multiparticle Dynamics, Providence, USA, 1999 (World Scientific, 
Singapore, to be published). 

[176] V.A. Khoze, S. Lupia and W. Ochs, Phys. Lett. B394 (1997) 179. 

[177] K. Konishi, A. Ukawa and G. Veneziano, Phys. Lett. B78 (1978) 243. 

[178] I.M. Dremin, V. Arena, G. Boca, G. Gianini et al., Phys. Lett. B336 
(1994) 119. 

[179] I.M. Dremin, V.A. Nechitailo, M. Biyajima and N. Suzuki, Phys. Lett. 
B403 (1997) 149. 

[180] R.K. Ellis, D.A. Ross and A.E. Terrano, Nucl. Phys. B178 (1981) 421. 

[181] G. Sterman and S. Weinberg, Phys. Rev. Lett. 39 (1977) 1436. 

[182] S. Bethke, Z. Kunszt, D.E. Soper and W.J. Stirhng, Nucl. Phys. B370 
(1992) 310. 

[183] L3 Collaboration, O. Adriani et al., Phys. Rep. 236 (1993) 1. 

[184] OPAL Collaboration, R. Akers et al., Z. Phys. C63 (1994) 363. 

[185] ALEPH Collaboration, D. Buskulic et al., Phys. Lett. B346 (1995) 389; 
ALEPH Collaboration, D. Barate et al., CERN-EP/98-016. 

[186] M.H. Seymour, Phys. Lett. B378 (1996) 279. 

[187] OPAL Collaboration, M.Z. Akrawy et al. Phys. Lett. B263 (1991) 311; 
OPAL Collaboration., G. Alexander et al, Phys. Lett. B265 (1991) 462. 

[188] J.W. Gary, Phys. Rev. D61 (2000) 114007. 



139 



[189] P. Eden, G. Gustafson and V. Khoze, Eur. Phys. J. Cll (1999) 345. 

[190] P. Eden and G. Gustafson, JHEP 9809 (1998) 015. 

[191] A. Bialas and R. Peschanski, Nucl. Phys. B273 (1986) 703. 

[192] I.M. Dremin, JETP Lett. 45 (1987) 643. 

[193] R.P. Fcynman, Proc. Ill Workshop on Current Problems in High En- 
ergy Particle Theory, Florence, Italy, 1979, eds. Casalbuoni et al. (John 
Hopkins University Press, Baltimore). 

[194] A. Giovannini, Nucl. Phys. B161 (1979) 429. 

[195] G. Veneziano, Proc. Ill Workshop on Current Problems in High Energy 
Particle Theory, Florence, Italy, 1979, eds. Casalbuoni et al. (John 
Hopkins University Press, Baltimore). 

[196] A. Bialas and R. Peschanski, Nucl. Phys. B308 (1988) 857. 

[197] W. Ochs and J. Wosiek, Phys. Lett. B214 (1988) 617. 

[198] I.M. Dremin, Mod. Phys. Lett. A4 (1989) 1027. 

[199] C.B. Chiu and R. Hwa, Phys. Lett. B236 (1990) 446. 

[200] G. Gustafson and A. Nilsson, Z. Phys. C52 (1991) 533; Nucl. Phys. 
B355 (1991) 106. 

[201] W. Ochs and J. Wosiek, Phys. Lett. B289 (1992) 159; B304 (1993) 144. 

[202] W. Ochs and J. Wosiek, Z. Phys. C68 (1995) 269. 

[203] Ph. Brax, J.L. Meunier and R. Peschanski, Z. Phys. C62 (1994) 649. 

[204] Yu.L. Dokshitzer, G. Marchesini and G. Oriani, Nucl. Phys. B387 
(1992) 675. 

[205] J.L. Meunier and R. Peschanski, Z. Phys. C72 (1996) 647. 

[206] I.V. Andrecv, M. Biyajima, I.M. Dremin and N. Suzuki, Int. J. Mod. 
Phys. AlO (1995) 3951. 



140 



L3 Collaboration, M. Acciarri et al., Phys. Lett. B428 (1998) 186. 

DELPHI Collaboration, P. Abreu et al, Phys. Lett. B457 (1999) 368. 

L3 Collaboration, M. Acciarri et al., Phys. Lett. B429 (1998) 375. 

S.V. Chekanov and V.I. Kuvshinov, Acta Phys. Polonica B25 
(1994) 1189. 

S.V. Chekanov, W. Kittel and V.I. Kuvshinov, Z. Phys. C74 (1997) 517. 
M.I. Polikarpov, Phys. Lett. B236 (1990) 61. 
I.M. Dremin and B.B. Levtchenko, Phys. Lett. B292 (1992) 155. 
I.M. Dremin, JETP Lett. 34 (1981) 594. 

I.M. Dremin, JETP Lett. 30 (1979) 140; Sov. J. Nucl. Phys. 33 
(1981) 726. 

A. V. Apanasenko, N.A. Dobrotin, I.M. Dremin and K.A. Kotelnikov, 
JETP Lett. 30 (1979) 145. 

B. A. Schumm, Yu.L. Dokshitzer, V.A. Khoze and D.S. Koetke, Phys. 
Rev. Lett. 69 (1992) 3025. 

Yu.L. Dokshitzer, V.A. Khoze and S.I. Troyan, J. Phys. G17 (1991) 
1481; 1602. 

Ya.I. Azimov, Yu.L. Dokshitzer and V.A. Khoze, Sov. J. Nucl. Phys. 
36 (1982) 878. 

A.V. Kisselev, V.A. Petrov and O.P. Yuschenko, Z. Phys. C41 
(1988) 521. 

V.A. Petrov and A.V. Kisselev, Z. Phys. C66 (1995) 453. 

DELPHI Collaboration, P. Abreu et al, Phys. Lett. B479 (2000) 118. 

I.M. Dremin and V.I. Yakovlev, Proc. XVII Int. Symp. on Multiparticle 
Dynamics, Austria, 1986, ed. M. Markityan (World Scientific, Singa- 
pore, 1987) p. 849. 



141 



